Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,327,121 C→T 51.7% intergenic (+97/‑196) PP_16SE → / → PP_23SE 16S ribosomal RNA/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,327,1210CT51.7% 27.8 / 140.6 107intergenic (+97/‑196)PP_16SE/PP_23SE16S ribosomal RNA/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (29/22);  new base T (32/22);  total (61/46)
Fisher's exact test for biased strand distribution p-value = 8.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGT  >  NC_002947/1326977‑1327227
                                                                                                                                                    |                                                                                                          
tCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATc                                                                                                         >  2:247673/1‑148 (MQ=17)
tCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATc                                                                                                         >  2:165537/1‑148 (MQ=32)
  gTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGaa                                                                                                      >  5:275446/1‑149 (MQ=17)
  gTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGaa                                                                                                      >  2:203993/1‑149 (MQ=17)
    aaCAAGGTAGCCGTAGGGGAACCTGCCGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCAGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGATTTTCCTGGTAAGAACCTGGAAATGAGCATTCGCATCGAATg                                                                                                    >  1:23377/1‑149 (MQ=255)
     aCAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGATGCGAATgc                                                                                                   <  3:164791/149‑2 (MQ=12)
      cAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGtt                                                                                                  <  6:301631/149‑1 (MQ=17)
      cAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGtt                                                                                                  <  1:203993/149‑1 (MQ=17)
       aaGGTAGCCGTAGGGGAACCTGCGGTTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGtt                                                                                                  >  4:144520/1‑148 (MQ=17)
       aaGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGtt                                                                                                  >  4:282993/1‑148 (MQ=17)
          gTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTAGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGAt                                                                                               >  1:263873/1‑149 (MQ=12)
          gTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGAtt                                                                                              >  6:44775/1‑149 (MQ=17)
          gTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGAtt                                                                                              >  1:257614/1‑149 (MQ=32)
           tatccGTAGGGGAACCTGCGGCTTGTTCACCTCCTTAATCGTCGACATTATCCCCCTGCTGATACCTTCTCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGa                                                                                                <  1:280094/146‑1 (MQ=11)
             gCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTc                                                                                            <  6:46224/148‑1 (MQ=32)
              ccGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTg                                                                                          <  1:165537/149‑1 (MQ=32)
               cGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTgg                                                                                         <  8:181839/149‑1 (MQ=17)
               cGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATCAACATTCGATACGAATGCTCATTTccag                                                                                         <  8:761/149‑4 (MQ=11)
                gTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGc                                                                                        >  7:192014/1‑149 (MQ=32)
                    gggAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCtttt                                                                                    <  5:127338/149‑1 (MQ=32)
                      gAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGt                                                                                  <  6:275446/149‑1 (MQ=17)
                       aaCCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTc                                                                                 <  2:40369/149‑1 (MQ=17)
                       aaCCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGt                                                                                  >  2:58104/1‑148 (MQ=32)
                         ccTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGt                                                                                  >  4:160746/1‑146 (MQ=32)
                         ccTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGt                                                                                  >  2:249667/1‑146 (MQ=32)
                               gCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGtt                                                                         <  1:247673/149‑1 (MQ=17)
                                cttgATCACCTCCTTAATCGCCGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTc                                                                        <  8:286785/146‑1 (MQ=17)
                                  ggATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCtt                                                                      >  7:249816/1‑149 (MQ=17)
                                      cACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaa                                                                   >  6:304488/1‑148 (MQ=17)
                                       aCCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaaaa                                                                 <  1:249667/149‑1 (MQ=32)
                                         ctcctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAAt                                                               >  1:162156/1‑149 (MQ=17)
                                         ctcctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTAGGGGCTTTTGTCAGATCGTTATTTAAAAAt                                                               >  8:51504/1‑149 (MQ=255)
                                         ctcctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATAATCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAAt                                                               >  1:267074/1‑149 (MQ=12)
                                          tcctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGAttt                                                                                             >  3:34582/1‑118 (MQ=32)
                                          tcctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGAttt                                                                                             <  4:34582/118‑1 (MQ=32)
                                            ctTAATCGTCGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTTGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc                                                             <  7:112160/148‑1 (MQ=12)
                                            ctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAATATTCg                                                            <  5:304488/149‑1 (MQ=17)
                                            ctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc                                                             <  2:303378/148‑1 (MQ=17)
                                            ctTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc                                                             <  2:257614/148‑1 (MQ=32)
                                             ttAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCg                                                            >  5:40978/1‑148 (MQ=17)
                                             ttAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGAttt                                                                                             >  4:209512/1‑115 (MQ=32)
                                               aaacgacgaA‑‑‑ATCAGCCTGCTGATGAGCTCCAACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGat                                                         <  8:135395/146‑1 (MQ=12)
                                               aaTCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCgg                                                           <  3:282993/147‑1 (MQ=17)
                                                aTCGACGAC‑‑‑CTCCGGCTGATGATGAGCTCCCACACTAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGat                                                         <  8:131000/148‑1 (MQ=11)
                                                aTCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTTTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGata                                                        <  5:237529/149‑1 (MQ=17)
                                                aTCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGat                                                         <  3:144520/148‑1 (MQ=32)
                                                 tCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGatat                                                       >  2:101747/1‑149 (MQ=17)
                                                 tCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGatat                                                       <  3:160746/149‑1 (MQ=32)
                                                   gacgac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtgt                                                     <  1:163775/149‑1 (MQ=17)
                                                   gacgac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtgt                                                     <  1:58104/149‑1 (MQ=32)
                                                   gacgac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtg                                                      <  2:35930/148‑1 (MQ=32)
                                                    acgac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtgt                                                     >  2:206887/1‑148 (MQ=32)
                                                     cgac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  3:35022/1‑149 (MQ=17)
                                                     cgac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  1:36341/1‑149 (MQ=32)
                                                      gac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgat                                                  >  3:302388/1‑149 (MQ=17)
                                                      gac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  4:175993/1‑148 (MQ=17)
                                                      gac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  2:281934/1‑148 (MQ=17)
                                                       ac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgata                                                 <  1:281934/149‑1 (MQ=17)
                                                       ac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgata                                                 <  8:249816/149‑1 (MQ=17)
                                                       ac‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCCTTAAAAATTCGGATATGTgata                                                 <  7:51504/149‑1 (MQ=17)
                                                            aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGCACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  5:323051/1‑145 (MQ=17)
                                                            aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  6:128859/1‑145 (MQ=17)
                                                            aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  8:249304/1‑145 (MQ=32)
                                                            aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                   >  2:269655/1‑145 (MQ=32)
                                                             tCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgataga                                               >  6:110154/1‑148 (MQ=17)
                                                             tCAGCAGGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagat                                              >  3:296443/1‑149 (MQ=17)
                                                               aGCCTGCTGATGAGCTCCCACACGTATTGCTTGATTCATTGTCGAAGACTTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTTTCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagatt                                             <  2:23377/148‑2 (MQ=255)
                                                               aGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagata                                             <  3:175993/148‑1 (MQ=17)
                                                               aGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagata                                             <  6:40978/148‑1 (MQ=17)
                                                                 ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAg                                          >  3:22339/1‑149 (MQ=32)
                                                                  ctgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGAAGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAg                                          >  1:69791/1‑148 (MQ=17)
                                                                   tgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGAc                                        >  4:247781/1‑149 (MQ=32)
                                                                       tgatgaGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                      >  1:290714/1‑148 (MQ=17)
                                                                             gCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAacac                                 >  1:33618/1‑147 (MQ=17)
                                                                             gCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa                               >  3:149820/1‑149 (MQ=32)
                                                                               tCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGt                             >  8:19192/1‑149 (MQ=32)
                                                                                 ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGtt                            >  3:62791/1‑148 (MQ=17)
                                                                                   acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                          >  6:102865/1‑148 (MQ=17)
                                                                                   acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                          >  4:311986/1‑148 (MQ=17)
                                                                                   acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                          >  1:237933/1‑148 (MQ=32)
                                                                                   acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa                         >  2:273838/1‑149 (MQ=32)
                                                                                     acGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa                         <  8:311758/147‑1 (MQ=32)
                                                                                       gAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgc                     >  7:72121/1‑149 (MQ=17)
                                                                                            gCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt                 <  2:36341/148‑1 (MQ=32)
                                                                                             cTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                >  4:283031/1‑148 (MQ=17)
                                                                                             cTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaa                                                                   >  4:51283/1‑97 (MQ=33)
                                                                                             cTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaa                                                                   <  3:51283/97‑1 (MQ=33)
                                                                                             cTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                >  3:168153/1‑148 (MQ=32)
                                                                                              ttGATCCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGa              >  6:323619/1‑149 (MQ=17)
                                                                                                  ttCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg          <  5:150001/149‑1 (MQ=17)
                                                                                                    cATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg         <  5:44775/148‑1 (MQ=17)
                                                                                                     aTTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGc        >  3:239209/1‑148 (MQ=17)
                                                                                                     aTTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGc        >  1:113635/1‑148 (MQ=32)
                                                                                                         tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAgg   <  2:174850/149‑1 (MQ=17)
                                                                                                         tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg    <  8:138696/148‑1 (MQ=17)
                                                                                                         tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg    >  7:138696/1‑148 (MQ=17)
                                                                                                         tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg    <  4:239209/148‑1 (MQ=17)
                                                                                                         tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAgg   >  5:237407/1‑149 (MQ=32)
                                                                                                          cGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAgg   >  1:52576/1‑148 (MQ=31)
                                                                                                           gAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGt  >  8:269709/1‑148 (MQ=11)
                                                                                                                                                    |                                                                                                          
TCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGT  >  NC_002947/1326977‑1327227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: