Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,609,381 A→T 5.0% K29M (AAG→ATG) ‡ PP_4079 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,609,3810AT 5.0% 994.6 / 22.8 290K29M (AAG→ATG) ‡PP_4079hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (134/141);  new base T (8/7);  total (142/148)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

GACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4609233‑4609529
                                                                                                                                                    |                                                                                                                                                    
gACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAg                                                                                                                                                     >  1:719885/1‑150 (MQ=255)
gACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAg                                                                                                                                                     <  2:720589/150‑1 (MQ=255)
  ccAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGt                                                                                                                                                   >  1:3920625/1‑150 (MQ=255)
       ttCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGgcgc                                                                                                                                              <  2:8379053/150‑1 (MQ=255)
         cgtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTg                                                                                                                                            <  2:5852319/150‑1 (MQ=255)
           tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          <  1:6610890/150‑1 (MQ=255)
            cgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGAc                                                                                                                                         <  1:2675045/150‑1 (MQ=255)
              aTGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCa                                                                                                                                       <  1:1808656/150‑1 (MQ=255)
               tGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCAc                                                                                                                                      <  1:4165149/150‑1 (MQ=255)
                 aCCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTc                                                                                                                                    >  1:2222758/1‑150 (MQ=255)
                   cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGt                                                                                                                                  >  1:5787598/1‑150 (MQ=255)
                      gATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTAc                                                                                                                               >  1:7263466/1‑150 (MQ=255)
                       aTGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGTCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCTGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCTCTGGACCACTCGTTAcc                                                                                                                              <  2:3316164/150‑1 (MQ=255)
                          gCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAAt                                                                                                                           >  1:2161299/1‑150 (MQ=255)
                           ccAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGTCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTGGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                          >  2:1339981/1‑150 (MQ=255)
                            cAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCg                                                                                                                         >  2:4982851/1‑150 (MQ=255)
                            cAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCg                                                                                                                         <  1:506360/150‑1 (MQ=255)
                            cAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCg                                                                                                                         <  1:861705/150‑1 (MQ=255)
                             aaaTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCgg                                                                                                                        <  1:2448080/150‑1 (MQ=255)
                             aaaTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCgg                                                                                                                        <  1:7345850/150‑1 (MQ=255)
                                   gTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGgccgcc                                                                                                                  <  2:6019762/150‑1 (MQ=255)
                                    tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGccc                                                                                                                 >  1:4352001/1‑150 (MQ=255)
                                     cccGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCa                                                                                                                <  1:309335/150‑1 (MQ=255)
                                       cGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAg                                                                                                              >  1:6835666/1‑150 (MQ=255)
                                       cGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAg                                                                                                              <  1:8448857/150‑1 (MQ=255)
                                               gCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaaccc                                                                                                      >  1:2758941/1‑150 (MQ=255)
                                                 ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    <  2:8670862/150‑1 (MQ=255)
                                                  tCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaac                                                                                                   >  2:5203774/1‑150 (MQ=255)
                                                     gggggCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAacccaaccca                                                                                                >  1:533102/1‑150 (MQ=255)
                                                     gggggCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAacccaaccca                                                                                                >  2:390164/1‑150 (MQ=255)
                                                           ggTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcgc                                                                                          <  2:6339973/150‑1 (MQ=255)
                                                                 gAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAg                                                                                    >  1:1411979/1‑150 (MQ=255)
                                                                   gccgccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCa                                                                                  >  1:1503529/1‑150 (MQ=255)
                                                                    ccgccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCAt                                                                                 >  1:382432/1‑150 (MQ=255)
                                                                       ccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCaa                                                                              <  2:3629552/150‑1 (MQ=255)
                                                                          ggCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGt                                                                           <  1:8156183/150‑1 (MQ=255)
                                                                                  ggATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCaa                                                                   <  2:1264067/150‑1 (MQ=255)
                                                                                   gATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAAc                                                                  >  1:7063881/1‑150 (MQ=255)
                                                                                   gATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAAc                                                                  >  2:7629501/1‑150 (MQ=255)
                                                                                   gATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAAc                                                                  >  1:2884148/1‑150 (MQ=255)
                                                                                   gATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAAc                                                                  >  2:1760522/1‑150 (MQ=255)
                                                                                     ttGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGc                                                                >  1:2320957/1‑150 (MQ=255)
                                                                                       gcgcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTc                                                              <  2:7387815/150‑1 (MQ=255)
                                                                                        cgcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCa                                                             >  1:6668736/1‑150 (MQ=255)
                                                                                          ccTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAAc                                                           <  1:6591245/150‑1 (MQ=255)
                                                                                          ccTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAAc                                                           <  1:6128334/150‑1 (MQ=255)
                                                                                           ctctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          <  1:4802342/150‑1 (MQ=255)
                                                                                            tctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAg                                                         >  1:7860621/1‑150 (MQ=255)
                                                                                            tctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAg                                                         >  2:4343659/1‑150 (MQ=255)
                                                                                            tctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAg                                                         >  2:5572871/1‑150 (MQ=255)
                                                                                              tgctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCa                                                       >  2:1452164/1‑150 (MQ=255)
                                                                                               gctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCaa                                                      <  1:1954303/150‑1 (MQ=255)
                                                                                                ctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAc                                                     <  1:7711410/150‑1 (MQ=255)
                                                                                                 tggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            >  2:7772302/1‑142 (MQ=255)
                                                                                                 tggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAcc                                                    >  2:3470539/1‑150 (MQ=255)
                                                                                                      gtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAAc                                               >  2:1058704/1‑150 (MQ=255)
                                                                                                       tgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACa                                              <  2:6807625/150‑1 (MQ=255)
                                                                                                       tgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACa                                              <  2:6826473/150‑1 (MQ=255)
                                                                                                         ctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAgt                                            <  1:8742789/150‑1 (MQ=255)
                                                                                                            gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         >  2:8743527/1‑150 (MQ=255)
                                                                                                            gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         >  2:7500949/1‑150 (MQ=255)
                                                                                                               ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGAc                                      <  2:4489189/150‑1 (MQ=255)
                                                                                                               ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGAc                                      >  2:2658138/1‑150 (MQ=255)
                                                                                                                 aCTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGt                                    >  1:7098252/1‑150 (MQ=255)
                                                                                                                   tGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTgg                                  >  1:2557150/1‑150 (MQ=255)
                                                                                                                    gCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggc                                 <  2:5781837/150‑1 (MQ=255)
                                                                                                                     ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcg                                <  2:4355085/150‑1 (MQ=255)
                                                                                                                     ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcg                                >  2:4668752/1‑150 (MQ=255)
                                                                                                                      cAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcgg                               <  2:3656013/150‑1 (MQ=255)
                                                                                                                        ggTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGcc                             <  1:4656959/150‑1 (MQ=255)
                                                                                                                         gTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCa                                                                                     >  2:7259666/1‑93 (MQ=255)
                                                                                                                         gTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCa                            <  2:3553989/150‑1 (MQ=255)
                                                                                                                              aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCt                       >  1:4074350/1‑150 (MQ=255)
                                                                                                                              aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCt                       >  1:3525056/1‑150 (MQ=255)
                                                                                                                              aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCt                       >  1:8032849/1‑150 (MQ=255)
                                                                                                                              aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCt                       <  1:7566053/150‑1 (MQ=255)
                                                                                                                              aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCt                       >  1:1001509/1‑150 (MQ=255)
                                                                                                                                gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     <  2:2052583/150‑1 (MQ=255)
                                                                                                                                  ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  2:481283/150‑1 (MQ=255)
                                                                                                                                  ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCAAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  1:8269923/150‑1 (MQ=255)
                                                                                                                                   gtggtccAGCGCCACCTTGGTGGCTCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAACCGAACAGTGTCGACGTGGCGGCCAAACCTGagag                  >  2:274702/7‑150 (MQ=255)
                                                                                                                                   cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGagag                  >  2:2804760/1‑150 (MQ=255)
                                                                                                                                    ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGc                 >  2:2158540/1‑150 (MQ=255)
                                                                                                                                    ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGc                 >  2:2694725/1‑150 (MQ=255)
                                                                                                                                     tGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCa                >  2:5952208/1‑150 (MQ=255)
                                                                                                                                       tccAGCGCCACCTTGGTGGCTCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATc              >  1:5471401/3‑150 (MQ=255)
                                                                                                                                       tccAGCGCCACCTTGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACCGTGTCGACGTGGCGGCCATCCCTGAGAGCATc              >  2:6209743/3‑150 (MQ=255)
                                                                                                                                        gCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATcc             >  1:1252728/1‑150 (MQ=255)
                                                                                                                                        ccAGCGCCACCTTGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATcc             >  1:355488/2‑150 (MQ=255)
                                                                                                                                         cAGCGCCACCTTGGTGGCTCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCt            >  1:5207454/1‑150 (MQ=255)
                                                                                                                                         cAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCt            <  1:9108979/150‑1 (MQ=255)
                                                                                                                                           gagCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTgg          <  2:2432053/150‑1 (MQ=255)
                                                                                                                                            agCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGa         >  1:4702955/1‑150 (MQ=255)
                                                                                                                                            agCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGa         <  1:314871/150‑1 (MQ=255)
                                                                                                                                               caccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAcc      >  2:9112134/1‑150 (MQ=255)
                                                                                                                                                accaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACCGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc     >  2:5060053/1‑150 (MQ=255)
                                                                                                                                                accaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc     >  2:700227/1‑150 (MQ=255)
                                                                                                                                                 ccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCa    <  2:6843843/150‑1 (MQ=255)
                                                                                                                                                  caAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAg   >  1:2055110/1‑150 (MQ=255)
                                                                                                                                                   aaGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGc  >  2:1969024/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4609233‑4609529

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