Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,422,855:1 +G 7.5% coding (970/2211 nt) PP_3038 → TraC domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,422,8551.G 7.5% 424.1 / 28.6 107coding (970/2211 nt)PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (49/50);  new base G (4/4);  total (53/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 = 4.19e-01

Only 100 of 111 total aligned reads displayed.

TGCTGCGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATGCGCCTTCG  >  NC_002947/3422707‑3423003
                                                                                                                                                     |                                                                                                                                                    
tgctgcGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTAc                                                                                                                                                     >  2:544461/1‑150 (MQ=255)
tgctgcGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTAc                                                                                                                                                     >  2:1492423/1‑150 (MQ=255)
   tgcGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCgg                                                                                                                                                  <  1:2006994/150‑1 (MQ=255)
      ggtggtTGCGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTgg                                                                                                                                               >  2:1932562/1‑150 (MQ=255)
         ggtTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGTAGCCGGTGGCtt                                                                                                                                            <  2:1931022/150‑2 (MQ=255)
          gtTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTgg                                                                                                                                           >  2:1740756/1‑150 (MQ=255)
            tgtgGGTGGCAAGGTAGTCGGCCCGATCTGCTCGGCTGAGAGTGGGGTTAAATGGACCGACTTCAACTACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGGTCAAGCCACGGGCTACCGGTGGCTGTaa                                                                                                                                         <  2:510768/150‑1 (MQ=255)
             gtgggtggCCAGGTAGTCGGCCAGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAg                                                                                                                                        <  1:599041/150‑1 (MQ=255)
              ggggtggCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAgg                                                                                                                                       <  1:1932562/149‑1 (MQ=255)
                  tggCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCCGGAAGGACAt                                                                                                                                   <  1:967738/150‑1 (MQ=255)
                         gTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATg                                                                                                                                  >  2:932469/1‑144 (MQ=255)
                         gTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATg                                                                                                                                  <  1:932469/144‑1 (MQ=255)
                         gTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGcc                                                                                                                            <  1:182064/150‑1 (MQ=255)
                             tCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGa                                                                                                                        <  2:1183478/150‑1 (MQ=255)
                              cGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGAc                                                                                                                       <  1:1740756/150‑1 (MQ=255)
                               ggCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGa                                                                                                                                      >  2:696177/1‑134 (MQ=255)
                               ggCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGa                                                                                                                                      <  1:696177/134‑1 (MQ=255)
                                gCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTg                                                                                                                     >  1:926177/1‑150 (MQ=255)
                                 cccGATATTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTAAACGCCCTCCATTGCGCTGACTGCCTTGAGGCCGTTCCCCGGGAAGCCACGGAAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCCGGCCAGGACTgg                                                                                                                    <  2:586987/150‑2 (MQ=255)
                                  ccGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTgag                                                                                                                   >  2:1340935/1‑150 (MQ=255)
                                    gATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGagag                                                                                                                 <  2:1318351/150‑1 (MQ=255)
                                      tcttctCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtg                                                                                                               <  1:1492423/150‑1 (MQ=255)
                                      tcttctCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtg                                                                                                               <  2:1099826/150‑1 (MQ=255)
                                        ttctCGGCTGAGCGTGAGTTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtggt                                                                                                             <  1:230254/150‑1 (MQ=255)
                                         tctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCCAGCCACCGGTAGCCGGTGGCt                                                                                                                                             >  1:898487/1‑117 (MQ=255)
                                         tctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCCAGCCACCGGTAGCCGGTGGCt                                                                                                                                             <  2:898487/117‑1 (MQ=255)
                                             ggCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCt                                                                                                        <  1:842728/150‑1 (MQ=255)
                                                tGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGAt                                                                                                     >  1:2284946/1‑150 (MQ=255)
                                                   gCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGc                                                                                                  <  1:1633678/150‑1 (MQ=255)
                                                    cGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGcc                                                                                                 >  1:307177/1‑150 (MQ=255)
                                                     gTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGccc                                                                                                >  2:924189/1‑150 (MQ=255)
                                                     gTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGccc                                                                                                >  1:1232145/1‑150 (MQ=255)
                                                     gTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCCTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGccc                                                                                                <  1:1005040/150‑1 (MQ=255)
                                                         ggTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTAc                                                                                                                                                     <  1:449392/93‑1 (MQ=255)
                                                         ggTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTAc                                                                                                                                                     >  2:449392/1‑93 (MQ=255)
                                                          gTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGAc                                                                                                                                     <  1:1207170/108‑1 (MQ=255)
                                                          gTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGAc                                                                                                                                     >  2:1207170/1‑108 (MQ=255)
                                                          gTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATg                                                                                           >  2:1563987/1‑150 (MQ=255)
                                                           ttAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGc                                                                                          >  2:684484/1‑150 (MQ=255)
                                                           ttAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGACTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGc                                                                                          >  2:561671/1‑150 (MQ=255)
                                                            tAAATGGACCGTCGTCAACGTCCTTCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGACCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCa                                                                                         <  2:444657/150‑1 (MQ=255)
                                                              aaTGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACAGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGa                                                                                       >  1:343851/1‑150 (MQ=255)
                                                               aTGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGaa                                                                                      >  1:1928131/1‑150 (MQ=255)
                                                                tGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAAc                                                                                     <  1:1498808/150‑1 (MQ=255)
                                                                tGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAAc                                                                                     <  2:926177/150‑1 (MQ=255)
                                                                  gACCGACTTCAACGACCTGCATTGCGCTGAAAGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCACGCCCAGGGATGTACACACCGCCAGTTACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCGCATATGCAGAACGt                                                                                   >  1:2280505/1‑150 (MQ=255)
                                                                   aCCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGGTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTc                                                                                  <  2:358559/150‑1 (MQ=255)
                                                                   aCCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTc                                                                                  >  1:1456372/1‑150 (MQ=255)
                                                                      gACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGaa                                                                               <  2:1689103/150‑1 (MQ=255)
                                                                       aCTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAAc                                                                              >  1:1052595/1‑150 (MQ=255)
                                                                         ttCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTg                                                                            >  1:1339919/1‑150 (MQ=255)
                                                                         ttCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGCTGCCCTGATCGCCTATATGCAGAACGTCGAACTg                                                                            >  1:1337213/1‑150 (MQ=255)
                                                                          tCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGa                                                                                                                                      >  1:1453758/1‑91 (MQ=255)
                                                                          tCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGa                                                                                                                                      <  2:1453758/91‑1 (MQ=255)
                                                                          tCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGa                                                                           >  2:2076947/1‑150 (MQ=255)
                                                                           cAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGAt                                                                          >  2:84183/1‑150 (MQ=255)
                                                                                aCCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTgg                                                                     <  2:1339919/150‑1 (MQ=255)
                                                                                  cTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGcc                                                                   >  2:1162335/1‑150 (MQ=255)
                                                                                    tcATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCaa                                                                 <  2:963155/149‑1 (MQ=255)
                                                                                           gcTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGt                                                                                                                                                 >  2:949387/1‑63 (MQ=255)
                                                                                           gcTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGt                                                                                                                                                 <  1:949387/63‑1 (MQ=255)
                                                                                           gcTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTg                                                                                                                               >  1:1272051/1‑81 (MQ=255)
                                                                                             tGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAgcg                                                        >  1:1437543/1‑150 (MQ=255)
                                                                                               aaGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCt                                                      >  2:2241577/1‑150 (MQ=255)
                                                                                                 ggCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGGGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTgg                                                    >  1:1429248/1‑150 (MQ=255)
                                                                                                  gCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGa                                                   >  1:971544/1‑150 (MQ=255)
                                                                                                  gCCTTCAGGCCGTCCTCCGGCAAGTCACGGCAGTGATCAAGCCAGAGGCTACCGGTGCCTGAAAGGATATGCTGGCCAGGACTGAGAGTGATGCCCTGCTTGTCCATATGCAGAACGTAGAACTTATCCTGGCCAATGATGAGCGCTGGa                                                   <  2:959030/150‑1 (MQ=255)
                                                                                                    cTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtg                                                                                                               <  1:302235/88‑1 (MQ=255)
                                                                                                    cTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtg                                                                                                               >  2:302235/1‑88 (MQ=255)
                                                                                                    cTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAg                                                          <  1:1250970/141‑1 (MQ=255)
                                                                                                    cTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAg                                                          >  2:1250970/1‑141 (MQ=255)
                                                                                                    cTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGc                                                 >  1:616011/1‑150 (MQ=255)
                                                                                                             gTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGaa                                                                                                                                         >  1:2290915/1‑53 (MQ=255)
                                                                                                             gTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGaa                                                                                                                                         <  2:2290915/53‑1 (MQ=255)
                                                                                                                  ccGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTa                                   <  2:1456372/150‑1 (MQ=255)
                                                                                                                      cAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAAc                               <  2:1437543/150‑1 (MQ=255)
                                                                                                                       aaGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACg                              >  2:1972101/1‑150 (MQ=255)
                                                                                                                         gTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGcc                            <  2:2280505/150‑1 (MQ=255)
                                                                                                                          tCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCt                           >  1:97849/1‑150 (MQ=255)
                                                                                                                                    cttccaGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCggt                                              >  1:643509/6‑120 (MQ=255)
                                                                                                                                    cttccaGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCggt                                              <  2:643509/116‑2 (MQ=255)
                                                                                                                                     ttccaGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCgg                                      <  2:1312533/124‑1 (MQ=255)
                                                                                                                                     ttccaGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACTTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCgg                                      >  1:1312533/5‑128 (MQ=255)
                                                                                                                                     aTCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatgat                                                             <  1:50887/105‑1 (MQ=255)
                                                                                                                                     aTCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatgat                                                             >  2:50887/1‑105 (MQ=255)
                                                                                                                                     aTCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAg                >  1:1731030/1‑150 (MQ=255)
                                                                                                                                         aGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCa            >  1:1501684/1‑150 (MQ=255)
                                                                                                                                         aGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGt                                     <  2:1588301/125‑1 (MQ=255)
                                                                                                                                         aGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGt                                     >  1:1588301/1‑125 (MQ=255)
                                                                                                                                           ccACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTcc                   >  1:1803976/1‑141 (MQ=255)
                                                                                                                                           ccACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTcc                   <  2:1803976/141‑1 (MQ=255)
                                                                                                                                           ccACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATg          <  1:2076947/150‑1 (MQ=255)
                                                                                                                                            cACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCaa                                                                 <  2:1045395/94‑1 (MQ=255)
                                                                                                                                            cACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCaa                                                                 >  1:1045395/1‑94 (MQ=255)
                                                                                                                                            cACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATgc         <  2:307177/150‑1 (MQ=255)
                                                                                                                                             aCCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATc                                                                         <  2:2151181/85‑1 (MQ=255)
                                                                                                                                              ccGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGCCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATgcgc       <  2:2284946/150‑1 (MQ=255)
                                                                                                                                                 gCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTg                                                     >  1:1932757/1‑101 (MQ=255)
                                                                                                                                                 gCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTg                                                     <  2:1932757/101‑1 (MQ=255)
                                                                                                                                                   tACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATGCGCCTTCg  >  2:1308841/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCTGCGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATGCGCCTTCG  >  NC_002947/3422707‑3423003

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: