Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,743,687 T→C 7.1% Y745C (TAC→TGC) ‡ glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,743,6870TC 7.1% 380.4 / 12.8 126Y745C (TAC→TGC) ‡glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (56/61);  new base C (2/7);  total (58/68)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 129 total aligned reads displayed.

GTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATG  >  NC_002947/5743541‑5743836
                                                                                                                                                  |                                                                                                                                                     
gtgccgtgTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTAcc                                                                                                                                                    <  2:818693/150‑1 (MQ=255)
  gccgtgTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCgg                                                                                                                                                  <  1:2316139/150‑1 (MQ=255)
   ccgtgTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCggg                                                                                                                                                 >  2:649587/1‑150 (MQ=255)
   ccgtgTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCggg                                                                                                                                                 <  1:1816776/150‑1 (MQ=255)
   ccgtgTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGCCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTAcc                                                                                                                                                    <  1:1036544/147‑1 (MQ=255)
   ccgtgTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGCCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTAcc                                                                                                                                                    >  2:1036544/1‑147 (MQ=255)
    cgtgTTCAATACCGCCTTTCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGa                                                                                                                                                >  2:1306354/1‑150 (MQ=255)
     ttgtTCAATACCACCTGGCGCCACCGCTCCTGCGGTGTGTGCCACAGATTCTCGACCCGCCGCTGAGCATCCCCGTAGTCATCGAACTCGGCACACACCAGGAAGTGATCATAGGGCACCACGCCGTCGATCAACCCGGTGTACCGGGaa                                                                                                                                               <  2:385890/149‑1 (MQ=255)
     gtgtTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGaa                                                                                                                                               <  2:2407605/150‑1 (MQ=255)
     gtgtTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCgggtt                                                                                                                                               <  2:2050202/150‑3 (MQ=255)
     gtgtTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCgggtg                                                                                                                                               <  1:85194/150‑3 (MQ=255)
      tgtTCAATACCGCCTTGCGCCACCACTCCTGCGGGGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCGACCAGGCCGTCGATCAACCCGGTACACCgggttg                                                                                                                                              <  2:719854/150‑4 (MQ=255)
        ttCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCgggttgat                                                                                                                                            <  2:816421/150‑6 (MQ=255)
          cAATACCGCCTTGCGCCACCCCCCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGgtc                                                                                                                                          <  1:591718/150‑1 (MQ=255)
          cAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGgtc                                                                                                                                          <  1:2384087/150‑1 (MQ=255)
           aaTACCGCCTTGCGCCACCGCTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGgtcg                                                                                                                                         >  2:917903/1‑150 (MQ=255)
                      gcgcCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGaa                                                                                                                              <  2:1899464/150‑1 (MQ=255)
                      gcgcCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGaa                                                                                                                              <  1:742849/150‑1 (MQ=255)
                       cgcCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAac                                                                                                                             <  1:649587/150‑1 (MQ=255)
                        gcCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAaca                                                                                                                            >  2:1569391/1‑150 (MQ=255)
                         ccaccaCTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAacac                                                                                                                           <  2:1224546/150‑1 (MQ=255)
                          caccacTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATAGGCAACAACCAGTAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAacaca                                                                                                                          >  2:694732/1‑150 (MQ=255)
                           accacTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCTGGAAGGGTCGTCTGGCGAGAacacac                                                                                                                         <  1:1506988/150‑1 (MQ=255)
                             cacTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATTTCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACAcacca                                                                                                                       >  1:560932/1‑150 (MQ=255)
                              acTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCgggttgatc                                                                                                                                           <  1:690651/129‑7 (MQ=255)
                              acTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCgggttgatc                                                                                                                                           >  2:690651/1‑123 (MQ=255)
                                tCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGGGAGAACACACCACTg                                                                                                                    >  2:40430/1‑150 (MQ=255)
                                       gtgtgtGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATgg                                                                                                             <  2:1509072/150‑1 (MQ=255)
                                          tgtgCCACAGATCCTCGCCCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCt                                                                                                          <  1:379977/150‑1 (MQ=255)
                                          tgtgCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCt                                                                                                          <  1:1306354/150‑1 (MQ=255)
                                              ccACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCgggtt                                                                                                                                               <  1:1114058/109‑3 (MQ=255)
                                              ccACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTACACCgggtt                                                                                                                                               >  2:1114058/1‑107 (MQ=255)
                                               cacaGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATTGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAac                                                                                                     >  1:2295063/1‑150 (MQ=255)
                                               cacaGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAacacac                                                                                                                         <  1:863830/130‑1 (MQ=255)
                                                 caGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAacac                                                                                                   >  1:1011406/1‑150 (MQ=255)
                                                  aGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAAcaca                                                                                                  <  2:842802/150‑1 (MQ=255)
                                                      tcTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAACCGTCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCgg                                                                                              <  1:1832466/149‑1 (MQ=255)
                                                      ccTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCgg                                                                                              <  2:1759817/150‑1 (MQ=255)
                                                       cTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGt                                                                                             >  2:2376278/1‑150 (MQ=255)
                                                        tCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTGTGGAGACAAGATATGATTGCGAATGGCCTGCAACACATAGGta                                                                                            >  2:1832466/1‑149 (MQ=255)
                                                        tCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCt                                                                                                          >  2:1314132/1‑136 (MQ=255)
                                                        tCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTc                                                                                            >  1:1278713/1‑150 (MQ=255)
                                                        tCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTc                                                                                            >  2:1188872/1‑150 (MQ=255)
                                                        tCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTc                                                                                            >  2:384216/1‑150 (MQ=255)
                                                        tCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTc                                                                                            >  2:292188/1‑150 (MQ=255)
                                                          gACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCaa                                                                                          >  1:994921/1‑150 (MQ=255)
                                                             ccgccgCTGCGCATCCCCGTAGTCGTCGAAATCGGCACACACCCGGAAGGGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGTAGGGTAGTCTTGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCg                                                                                       >  1:557495/1‑150 (MQ=255)
                                                             ccgccgCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCTGCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCg                                                                                       >  1:1374632/1‑150 (MQ=255)
                                                             ccgccgCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCg                                                                                               <  2:533581/142‑1 (MQ=255)
                                                             ccgccgCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGAGAATGGCCTGCAACACATCg                                                                                               >  1:533581/1‑142 (MQ=255)
                                                                ccgcTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGtt                                                                                    <  1:2376278/150‑1 (MQ=255)
                                                                ccgcTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGtt                                                                                    <  1:1744841/150‑1 (MQ=255)
                                                                        cATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCgg                                                                            <  2:936683/150‑1 (MQ=255)
                                                                            ccAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTAt                                                                        <  1:1188872/150‑1 (MQ=255)
                                                                            ccAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTAt                                                                        <  1:384216/150‑1 (MQ=255)
                                                                             cAGTAGTCAGCGACATCGGCACACACCACGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATg                                                                       <  2:560932/150‑1 (MQ=255)
                                                                              agtagtCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATgg                                                                      >  1:1127994/1‑150 (MQ=255)
                                                                                tagtCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGTGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGgca                                                                    <  1:292188/150‑1 (MQ=255)
                                                                                   tCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGgcagca                                                                 >  1:1986411/1‑150 (MQ=255)
                                                                                    cATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGACGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAATCTTGCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCAt                                                                <  1:1398815/150‑1 (MQ=255)
                                                                                     aTCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCAtt                                                               >  1:505553/1‑150 (MQ=255)
                                                                                     aTCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCAtt                                                               >  2:2164866/1‑150 (MQ=255)
                                                                                      tCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTg                                                              >  1:330155/1‑150 (MQ=255)
                                                                                      tCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTg                                                              >  2:1856464/1‑150 (MQ=255)
                                                                                        gAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGgc                                                            <  2:330155/150‑1 (MQ=255)
                                                                                        gAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGgc                                                            <  1:917903/150‑1 (MQ=255)
                                                                                        gAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGgc                                                            <  2:156927/150‑1 (MQ=255)
                                                                                         aaaTCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGgcg                                                           >  1:314793/1‑150 (MQ=255)
                                                                                           aTCGGCACACACCAGGAAGCGATCATAGGCCACCTGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGcc                                                         <  2:1207580/150‑1 (MQ=255)
                                                                                                cacacaCCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTATGGCGAGGACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGACGCGGCTGTGGCAGCATTGGCGCAAAGAt                                                    >  2:963662/1‑150 (MQ=255)
                                                                                                   acacCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCaa                                                                                          >  2:1552070/1‑109 (MQ=255)
                                                                                                   acacCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCaa                                                                                          <  1:1552070/109‑1 (MQ=255)
                                                                                                   acacCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCAc                                                 >  1:1566874/1‑150 (MQ=255)
                                                                                                       cAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGc                                             >  1:1703552/1‑150 (MQ=255)
                                                                                                          gAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACg                                          >  1:2262748/1‑150 (MQ=255)
                                                                                                            aGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACAcacca                                                                                                                       <  2:248877/71‑1 (MQ=255)
                                                                                                             gCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCtt                                       >  1:1504907/1‑150 (MQ=255)
                                                                                                               gATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGtt                                                                                    >  1:1070313/1‑103 (MQ=255)
                                                                                                                    tAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCt                                >  1:1964716/1‑150 (MQ=255)
                                                                                                                     aGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCAc                                           <  1:1059673/138‑1 (MQ=255)
                                                                                                                     aGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCAc                                           >  2:1059673/1‑138 (MQ=255)
                                                                                                                       gCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCa                             >  2:1321073/1‑150 (MQ=255)
                                                                                                                         caccaGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAAc                           <  2:107619/150‑1 (MQ=255)
                                                                                                                         caccaGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAAc                           <  1:1569391/150‑1 (MQ=255)
                                                                                                                         caccaGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAAc                           <  2:1127994/150‑1 (MQ=255)
                                                                                                                         caccaGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAAc                           <  2:505553/150‑1 (MQ=255)
                                                                                                                         caccaGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAAc                           <  2:2295063/150‑1 (MQ=255)
                                                                                                                                   tCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAg                 >  1:1263046/1‑150 (MQ=255)
                                                                                                                                    cGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAgg                <  2:314793/150‑1 (MQ=255)
                                                                                                                                     gATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATGGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATGAACGGCACGCTTGCGCGCTTCAACCTGCTGCGAgat               >  1:1600277/1‑148 (MQ=255)
                                                                                                                                      aTCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGgca                                                                    >  1:1452600/1‑96 (MQ=255)
                                                                                                                                      aTCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTc              >  1:1262343/1‑150 (MQ=255)
                                                                                                                                        cAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAg            <  2:1263046/150‑1 (MQ=255)
                                                                                                                                         aaCCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAgg           <  2:1374632/150‑1 (MQ=255)
                                                                                                                                          aCCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCg                                               <  1:288062/113‑1 (MQ=255)
                                                                                                                                          aCCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCg                                               >  2:288062/1‑113 (MQ=255)
                                                                                                                                          aCCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGc          <  2:1566874/150‑1 (MQ=255)
                                                                                                                                                 gtACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGAt   <  2:1986411/150‑1 (MQ=255)
                                                                                                                                                 gtACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGCCGCGCTTTCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGAt   <  1:1856464/150‑1 (MQ=255)
                                                                                                                                                  tACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATg  <  2:1504409/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGCCGCTGAGCATCCCAGTAGTCATCGAAATCGGCACACACCAGGAAGCGATCATAGGCCACCAGGCCGTCGATCAACCCGGTGTACCGGGAAGGGTCGTCTGGCGAGAACACACCACTGCGAATGGCCTGCAACACATCGGTCAAGCGGTTGGAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATG  >  NC_002947/5743541‑5743836

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