Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,607,466 G→T 5.9% A31A (GCG→GCT PP_4079 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,607,4660GT 5.9% 326.2 / 5.3 102A31A (GCG→GCTPP_4079hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (45/51);  new base T (3/3);  total (48/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 105 total aligned reads displayed.

TTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGCCCCT  >  NC_002947/4607318‑4607611
                                                                                                                                                    |                                                                                                                                                 
   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                                                                                                                                                  <  2:221157/139‑1 (MQ=255)
           ccACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTCGAGCAGAGCCACCAAGGTGgcgc                                                                                                                                                  >  1:221157/1‑139 (MQ=255)
            cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgctcaa                                                                                                                                           >  1:630184/1‑140 (MQ=255)
            cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgctcaa                                                                                                                                           <  2:630184/145‑6 (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                                                                                                                                           >  1:2021729/1‑132 (MQ=255)
                    ccAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                              >  2:205197/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                                                                                                           >  2:1283998/1‑150 (MQ=255)
                                       cGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacc                                                                                                           >  1:1994647/1‑150 (MQ=255)
                                          ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                        <  1:607685/150‑1 (MQ=255)
                                          ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                        <  1:1283998/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                                                                                                               >  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                                                                                              <  2:1350117/120‑1 (MQ=255)
                                                                                  gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcgc                                                                                              >  1:1350117/1‑120 (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                                                                                                                                   >  2:1479685/1‑64 (MQ=255)
                                                                                                     gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTAc                                                                                                                                   <  1:1479685/64‑1 (MQ=255)
                                                                                                     gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                        <  1:1272374/91‑1 (MQ=255)
                                                                                                     gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                        >  2:1272374/1‑91 (MQ=255)
                                                                                                        ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                              >  2:1352057/1‑66 (MQ=255)
                                                                                                        ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                              <  1:1352057/66‑1 (MQ=255)
                                                                                                        ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                             <  2:237131/147‑1 (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                                                                                                        <  2:911641/68‑1 (MQ=255)
                                                                                                                            cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                        >  1:911641/1‑68 (MQ=255)
                                                                                                                            cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGagag                      >  1:1059215/1‑150 (MQ=255)
                                                                                                                            cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGagag                      >  1:1059247/1‑150 (MQ=255)
                                                                                                                             ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                                >  1:52710/1‑107 (MQ=255)
                                                                                                                             ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                                <  2:52710/107‑1 (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)
                                                                                                                                         accttggTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc         <  1:198465/145‑1 (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)
                                                                                                                                              ggTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGccc    <  1:2017785/150‑1 (MQ=255)
                                                                                                                                               gTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaac                                                                                                            >  2:315204/1‑45 (MQ=255)
                                                                                                                                               gTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaac                                                                                                            <  1:315204/45‑1 (MQ=255)
                                                                                                                                                tGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                         >  1:868782/1‑127 (MQ=255)
                                                                                                                                                tGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                         <  2:868782/127‑1 (MQ=255)
                                                                                                                                                tGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGCCCCt  <  2:1059215/150‑1 (MQ=255)
                                                                                                                                                tGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGCCCCt  <  2:1059247/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGCCCCT  >  NC_002947/4607318‑4607611

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: