Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,402 G→A 18.1% R274H (CGC→CAC)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,441,4020G→A18.1% 138.8 ‑0.4 ‑0.0287R274H (CGC→CAC) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (27/25):  ref base (113/122):  total (140/147)
Fisher's exact test strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.91e-01

Only 100 of 292 total aligned reads displayed.

GGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGC  >  NC_007779/3441302‑3441501
                                                                                                    |                                                                                                   
ggCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCg                                                                                                     <  1:899729/101‑1 (MQ=255)
 gCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGGTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGc                                                                                                    >  1:12101124/1‑101 (MQ=255)
  cAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCt                                                                                                   <  1:2803231/101‑1 (MQ=255)
    aCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtat                                                                                                 <  1:13573368/101‑1 (MQ=255)
    aCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtat                                                                                                 <  1:10776108/101‑1 (MQ=255)
    aCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtat                                                                                                 <  1:7420000/101‑1 (MQ=255)
     cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                                                                <  1:3565201/101‑1 (MQ=255)
     cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                                                                <  1:7537819/101‑1 (MQ=255)
        ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                                             >  1:12445359/1‑101 (MQ=255)
        ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                                             >  1:10881557/1‑101 (MQ=255)
        ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                                             >  1:7094417/1‑101 (MQ=255)
         gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                                                                            >  1:11450826/1‑101 (MQ=255)
          cggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCAc                                                                                           <  1:5006198/101‑1 (MQ=255)
              cACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGaaa                                                                                       <  1:3636979/101‑1 (MQ=255)
               aCCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAc                                                                                      <  1:10342766/101‑1 (MQ=255)
               aCCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAc                                                                                      >  1:10399248/1‑101 (MQ=255)
                 ctctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAccc                                                                                    <  1:1155442/101‑1 (MQ=255)
                 ctctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAccc                                                                                    <  1:11201106/101‑1 (MQ=255)
                  tctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcg                                                                                   >  1:922457/1‑101 (MQ=255)
                   ctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgc                                                                                  >  1:1021243/1‑101 (MQ=255)
                    ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                                                                 >  1:6169851/1‑101 (MQ=255)
                       catctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAgg                                                                              >  1:13738428/8‑100 (MQ=255)
                       cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                                                              <  1:13679159/101‑1 (MQ=255)
                         tctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTgg                                                                            >  1:11105855/6‑101 (MQ=255)
                          cGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGa                                                                           <  1:6616350/101‑1 (MQ=255)
                            ccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAtt                                                                         >  1:3880483/3‑101 (MQ=255)
                            aagaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtt                                                                         >  1:3586474/1‑101 (MQ=255)
                            aagaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtt                                                                         >  1:12882687/1‑101 (MQ=255)
                             agaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAttt                                                                        <  1:4346110/101‑1 (MQ=255)
                              gaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtttt                                                                       >  1:8907280/1‑101 (MQ=255)
                               aagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTg                                                                      >  1:4635824/1‑101 (MQ=255)
                                agaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGc                                                                     >  1:161751/1‑101 (MQ=255)
                                 gaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCg                                                                    <  1:1219757/98‑1 (MQ=255)
                                 gaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCg                                                                    >  1:2434578/4‑101 (MQ=255)
                                  aggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGt                                                                   <  1:5116628/99‑1 (MQ=255)
                                   ggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTa                                                                  <  1:8894743/100‑1 (MQ=255)
                                   ggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTa                                                                  >  1:1998057/2‑101 (MQ=255)
                                      ggAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACa                                                               <  1:3766939/101‑1 (MQ=255)
                                       gAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAg                                                              <  1:9874337/101‑1 (MQ=255)
                                        aaaaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGc                                                             >  1:10003788/1‑101 (MQ=255)
                                          aaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGccc                                                           >  1:6581786/1‑101 (MQ=255)
                                           agatGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCg                                                          <  1:4452833/98‑1 (MQ=255)
                                           aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg                                                          >  1:5979394/1‑101 (MQ=255)
                                           aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg                                                          >  1:4187870/1‑101 (MQ=255)
                                            gatGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCgg                                                         <  1:7271961/99‑1 (MQ=255)
                                             atGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGt                                                        >  1:8063158/2‑101 (MQ=255)
                                              tGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTg                                                       >  1:12028769/1‑101 (MQ=255)
                                               ggTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGc                                                      <  1:8548026/101‑1 (MQ=255)
                                               ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGc                                                      >  1:5778246/1‑101 (MQ=255)
                                               ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGATCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGc                                                      <  1:10294207/101‑1 (MQ=255)
                                                 tATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg                                                    >  1:1408884/1‑101 (MQ=255)
                                                     gTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCg                                                >  1:6448514/1‑101 (MQ=255)
                                                        gacccgttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTa                                             <  1:7699529/94‑1 (MQ=255)
                                                        gacccgttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTa                                             <  1:5298375/94‑1 (MQ=255)
                                                        gacccgttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTa                                             >  1:90806/8‑101 (MQ=255)
                                                         acccgttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTAc                                            <  1:10500374/95‑1 (MQ=255)
                                                         aCGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTAc                                            >  1:1265006/1‑101 (MQ=255)
                                                           gcgcTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAg                                          >  1:12774242/1‑101 (MQ=255)
                                                           gcgcTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAg                                          <  1:10749780/101‑1 (MQ=255)
                                                             gttGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTACGGTg                                        >  1:2676414/3‑101 (MQ=255)
                                                             gcTNGGGAGCGGATACGGTGATGGATCTCTCCACCGGTNGCNATATTCACGAAACCCGCGAGTGGATNTTGCGTAACAGCCCGNTGCNGANCGGTACAGTg                                        >  1:9500783/1‑101 (MQ=255)
                                                             gcTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTg                                        >  1:1772293/1‑101 (MQ=255)
                                                              cTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGc                                       >  1:762147/1‑101 (MQ=255)
                                                               tGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGcc                                      >  1:475710/1‑101 (MQ=255)
                                                                ggggCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTACGGTGCCg                                     <  1:789337/101‑1 (MQ=255)
                                                                ggggAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCg                                     >  1:5644541/1‑101 (MQ=255)
                                                                 gggAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGa                                    <  1:6331395/101‑1 (MQ=255)
                                                                     gCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTa                                <  1:6847453/101‑1 (MQ=255)
                                                                      cGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAc                               <  1:13085594/101‑1 (MQ=255)
                                                                       ggATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAcc                              <  1:2154213/101‑1 (MQ=255)
                                                                       ggATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTAcc                              >  1:5218129/1‑101 (MQ=255)
                                                                        gATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCa                             <  1:6984204/101‑1 (MQ=255)
                                                                        gATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCa                             <  1:8632601/101‑1 (MQ=255)
                                                                        gATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCa                             <  1:13996982/101‑1 (MQ=255)
                                                                          tACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAgg                           <  1:11778545/101‑1 (MQ=255)
                                                                          tACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAgg                           <  1:10714373/101‑1 (MQ=255)
                                                                           aCGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGgc                          <  1:5164988/101‑1 (MQ=255)
                                                                            cGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGgcg                         <  1:5872438/101‑1 (MQ=255)
                                                                             catcaTGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTACGGTGCCGATTTACCAGgcgc                        <  1:11625115/97‑1 (MQ=255)
                                                                                  tGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGgag                   >  1:11376632/1‑101 (MQ=255)
                                                                                    gATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGaa                 <  1:8353938/101‑1 (MQ=255)
                                                                                     aTCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAg                >  1:10831347/1‑101 (MQ=255)
                                                                                     aTCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAg                >  1:13397344/1‑101 (MQ=255)
                                                                                      tctctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAgg               >  1:7832366/1‑101 (MQ=255)
                                                                                      tctctcCACCGGGCGCTATATTCACGAAACCCGCGAGTGGATTTTTCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAgg               >  1:970796/1‑101 (MQ=255)
                                                                                       ctctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGt              >  1:10610233/1‑101 (MQ=255)
                                                                                       ctctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGt              <  1:8713569/101‑1 (MQ=255)
                                                                                         ctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTa            <  1:1546244/101‑1 (MQ=255)
                                                                                         ctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTa            <  1:13834055/101‑1 (MQ=255)
                                                                                         ctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTa            <  1:3020972/101‑1 (MQ=255)
                                                                                         ctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTa            >  1:2719132/1‑101 (MQ=255)
                                                                                         ctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTa            >  1:1488879/1‑101 (MQ=255)
                                                                                          tcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTaa           <  1:9284186/101‑1 (MQ=255)
                                                                                            cACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACg         <  1:6902246/101‑1 (MQ=255)
                                                                                            cACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACg         <  1:13525885/101‑1 (MQ=255)
                                                                                              ccGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACggg       >  1:3352177/1‑101 (MQ=255)
                                                                                                 gTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATc    <  1:3948246/101‑1 (MQ=255)
                                                                                                  tCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCg   >  1:1648759/1‑101 (MQ=255)
                                                                                                   cGCTATATTCACGAAACCCGCGAGTGGATTTTTCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGc  >  1:1698901/1‑101 (MQ=255)
                                                                                                   cGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGc  >  1:12248314/1‑101 (MQ=255)
                                                                                                    |                                                                                                   
GGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGC  >  NC_007779/3441302‑3441501

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.