Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,607,459 A→T 7.1% 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,607,4590AT 7.1% 312.9 / 8.9 98K29M (AAG→ATG) ‡PP_4079hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/48);  new base T (4/3);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4607311‑4607607
                                                                                                                                                    |                                                                                                                                                    
gACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAg                                                                                                                                                     <  1:116043/150‑1 (MQ=255)
gACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAg                                                                                                                                                     >  2:1428876/1‑150 (MQ=255)
       ttCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGgcgc                                                                                                                                              <  1:1428876/150‑1 (MQ=255)
          gtataTGACCACTATGGCCAAAGCTGTCCCGTTGGCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTgg                                                                                                                                           <  1:2313691/146‑1 (MQ=255)
           tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          >  2:1304487/1‑150 (MQ=255)
           tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACGGCCGGGTTGCACGTCTGTTAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          >  2:1600222/1‑150 (MQ=255)
            cgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGAc                                                                                                                                         <  1:1986991/150‑1 (MQ=255)
            cgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTAGAGCAGCGCCACCTTGGTGGCTCTgctc                                                                                                                                         <  1:1052482/150‑4 (MQ=255)
             gATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGACCGTCTGGATTTCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGAcc                                                                                                                                        >  2:1081220/1‑150 (MQ=255)
                gACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACt                                                                                                                                     <  1:1304487/150‑1 (MQ=255)
                 aCCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTc                                                                                                                                    <  1:356053/150‑1 (MQ=255)
                  ccACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTCGAGCAGAGCCACCAAGGTGgcgc                                                                                                                                              >  1:221157/1‑139 (MQ=255)
                  ccACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTCGAGCAGAGCCACCAAGGTGgcgc                                                                                                                                              <  2:221157/139‑1 (MQ=255)
                   cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgctcaa                                                                                                                                       <  2:630184/145‑6 (MQ=255)
                   cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgctcaa                                                                                                                                       >  1:630184/1‑140 (MQ=255)
                   cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACt                                                                                                                                     >  2:1530955/1‑147 (MQ=255)
                   cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACt                                                                                                                                     <  1:1530955/147‑1 (MQ=255)
                     cGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCTTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTa                                                                                                                                <  1:1184553/150‑1 (MQ=255)
                          gCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAAt                                                                                                                           >  1:1662424/1‑150 (MQ=255)
                           ccAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCTTGGTGGCTCTgctcaa                                                                                                                                       <  2:2021729/137‑6 (MQ=255)
                           ccAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCTTGGTGGCTCTgctcaa                                                                                                                                       >  1:2021729/1‑132 (MQ=255)
                           ccAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                          >  2:205197/1‑150 (MQ=255)
                            cAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCg                                                                                                                         >  2:945566/1‑150 (MQ=255)
                            cAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCg                                                                                                                         >  2:740716/1‑150 (MQ=255)
                              aaTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCggg                                                                                                                       >  2:607685/1‑150 (MQ=255)
                                   gTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGTCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCATTGAGCAGAACCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGgccgcc                                                                                                                  >  1:2298711/1‑150 (MQ=255)
                                    tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCt                                                                                                                                             <  2:1245833/122‑1 (MQ=255)
                                    tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCt                                                                                                                                             >  1:1245833/1‑122 (MQ=255)
                                    tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGccc                                                                                                                 >  1:2204442/1‑150 (MQ=255)
                                          tGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAc                                                                                                           <  1:101868/150‑1 (MQ=255)
                                              cGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacc                                                                                                       >  1:1994647/1‑150 (MQ=255)
                                              cGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacc                                                                                                       >  2:1283998/1‑150 (MQ=255)
                                                 ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    <  1:1283998/150‑1 (MQ=255)
                                                 ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    <  1:607685/150‑1 (MQ=255)
                                                   ccGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaacc                                                                                                  <  1:945566/150‑1 (MQ=255)
                                                      ggggCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAAcccaacccaa                                                                                               <  1:205197/150‑1 (MQ=255)
                                                          cGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcg                                                                                           <  2:2315339/150‑1 (MQ=255)
                                                            gTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcgcg                                                                                         >  2:362236/1‑150 (MQ=255)
                                                             tCATGAGCCGCCTGGCGTTCTGGATTGCGCCTCTGCTGGTCGTGGTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGt                                                                                        >  1:1891464/1‑150 (MQ=255)
                                                             tCATGAGCCGCCTGGCCGTCTGGATTGCGCAACTGTTGGGCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGt                                                                                        <  2:1397471/150‑1 (MQ=255)
                                                                tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCa                                                                                     <  1:1600222/150‑1 (MQ=255)
                                                                tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCa                                                                                     <  1:1477741/150‑1 (MQ=255)
                                                                 taGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAg                                                                                    <  1:1097096/149‑1 (MQ=255)
                                                                     cgccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATc                                                                                <  1:618843/150‑1 (MQ=255)
                                                                     cgccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATc                                                                                <  2:1627387/150‑1 (MQ=255)
                                                                       ccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCaa                                                                              <  2:1448676/150‑1 (MQ=255)
                                                                        cTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAAc                                                                             <  2:1994647/150‑1 (MQ=255)
                                                                               tCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGa                                                                                                            <  1:1308709/112‑1 (MQ=255)
                                                                               tCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGa                                                                                                            >  2:1308709/1‑112 (MQ=255)
                                                                               tCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAc                                                                                                           <  2:1115454/113‑1 (MQ=255)
                                                                               tCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAc                                                                                                           >  1:1115454/1‑113 (MQ=255)
                                                                                 tGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAAt                                                                                                                           <  2:1634516/95‑1 (MQ=255)
                                                                                 tGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAAt                                                                                                                           >  1:1634516/1‑95 (MQ=255)
                                                                                         gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGCGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            >  2:144669/1‑150 (MQ=255)
                                                                                         gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcgc                                                                                          >  1:1350117/1‑120 (MQ=255)
                                                                                         gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcgc                                                                                          <  2:1350117/120‑1 (MQ=255)
                                                                                         gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            >  2:1502444/1‑150 (MQ=255)
                                                                                          ccTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAAc                                                           <  1:740716/150‑1 (MQ=255)
                                                                                           ctctGCTTGTCGTGCTGGTCTGACTGCCAGGTTGCACTTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCATGTCGCTTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          <  2:2298711/150‑1 (MQ=255)
                                                                                                      gtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCAACAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAAc                                               <  1:264449/150‑1 (MQ=255)
                                                                                                       tgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACa                                              <  2:258156/150‑1 (MQ=255)
                                                                                                            gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTAc                                                                                                                               <  1:1479685/64‑1 (MQ=255)
                                                                                                            gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTAc                                                                                                                               >  2:1479685/1‑64 (MQ=255)
                                                                                                            gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    >  2:1272374/1‑91 (MQ=255)
                                                                                                            gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    <  1:1272374/91‑1 (MQ=255)
                                                                                                               ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                          <  1:1352057/66‑1 (MQ=255)
                                                                                                               ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                          >  2:1352057/1‑66 (MQ=255)
                                                                                                               ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         <  2:237131/147‑1 (MQ=255)
                                                                                                               ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         >  1:237131/1‑147 (MQ=255)
                                                                                                               ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGAc                                      <  1:871432/150‑1 (MQ=255)
                                                                                                                gACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACg                                     >  2:2386922/1‑150 (MQ=255)
                                                                                                                   tGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTgg                                  <  2:2204442/150‑1 (MQ=255)
                                                                                                                     ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAAcccaacccaa                                                                                               <  2:1291368/87‑1 (MQ=255)
                                                                                                                     ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAAcccaacccaa                                                                                               >  1:1291368/1‑87 (MQ=255)
                                                                                                                     ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCCGCCCAACGCGCAACAGCa                                                       <  2:102018/127‑1 (MQ=255)
                                                                                                                     ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCa                                                       >  1:102018/1‑127 (MQ=255)
                                                                                                                     ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcg                                >  1:861010/1‑150 (MQ=255)
                                                                                                                      cAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcgg                               >  2:310928/1‑150 (MQ=255)
                                                                                                                                gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAgtgt                                          >  2:1163988/1‑129 (MQ=255)
                                                                                                                                gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAgtgt                                          <  1:1163988/129‑1 (MQ=255)
                                                                                                                                 tCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCCCTCTTGATCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgag                    <  2:2371125/150‑1 (MQ=255)
                                                                                                                                 tCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGa                                                                                                            >  2:1112105/1‑62 (MQ=255)
                                                                                                                                 tCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGa                                                                                                            <  1:1112105/62‑1 (MQ=255)
                                                                                                                                  ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  1:2220464/150‑1 (MQ=255)
                                                                                                                                   cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    >  1:911641/1‑68 (MQ=255)
                                                                                                                                   cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    <  2:911641/68‑1 (MQ=255)
                                                                                                                                   cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGagag                  >  1:1059215/1‑150 (MQ=255)
                                                                                                                                   cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGagag                  >  1:1059247/1‑150 (MQ=255)
                                                                                                                                    ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            <  2:52710/107‑1 (MQ=255)
                                                                                                                                    ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            >  1:52710/1‑107 (MQ=255)
                                                                                                                                    ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGc                 >  1:1280851/1‑150 (MQ=255)
                                                                                                                                       tccAGCGCCACCTTGGTGGCTCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGAGCTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATc              >  2:1067787/3‑150 (MQ=255)
                                                                                                                                       aGCATAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATc              <  1:275056/150‑1 (MQ=255)
                                                                                                                                         cAGCGCCACCTTGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCt            >  1:1987625/1‑150 (MQ=255)
                                                                                                                                               caccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAccc      >  2:806549/1‑150 (MQ=255)
                                                                                                                                                accaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc     >  1:1048051/1‑150 (MQ=255)
                                                                                                                                                accaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc     >  2:2017785/1‑150 (MQ=255)
                                                                                                                                                  caAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAg   <  1:1502444/150‑1 (MQ=255)
                                                                                                                                                   aaGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGc  >  2:352350/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4607311‑4607607

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