Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 224,937 A→T 39.5% noncoding (1167/1542 nt) rrsH → 16S ribosomal RNA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_007779224,9370A→T39.5% 260.9 ‑5.3 0.0231noncoding (1167/1542 nt)rrsH16S ribosomal RNA
Reads supporting (aligned to +/- strand):  new base (65/25):  ref base (69/62):  total (134/97)
Fisher's exact test strand distribution p-value = 4.92e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

AGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAG  >  NC_007779/224837‑224990
                                                                                                        |                                                      
aGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGa                                                        >  1:4169562/1‑101 (MQ=32)
aGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGa                                                        >  1:1536829/1‑101 (MQ=32)
 gCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGAt                                                       >  1:4018841/1‑101 (MQ=32)
  cTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATa                                                      <  1:3456723/101‑1 (MQ=32)
  cTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATa                                                      <  1:11373322/101‑1 (MQ=32)
  cTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATa                                                      >  1:6995967/1‑101 (MQ=32)
    cGTGTTGTGAAATGTTGGGTTAGGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATaaa                                                    >  1:4144662/1‑101 (MQ=17)
    cGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATaaa                                                    >  1:13938631/1‑101 (MQ=32)
    cGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATaaa                                                    <  1:3025758/101‑1 (MQ=32)
      tgttgtGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACt                                                  <  1:14967489/101‑1 (MQ=32)
      tgttgtGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACt                                                  <  1:12295779/101‑1 (MQ=32)
      tgttgtGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACt                                                  <  1:8723611/101‑1 (MQ=32)
       gttgtGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTg                                                <  1:2695056/101‑1 (MQ=32)
       gttgtGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTg                                                >  1:14083849/1‑101 (MQ=32)
        ttgtGAAATGTTGGGTTAAGTCCCGCAACGAGAGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGCACTCAAAGGAGACTGCCAGTGATAAACTgg                                               >  1:6062233/1‑101 (MQ=12)
            gAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                           <  1:11808385/101‑1 (MQ=32)
             aaaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGaa                                          >  1:7118751/1‑101 (MQ=32)
                 gTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTg                                      >  1:1368628/1‑101 (MQ=33)
                   tGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTggg                                    >  1:13215874/1‑101 (MQ=33)
                   tGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTggg                                    <  1:4448702/101‑1 (MQ=32)
                     ggTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGa                                  <  1:7011487/101‑1 (MQ=32)
                      gTTAAGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGa                                  <  1:1967432/101‑1 (MQ=2)
                      gTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGAt                                 >  1:4335639/1‑101 (MQ=33)
                       ttAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATg                                >  1:9541269/1‑101 (MQ=32)
                         aaGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGa                               <  1:8408253/101‑1 (MQ=2)
                         aaGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGa                               >  1:7139193/1‑101 (MQ=2)
                         aaGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGACCTCAAAGGGGACTGCCAGGGATAAACTGGAGAAAGGTGGGGATGag                              >  1:12384642/1‑100 (MQ=11)
                         aaGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGAc                              >  1:6557497/1‑101 (MQ=33)
                          aGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGAc                              >  1:1416745/1‑101 (MQ=2)
                          aGTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGAc                              >  1:13255388/1‑101 (MQ=2)
                           gTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACg                             >  1:6322633/1‑101 (MQ=2)
                           gTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACg                             >  1:1870943/1‑101 (MQ=2)
                           gTCCCGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACg                             >  1:2923112/1‑101 (MQ=2)
                            tCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTc                           >  1:14829211/1‑101 (MQ=33)
                            tCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTc                           <  1:8155690/101‑1 (MQ=32)
                             cccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTc                           >  1:3369832/1‑101 (MQ=2)
                             cccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTc                           >  1:1705614/1‑101 (MQ=2)
                             cccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTc                           >  1:10267797/1‑101 (MQ=2)
                             cccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTc                           >  1:1757437/1‑101 (MQ=2)
                             cccGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCa                          >  1:8508340/1‑101 (MQ=33)
                              ccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCa                          <  1:9812015/101‑1 (MQ=2)
                              ccGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCa                          >  1:12621661/1‑101 (MQ=2)
                              ccGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCaa                         >  1:7525210/1‑101 (MQ=33)
                               cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         >  1:1665840/1‑101 (MQ=2)
                               cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         <  1:7573479/101‑1 (MQ=2)
                               cGCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaa                         >  1:11689713/1‑101 (MQ=2)
                               cGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACACAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAg                        >  1:5028699/1‑101 (MQ=18)
                                gCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGTGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaa                        >  1:14625844/1‑100 (MQ=2)
                                gCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaa                        >  1:14425244/1‑100 (MQ=2)
                                gCAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaa                        >  1:6372330/1‑100 (MQ=2)
                                gCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGt                       >  1:7626131/1‑101 (MQ=33)
                                gCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGt                       <  1:3938455/101‑1 (MQ=32)
                                gCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGt                       <  1:730919/101‑1 (MQ=32)
                                 cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       >  1:12800556/1‑99 (MQ=2)
                                 cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       <  1:10823720/101‑3 (MQ=2)
                                 cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       >  1:11594644/1‑99 (MQ=2)
                                 cAACGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCaaat                       <  1:3534264/101‑3 (MQ=2)
                                 cAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtc                      >  1:3339860/1‑101 (MQ=33)
                                  aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:7722506/1‑101 (MQ=2)
                                  aaCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtc                      >  1:3207203/1‑101 (MQ=255)
                                  aaCGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtca                     >  1:312224/1‑101 (MQ=33)
                                  aaCGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtca                     <  1:13570961/101‑1 (MQ=32)
                                  aaCGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtca                     <  1:13123209/101‑1 (MQ=32)
                                   aCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGAGGTCAAAtca                     <  1:2995410/101‑1 (MQ=2)
                                   aCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtca                     >  1:10537321/1‑101 (MQ=255)
                                   aCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtca                     <  1:5689020/101‑1 (MQ=2)
                                   aCGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtca                     >  1:11824766/1‑101 (MQ=2)
                                   aCGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcat                    >  1:3313411/1‑101 (MQ=33)
                                   aCGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcat                    >  1:7896645/1‑101 (MQ=33)
                                    cGAGCGCAACCCTTGTCTTATGTTGCCAGCACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:9488718/1‑101 (MQ=2)
                                    cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:2100405/1‑101 (MQ=2)
                                    cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:5490313/1‑101 (MQ=2)
                                    cGAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    >  1:9699476/1‑101 (MQ=2)
                                    cGAGCGCAACCCTTGTATTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcat                    <  1:9390749/101‑1 (MQ=2)
                                    cGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatc                   <  1:5644108/101‑1 (MQ=32)
                                     gAGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatc                   >  1:11260155/1‑101 (MQ=2)
                                     gAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                  >  1:2611763/1‑101 (MQ=33)
                                     gAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                  >  1:1514333/1‑101 (MQ=33)
                                     gAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                  >  1:3792444/1‑101 (MQ=33)
                                      aGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatca                  >  1:8021381/1‑101 (MQ=2)
                                      aGCGCAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAAtcatca                  <  1:528770/101‑1 (MQ=2)
                                      aGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTcatcat                 <  1:8814654/101‑1 (MQ=32)
                                       gcgcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATg                <  1:11187351/101‑1 (MQ=32)
                                        cgcATCCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                <  1:9091600/101‑1 (MQ=2)
                                        cgcAACCCTTGTCTTATGTTGCCAG‑‑CACGTTGTGGTGGGGACTCGTAAGAAACTGCCGGGGTTAACTCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                >  1:8339548/1‑101 (MQ=2)
                                        cgcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg               <  1:10668324/101‑1 (MQ=32)
                                         gcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGc              >  1:13742166/1‑101 (MQ=32)
                                         gcAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGc              >  1:10762940/1‑101 (MQ=32)
                                          cAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGcc             >  1:8346063/1‑101 (MQ=32)
                                           aaCCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGccc            >  1:7675972/1‑101 (MQ=32)
                                             cccTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCtt          >  1:4870501/1‑101 (MQ=32)
                                             cccTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCtt          <  1:1478896/101‑1 (MQ=32)
                                             cccTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGGGGCGATGACGTCAAGTCATCATGGCCCtt          >  1:12095070/1‑101 (MQ=12)
                                              ccTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTa         <  1:11877573/101‑1 (MQ=32)
                                               cTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc        >  1:2144383/1‑101 (MQ=32)
                                               cTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc        <  1:4601731/101‑1 (MQ=32)
                                                ttATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACg       <  1:5884058/101‑1 (MQ=32)
                                                  aTCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAc     <  1:5550027/101‑1 (MQ=32)
                                                  aTCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAc     >  1:2299803/1‑101 (MQ=32)
                                                     ttttGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAg  >  1:4774339/1‑101 (MQ=11)
                                                                                                        |                                                      
AGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAG‑‑CGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAG  >  NC_007779/224837‑224990

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ 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.