Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 127,123 G→C 7.2% intergenic (‑33/+37) thrB ← / ← PP_0122 homoserine kinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947127,1230GC 7.2% 747.6 / 35.6 223intergenic (‑33/+37)thrB/PP_0122homoserine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (96/111);  new base C (10/6);  total (106/117)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 240 total aligned reads displayed.

GAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTTTCG  >  NC_002947/126974‑127272
                                                                                                                                                     |                                                                                                                                                     
gaaATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTccgccg                                                                                                                                                       <  2:5085528/150‑1 (MQ=255)
 aaaTTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCgg                                                                                                                                                      <  2:2650288/150‑1 (MQ=255)
   aTTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGa                                                                                                                                                    <  2:4445850/150‑1 (MQ=255)
      gCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGt                                                                                                                                                 <  1:2576410/150‑1 (MQ=255)
        tGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGc                                                                                                                                               <  2:3610242/150‑1 (MQ=255)
          ttCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCg                                                                                                                                             >  2:1368974/1‑150 (MQ=255)
          ttCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCg                                                                                                                                             <  1:2737253/150‑1 (MQ=255)
           tctcGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCgg                                                                                                                                            <  2:2253495/150‑1 (MQ=255)
             tcGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGt                                                                                                                                          >  1:2854151/1‑150 (MQ=255)
              cGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGtt                                                                                                                                         <  1:5397464/150‑1 (MQ=255)
              cGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGtt                                                                                                                                         <  1:3048412/150‑1 (MQ=255)
               ggTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTg                                                                                                                                        <  2:1207110/150‑1 (MQ=255)
               ggTACCGGCGGCAATGCCCTGGAAGTCAAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGcccgtatg                                                                                                                                        <  2:3831459/150‑7 (MQ=255)
                   ccGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGc                                                                                                                                               >  2:1553944/1‑139 (MQ=255)
                   ccGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGc                                                                                                                                               >  2:4390396/1‑139 (MQ=255)
                   ccGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAtgt                                                                                                                                    <  2:4833078/150‑1 (MQ=255)
                      gcggcAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGaa                                                                                                                                 <  2:1050085/150‑1 (MQ=255)
                          cAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGt                                                                                                                             <  1:4887974/150‑1 (MQ=255)
                          cAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGt                                                                                                                             <  1:3121606/150‑1 (MQ=255)
                           aaTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTg                                                                                                                            >  1:2308586/1‑150 (MQ=255)
                           aaTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTg                                                                                                                            >  2:4950711/1‑150 (MQ=255)
                            aTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACCTGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTgc                                                                                                                           <  2:1839706/150‑1 (MQ=255)
                             tGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTgcg                                                                                                                          <  1:3792260/150‑1 (MQ=255)
                              gCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTgcgc                                                                                                                         <  2:1434616/150‑1 (MQ=255)
                                 cTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGccc                                                                                                                                             <  1:624050/127‑2 (MQ=255)
                                  tGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCa                                                                                                                     >  2:146858/1‑150 (MQ=255)
                                    gAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGc                                                                                                                   <  2:1725702/150‑1 (MQ=255)
                                       gTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTa                                                                                                                <  1:897769/150‑1 (MQ=255)
                                        tCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAc                                                                                                               <  2:2513758/150‑1 (MQ=255)
                                         cGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAcc                                                                                                              <  2:1473228/150‑1 (MQ=255)
                                         cGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAcc                                                                                                              <  2:2109186/150‑1 (MQ=255)
                                         cGAGCAGACGGCCCAGCTCGAACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAcc                                                                                                              <  2:1760995/150‑1 (MQ=255)
                                          gAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCa                                                                                                             >  2:4515552/1‑150 (MQ=255)
                                          gAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCa                                                                                                             >  1:1727493/1‑150 (MQ=255)
                                            gCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACt                                                                                                           >  2:2966905/1‑150 (MQ=255)
                                            gCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACt                                                                                                           <  1:4544564/150‑1 (MQ=255)
                                             cAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTc                                                                                                          <  2:57215/150‑1 (MQ=255)
                                                  ggCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGa                                                                                                     <  2:1930946/150‑1 (MQ=255)
                                                  ggCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGa                                                                                                     >  2:3196720/1‑150 (MQ=255)
                                                   gCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGAt                                                                                                    >  1:746237/1‑150 (MQ=255)
                                                    cccAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATc                                                                                                   <  2:3361776/150‑1 (MQ=255)
                                                     ccAGCTCGTACGGCGCCAGACAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCt                                                                                                  <  1:1707834/150‑1 (MQ=255)
                                                      cAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCtt                                                                                                 >  2:3778701/1‑150 (MQ=255)
                                                       aGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTc                                                                                                <  1:1340850/150‑1 (MQ=255)
                                                       aGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTc                                                                                                <  1:655341/150‑1 (MQ=255)
                                                         cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa                                                                                              <  1:4499492/150‑1 (MQ=255)
                                                         cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa                                                                                              <  1:2111821/150‑1 (MQ=255)
                                                         cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa                                                                                              >  2:2938417/1‑150 (MQ=255)
                                                             tACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGAc                                                                                          <  2:5165181/150‑1 (MQ=255)
                                                              aCGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACg                                                                                         <  1:486450/150‑1 (MQ=255)
                                                                 gcgcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGaa                                                                                      >  1:3620807/1‑150 (MQ=255)
                                                                   gcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATc                                                                                    <  2:3075250/150‑1 (MQ=255)
                                                                    ccAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCa                                                                                   <  1:1581102/150‑1 (MQ=255)
                                                                    ccAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCa                                                                                   <  1:1108828/150‑1 (MQ=255)
                                                                        aaaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCAt                                                                               >  1:1106022/1‑150 (MQ=255)
                                                                        aaaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCAt                                                                               <  1:2254212/150‑1 (MQ=255)
                                                                             tttCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCg                                                                          >  1:422829/1‑150 (MQ=255)
                                                                                 cAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTg                                                                      <  1:5321917/150‑1 (MQ=255)
                                                                                  aGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTgg                                                                     >  2:3550250/1‑150 (MQ=255)
                                                                                    cTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTc                                                                   >  2:3819107/1‑150 (MQ=255)
                                                                                    cTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTc                                                                   <  1:4321981/150‑1 (MQ=255)
                                                                                     tCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCg                                                                  >  1:1038331/1‑150 (MQ=255)
                                                                                     tCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCg                                                                  >  2:208106/1‑150 (MQ=255)
                                                                                       aGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGa                                                                <  2:3937755/150‑1 (MQ=255)
                                                                                       aGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGa                                                                >  1:2012877/1‑150 (MQ=255)
                                                                                       aGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGa                                                                >  2:5278582/1‑150 (MQ=255)
                                                                                         gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt                                                                                                                                      >  1:3325243/1‑78 (MQ=255)
                                                                                         gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt                                                                                                                                      <  2:3325243/78‑1 (MQ=255)
                                                                                         gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATc                                                              <  1:2540452/150‑1 (MQ=255)
                                                                                            gggTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGaa                                                           <  1:4064343/150‑1 (MQ=255)
                                                                                             ggTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAAt                                                          >  2:826028/1‑150 (MQ=255)
                                                                                               tCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGa                                                        >  1:1527119/1‑150 (MQ=255)
                                                                                                   gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                    >  1:5495108/1‑150 (MQ=255)
                                                                                                   gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                    <  2:3330000/150‑1 (MQ=255)
                                                                                                      ggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGc                                                 <  2:3649977/150‑1 (MQ=255)
                                                                                                      ggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGc                                                 <  1:3514970/150‑1 (MQ=255)
                                                                                                         gAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGCGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATc                                              >  2:5492088/1‑150 (MQ=255)
                                                                                                                 aCATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCAt                                      <  1:826028/150‑1 (MQ=255)
                                                                                                                   atgatgAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCa                                    <  1:3791047/150‑1 (MQ=255)
                                                                                                                       tgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcgc                                <  2:1518169/150‑1 (MQ=255)
                                                                                                                              ttGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAGTCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCCTCAGCGCGTACCAg                         >  2:2895166/1‑150 (MQ=255)
                                                                                                                              ttGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAg                         >  2:1942692/1‑150 (MQ=255)
                                                                                                                               tGCCCGTATGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGa                        >  1:3359610/1‑150 (MQ=255)
                                                                                                                                  ccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCgg                                                                 <  1:5037976/106‑1 (MQ=255)
                                                                                                                                  ccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaaaa                     >  1:4944264/1‑150 (MQ=255)
                                                                                                                                    aTACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATa                   <  2:4083516/150‑1 (MQ=255)
                                                                                                                                      cccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAg                 >  1:5519232/4‑150 (MQ=255)
                                                                                                                                       ccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg                >  1:4980491/3‑150 (MQ=255)
                                                                                                                                       cGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg                <  2:3525601/150‑1 (MQ=255)
                                                                                                                                       cGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg                <  1:5346868/150‑1 (MQ=255)
                                                                                                                                       cGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg                <  1:146858/150‑1 (MQ=255)
                                                                                                                                        cggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGc               >  1:2235843/2‑150 (MQ=255)
                                                                                                                                           cACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTg            >  1:5456826/1‑150 (MQ=255)
                                                                                                                                              ttCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGa                        >  2:4960193/1‑135 (MQ=255)
                                                                                                                                              ttCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATCAGGCTTGCCg         >  1:3210723/1‑150 (MQ=255)
                                                                                                                                              ttCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCg         >  2:5352390/1‑150 (MQ=255)
                                                                                                                                                 cgccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGttt      <  1:2033522/150‑1 (MQ=255)
                                                                                                                                                  gccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGtttt     >  1:2064543/1‑150 (MQ=255)
                                                                                                                                                   ccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGttttt    >  1:3556573/1‑150 (MQ=255)
                                                                                                                                                    cgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTTTc   <  2:1110222/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTTTCG  >  NC_002947/126974‑127272

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