Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,631,135 T→C 7.7% intergenic (+27/‑154) PP_1430 → / → lepA DegP‑like serine endoprotease/elongation factor 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,631,1350TC 7.7% 250.2 / 7.2 87intergenic (+27/‑154)PP_1430/lepADegP‑like serine endoprotease/elongation factor 4
Reads supporting (aligned to +/- strand):  ref base T (36/44);  new base C (4/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCCATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGTT  >  NC_002947/1630986‑1631283
                                                                                                                                                     |                                                                                                                                                    
gCCATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCt                                                                                                                                                      >  1:1017384/1‑150 (MQ=255)
    tCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCgg                                                                                                                                                  <  1:1940717/150‑1 (MQ=255)
    tCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCgg                                                                                                                                                  <  2:1375121/150‑1 (MQ=255)
    tCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCgg                                                                                                                                                  >  1:1375121/1‑150 (MQ=255)
       gTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTg                                                                                                                                               >  2:1114682/1‑150 (MQ=255)
         tCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGcc                                                                                                                                             <  1:671895/150‑1 (MQ=255)
         tCGGACAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGcc                                                                                                                                             <  1:478863/150‑1 (MQ=255)
             aaaggctgcgagccggacgcccgccgaaaaatagccgggaattgtaccttaagctgctgggttatgggagattcctctgtccgactgtacacagcaatcattccggcaattgtgaagattcatgaaAAAGGGCAGCTTCGGCTGCCCttt                                                                                                                                         >  1:322619‑M1/127‑150 (MQ=255)
              cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCtttt                                                                                                                                        <  1:1114682/150‑1 (MQ=255)
              cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCtttt                                                                                                                                        >  1:314861/1‑150 (MQ=255)
              cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCtttt                                                                                                                                        >  2:273441/1‑150 (MQ=255)
                  gAATTCACTGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCAt                                                                                                                                    <  1:1846393/150‑1 (MQ=255)
                    aTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGa                                                                                                                                  <  2:1206953/150‑1 (MQ=255)
                        aCCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCt                                                                                                                              <  2:124894/150‑1 (MQ=255)
                        aCCGAAATCGCCAAGGAACTGCCGAAGAATCGTTACGTGTCGATGCGCGTACTGCGTCAGGGACGTGCCAGCTTCATCACTTTCAAGCTTGCTAAATAAGCAGTTTCGCAAGGCGAAATGGCAGCTTCGGCTGCCCTTTTTCATCGATCt                                                                                                                              >  2:725908/1‑150 (MQ=255)
                          cGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTATCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTc                                                                                                                            <  2:142922/150‑1 (MQ=255)
                            aaaTCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTcac                                                                                                                          <  1:1083324/150‑1 (MQ=255)
                                cGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAAtt                                                                                                                      <  1:429458/150‑1 (MQ=255)
                                 gCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTg                                                                                                                     >  2:519861/1‑150 (MQ=255)
                                  ccAAGGAACTGCCGAAGCATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGc                                                                                                                    >  2:989307/1‑150 (MQ=255)
                                   gccgaaaaatagccgggaattgtaccttaagctgctgggttatgggagattcctctgtccgactgtacacagcaatcattccggcaattgtgaagattcatgaaAAAGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTGcc                                                                                                                   >  2:357065‑M1/105‑150 (MQ=255)
                                   cAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGcc                                                                                                                   <  1:1397934/150‑1 (MQ=255)
                                           tGCCGAAGAATCGTTCCGTGTCGATGCGCGTAATGAATCAGGGCCGTGCCTCCTTTATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGAt                                                                                                           >  1:339480/1‑150 (MQ=255)
                                             ccGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTg                                                                                                         >  2:1503028/1‑150 (MQ=255)
                                              cGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGc                                                                                                        <  1:273441/150‑1 (MQ=255)
                                               gaagaaTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCt                                                                                                       >  1:1873766/1‑150 (MQ=255)
                                                  gaaTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCACCCAAAGCTCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    <  1:118757/150‑1 (MQ=255)
                                                     tCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTAc                                                                                                 >  2:1620579/1‑150 (MQ=255)
                                                       gTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAg                                                                                               <  2:314861/150‑1 (MQ=255)
                                                            gtgtCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAg                                                                                               >  2:1167248/1‑145 (MQ=255)
                                                            gtgtCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAg                                                                                               <  1:1167248/145‑1 (MQ=255)
                                                               tCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTc                                                                                             <  1:292219/144‑1 (MQ=255)
                                                               tCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTc                                                                                             >  2:292219/1‑144 (MQ=255)
                                                                cGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACaga                                                                                      <  1:989307/150‑1 (MQ=255)
                                                                   tGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCCAGCTTGCTGAATAAGCAGTTTCCCAAGGCGGAAGTGTAGCCTCTGCTGCCCTTTTTTATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGa                                                                                   <  2:708156/150‑1 (MQ=255)
                                                                         tACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTcc                                                                             <  1:1700344/150‑1 (MQ=255)
                                                                          aCTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTccc                                                                            <  1:1053727/150‑1 (MQ=255)
                                                                          aCTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTccc                                                                            >  2:815799/1‑150 (MQ=255)
                                                                           cTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGTCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCa                                                                           >  1:678875/1‑150 (MQ=255)
                                                                            tGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGcc                                                                                                                                             <  1:1165362/83‑1 (MQ=255)
                                                                            tGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGcc                                                                                                                                             >  2:1165362/1‑83 (MQ=255)
                                                                             gCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                              <  2:888310/145‑1 (MQ=255)
                                                                             gCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                              >  1:888310/1‑145 (MQ=255)
                                                                             gCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttc                                                                                                                                        >  2:110091/1‑86 (MQ=38)
                                                                              ggggggggccgggccAGCTTCATCACGTTCAAGCTTGGTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATaa                                                                        <  2:339480/145‑1 (MQ=255)
                                                                              cGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATaa                                                                        <  1:519861/150‑1 (MQ=255)
                                                                              cGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATaa                                                                        <  2:261690/150‑1 (MQ=255)
                                                                                tCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAAcc                                                                      >  1:1852789/1‑150 (MQ=255)
                                                                                  aGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCca                                                                    <  2:1852789/150‑1 (MQ=255)
                                                                                      ccGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagcag                                                                >  2:495478/1‑150 (MQ=255)
                                                                                       cGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCagcagc                                                               <  2:880848/150‑1 (MQ=255)
                                                                                           ccAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                              >  1:203454/1‑131 (MQ=255)
                                                                                           ccAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                              <  2:203454/131‑1 (MQ=255)
                                                                                            cAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCa                                                                                                                                     <  2:28636/75‑1 (MQ=255)
                                                                                            cAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCa                                                                                                                                     >  1:28636/1‑75 (MQ=255)
                                                                                            cAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGa                                                                                          <  1:1591979/118‑1 (MQ=255)
                                                                                            cAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGa                                                                                          >  2:1591979/1‑118 (MQ=255)
                                                                                               cTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTa                                                       <  2:1873766/150‑1 (MQ=255)
                                                                                               cTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTa                                                       <  1:815799/150‑1 (MQ=255)
                                                                                                ttCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGa                                                                                                            >  2:648170/1‑96 (MQ=255)
                                                                                                ttCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGa                                                                                                            <  1:648170/96‑1 (MQ=255)
                                                                                                 tcatcaCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACa                                                     <  1:900195/150‑1 (MQ=255)
                                                                                                   atcaCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAAt                                                   >  2:544368/1‑150 (MQ=255)
                                                                                                   atcaCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAAt                                                   >  1:1647987/1‑150 (MQ=255)
                                                                                                    tcaCTTTCAAGCTTGCTGAATAAGCAGTTTCGCCAGGCGAAAGGGCAGCTTTGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCgg                                                                                           >  1:1910716/1‑109 (MQ=255)
                                                                                                    tcaCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTTGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCgg                                                                                           <  2:1910716/109‑1 (MQ=255)
                                                                                                        tttCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCg                                              >  2:1265924/1‑150 (MQ=255)
                                                                                                        tttCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCg                                              >  2:1227144/1‑150 (MQ=255)
                                                                                                         ttCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCgg                                             <  1:1265924/150‑1 (MQ=255)
                                                                                                         ttCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCgg                                             <  1:544368/150‑1 (MQ=255)
                                                                                                          tCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAAACTTCACAATTGCCGGAATGCTTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGc                                            >  1:87198/1‑150 (MQ=255)
                                                                                                            aaGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTa                                          <  2:1001803/150‑1 (MQ=255)
                                                                                                              gCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACGATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTAtt                                        >  2:531343/1‑150 (MQ=255)
                                                                                                                  tccggcaattgtgaagattcatgaaAAAGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACaga                                                                                      <  1:864561‑M1/75‑1 (MQ=255)
                                                                                                                  tccggcaattgtgaagattcatgaaAAAGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACaga                                                                                      >  2:864561‑M1/26‑100 (MQ=255)
                                                                                                                    tGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcgg                                  >  1:371128/1‑150 (MQ=255)
                                                                                                                      aaTAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcggcg                                <  2:1631740/150‑1 (MQ=255)
                                                                                                                       agaaGCAGGCTCGCAAGGCGAAAGGGCAGCTGCGGCTGCGCTTTTGCATGAATCTGCACAATTGCCAGAATGATTGCTGTGTACAGTCGTACAGAGGAAGCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcggcgg                               <  2:87198/148‑1 (MQ=255)
                                                                                                                            cAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTc                          <  2:1647987/150‑1 (MQ=255)
                                                                                                                             aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTc                                    >  1:800279/1‑139 (MQ=255)
                                                                                                                             aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTc                                    <  2:800279/139‑1 (MQ=255)
                                                                                                                              gaagattcatgaaAAAGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAAccc                                                                     <  2:1316654‑M1/92‑1 (MQ=255)
                                                                                                                              gaagattcatgaaAAAGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAAccc                                                                     >  1:1316654‑M1/14‑105 (MQ=255)
                                                                                                                                     aaGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                              <  2:213988/89‑1 (MQ=255)
                                                                                                                                     aaGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                              >  1:213988/1‑89 (MQ=255)
                                                                                                                                          gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCtt            >  1:1230502/1‑150 (MQ=255)
                                                                                                                                           AAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCttt           <  2:1760683‑M1/150‑1 (MQ=255)
                                                                                                                                                    CTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGtt  >  1:495478‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGTT  >  NC_002947/1630986‑1631283

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: