Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 224,933 A→G 44.3% noncoding (1163/1542 nt) rrsH → 16S ribosomal RNA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_007779224,9330A→G44.3% 0.0 ‑3.3 ‑0.0272noncoding (1163/1542 nt)rrsH16S ribosomal RNA
Reads supporting (aligned to +/- strand):  new base (77/36):  ref base (76/62):  total (153/119)
Fisher's exact test strand distribution p-value = 3.81e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.90e-01

Only 100 of 279 total aligned reads displayed.

CGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAG  >  NC_007779/224833‑224990
                                                                                                        |                                                          
 gtcAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAg                                                           >  1:10468936/1‑101 (MQ=32)
 gtcAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAg                                                           <  1:8367958/101‑1 (MQ=32)
 gtcAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAg                                                           >  1:4764655/1‑101 (MQ=32)
  tcAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGt                                                          >  1:14311867/1‑101 (MQ=32)
   cAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTg                                                         >  1:11046264/1‑101 (MQ=32)
    aGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGa                                                        <  1:8591216/101‑1 (MQ=32)
    aGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGa                                                        <  1:5475374/101‑1 (MQ=32)
     gCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGAt                                                       >  1:2194022/1‑101 (MQ=32)
       tCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATaa                                                     <  1:13712898/101‑1 (MQ=32)
         gtgtTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAAc                                                   >  1:7832544/1‑101 (MQ=32)
           gttgtGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTg                                                <  1:12860722/101‑1 (MQ=32)
             tgtgAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTgga                                              >  1:969084/1‑101 (MQ=33)
             tgtgAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTgga                                              >  1:7588289/1‑101 (MQ=33)
              gtgAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggag                                             <  1:14591708/101‑1 (MQ=32)
               tgAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagg                                            >  1:10869587/1‑101 (MQ=32)
                gAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                           >  1:2950821/1‑101 (MQ=32)
                gAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                           <  1:5607736/101‑1 (MQ=32)
                 aaaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGaa                                          >  1:9261956/1‑101 (MQ=33)
                 aaaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGaa                                          >  1:7756259/1‑101 (MQ=33)
                 aaaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGaa                                          <  1:1489724/101‑1 (MQ=32)
                  aaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAg                                         <  1:1238811/101‑1 (MQ=32)
                   agtttGGGTAAAGCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGt                                       <  1:5914546/98‑1 (MQ=255)
                   aTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAgg                                        <  1:1862033/101‑1 (MQ=32)
                        gggTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTgggg                                   >  1:11550575/1‑101 (MQ=33)
                           ttAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGAt                                 >  1:3810199/1‑101 (MQ=2)
                           ttAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGAt                                 >  1:11667930/1‑101 (MQ=2)
                           ttAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGAt                                 >  1:13896396/1‑101 (MQ=2)
                           ttAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATg                                <  1:9576679/101‑1 (MQ=32)
                            tAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATg                                >  1:6203570/1‑101 (MQ=2)
                            tAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATg                                >  1:12232247/1‑101 (MQ=2)
                            tAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATg                                >  1:6746680/1‑101 (MQ=2)
                            tAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATg                                >  1:1273910/1‑101 (MQ=2)
                            tAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGa                               <  1:5570598/101‑1 (MQ=32)
                             aaGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGa                               >  1:12638837/1‑101 (MQ=2)
                               gTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACg                             >  1:3627414/1‑101 (MQ=2)
                               gTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACg                             >  1:4496971/1‑101 (MQ=2)
                               gTCCCCCAGCGCGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGt                            <  1:7894242/101‑1 (MQ=255)
                                tCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGt                            <  1:7565641/101‑1 (MQ=2)
                                tCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGt                            >  1:14802712/1‑101 (MQ=2)
                                tCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTc                           >  1:1411362/1‑101 (MQ=32)
                                 cccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTc                           >  1:12474948/1‑101 (MQ=2)
                                 cccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGGGGGGACCCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGAGGTc                           >  1:5373648/1‑101 (MQ=2)
                                 cccGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCa                          <  1:283132/101‑1 (MQ=32)
                                  ccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCa                          <  1:7112975/101‑1 (MQ=2)
                                  ccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCa                          <  1:12991142/101‑1 (MQ=2)
                                  ccGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCaa                         >  1:4255874/1‑101 (MQ=33)
                                   cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         <  1:11880653/101‑1 (MQ=2)
                                   cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         >  1:4870823/1‑101 (MQ=2)
                                   cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         <  1:8019631/101‑1 (MQ=2)
                                   cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         <  1:8839515/101‑1 (MQ=2)
                                   cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         <  1:10663007/101‑1 (MQ=2)
                                    gCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaa                        <  1:11233666/101‑2 (MQ=255)
                                    gCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaa                        >  1:14788583/1‑100 (MQ=2)
                                    gCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaa                        <  1:305580/101‑2 (MQ=2)
                                    gCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGt                       >  1:332245/1‑101 (MQ=33)
                                     cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       <  1:4641972/101‑3 (MQ=2)
                                     cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       <  1:9922539/101‑3 (MQ=2)
                                     cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       <  1:6976703/101‑3 (MQ=255)
                                     cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       <  1:11262532/101‑3 (MQ=255)
                                     cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGGGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       >  1:1673896/1‑99 (MQ=2)
                                      aaCGAGCGCAACCCTTGTCTTATGTTGCCAGCACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:15187159/1‑101 (MQ=255)
                                      aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      <  1:4025118/101‑1 (MQ=2)
                                      aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:266779/1‑101 (MQ=2)
                                      aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:14662449/1‑101 (MQ=255)
                                      aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:14857924/1‑101 (MQ=255)
                                      aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:5764163/1‑101 (MQ=2)
                                      aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:3431069/1‑101 (MQ=2)
                                       aCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtca                     >  1:2352245/1‑101 (MQ=255)
                                       aCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtca                     <  1:10540891/101‑1 (MQ=255)
                                       aCGAGCACAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtca                     >  1:10777027/1‑101 (MQ=255)
                                        cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:6136120/1‑101 (MQ=2)
                                        cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:5085808/1‑101 (MQ=2)
                                        cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    <  1:9465607/101‑1 (MQ=2)
                                        cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    <  1:4142908/101‑1 (MQ=2)
                                        cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:11343616/1‑101 (MQ=2)
                                        cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:7659766/1‑101 (MQ=2)
                                        cGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatc                   <  1:12059981/101‑1 (MQ=32)
                                         gAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatc                   >  1:12000099/1‑101 (MQ=2)
                                         gAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                  >  1:8177441/1‑101 (MQ=33)
                                          aGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatca                  <  1:3336296/101‑1 (MQ=2)
                                          aGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatca                  >  1:1345884/1‑101 (MQ=2)
                                          aGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatca                  >  1:13062136/1‑101 (MQ=2)
                                          aGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatca                  <  1:13567528/101‑1 (MQ=2)
                                           gcgcAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATcatcat                 <  1:1931521/101‑1 (MQ=2)
                                           gcgcAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATcatcat                 <  1:10668284/101‑1 (MQ=2)
                                           gcgcAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATcatcat                 <  1:14260874/101‑1 (MQ=2)
                                           gcgcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCCAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATg                >  1:8712559/1‑101 (MQ=17)
                                            cgcAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                <  1:3688505/101‑1 (MQ=2)
                                            cgcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg               <  1:14461225/101‑1 (MQ=32)
                                            cgcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg               <  1:2591050/101‑1 (MQ=32)
                                            cgcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg               <  1:10342023/101‑1 (MQ=32)
                                             gcAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc               <  1:13216332/101‑2 (MQ=2)
                                             gcAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc               >  1:6422932/1‑100 (MQ=2)
                                                 cccTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCtt          <  1:9798830/101‑1 (MQ=32)
                                                  ccTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTa         >  1:2088519/1‑101 (MQ=32)
                                                   cTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc        >  1:8677855/1‑101 (MQ=32)
                                                   cTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc        <  1:963755/101‑1 (MQ=32)
                                                     tATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGa      <  1:10105404/101‑1 (MQ=32)
                                                      aTCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAc     >  1:6539879/1‑101 (MQ=32)
                                                         ttttGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAg  <  1:6701383/101‑1 (MQ=11)
                                                                                                        |                                                          
CGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAG  >  NC_007779/224833‑224990

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.