Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 127,121 C→G 7.1% intergenic (‑31/+39) 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,1210CG 7.1% 755.5 / 36.1 226intergenic (‑31/+39)thrB/PP_0122homoserine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (96/114);  new base G (10/6);  total (106/120)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 241 total aligned reads displayed.

AAGAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTTT  >  NC_002947/126972‑127270
                                                                                                                                                     |                                                                                                                                                     
aagaaATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTccgc                                                                                                                                                       <  1:2407624/150‑1 (MQ=255)
  gaaATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTccgccg                                                                                                                                                     <  2:5085528/150‑1 (MQ=255)
     aTTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGa                                                                                                                                                  <  2:4445850/150‑1 (MQ=255)
         cTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTg                                                                                                                                              >  2:142470/1‑150 (MQ=255)
            ttCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCg                                                                                                                                           <  1:2668532/150‑1 (MQ=255)
            ttCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCg                                                                                                                                           >  2:1368974/1‑150 (MQ=255)
            ttCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGccc                                                                                                                                           <  2:1065847/150‑2 (MQ=255)
            gtctcGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGccc                                                                                                                                           <  1:3803234/149‑2 (MQ=255)
                cGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTGTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGtt                                                                                                                                       <  1:3048470/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)
                            cAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGt                                                                                                                           <  1:4888500/150‑1 (MQ=255)
                             aaTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTg                                                                                                                          >  2:4950711/1‑150 (MQ=255)
                             aaTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTg                                                                                                                          >  1:2308586/1‑150 (MQ=255)
                              aTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACCTGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTgc                                                                                                                         <  2:1839706/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)
                                           cGAGCAGACGGCCCAGCTCGAACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAcc                                                                                                            <  2:1760995/150‑1 (MQ=255)
                                            gAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCa                                                                                                           >  1:1727493/1‑150 (MQ=255)
                                            gAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCa                                                                                                           >  2:4515552/1‑150 (MQ=255)
                                              gCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACt                                                                                                         >  2:1129702/1‑150 (MQ=255)
                                              gCAGACGGCCCAGCTCGGACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGACGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACt                                                                                                         >  2:2966928/1‑150 (MQ=255)
                                               cAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTc                                                                                                        <  2:57215/150‑1 (MQ=255)
                                                 gACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGa                                                                                                      <  1:458358/150‑1 (MQ=255)
                                                    ggCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGa                                                                                                   >  2:2355899/1‑150 (MQ=255)
                                                    ggCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGa                                                                                                   <  2:3933407/150‑1 (MQ=255)
                                                    ggCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGa                                                                                                   <  2:1930924/150‑1 (MQ=255)
                                                     gCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGAt                                                                                                  >  2:3321027/1‑150 (MQ=255)
                                                      cccAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATc                                                                                                 <  2:3361776/150‑1 (MQ=255)
                                                       ccAGCTCGTACGGCGCCAGACAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCt                                                                                                <  1:1707834/150‑1 (MQ=255)
                                                       ccAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCt                                                                                                >  1:4389627/1‑150 (MQ=255)
                                                        cAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCtt                                                                                               >  2:3778701/1‑150 (MQ=255)
                                                         aGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTc                                                                                              <  1:1340850/150‑1 (MQ=255)
                                                          gCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTcc                                                                                             <  1:5364413/150‑1 (MQ=255)
                                                           cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa                                                                                            <  1:2111821/150‑1 (MQ=255)
                                                           cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa                                                                                            <  1:4499492/150‑1 (MQ=255)
                                                           cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa                                                                                            >  2:2938543/1‑150 (MQ=255)
                                                               tACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGAc                                                                                        <  2:5165181/150‑1 (MQ=255)
                                                                aCGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACg                                                                                       <  1:486450/150‑1 (MQ=255)
                                                                    cgcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAAt                                                                                   >  1:705326/1‑150 (MQ=255)
                                                                     gcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATc                                                                                  <  2:3075250/150‑1 (MQ=255)
                                                                      ccAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCa                                                                                 <  1:1108828/150‑1 (MQ=255)
                                                                      ccAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCa                                                                                 <  1:1581102/150‑1 (MQ=255)
                                                                      ccAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCa                                                                                 >  2:3233778/1‑150 (MQ=255)
                                                                        aGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGc                                                                               <  2:4685539/150‑1 (MQ=255)
                                                                         gAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCa                                                                              <  1:4950711/150‑1 (MQ=255)
                                                                          aaaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCAt                                                                             >  1:1106022/1‑150 (MQ=255)
                                                                          aaaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCAt                                                                             <  1:2254212/150‑1 (MQ=255)
                                                                              gTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTcc                                                                         <  1:3371675/150‑1 (MQ=255)
                                                                                 tCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGc                                                                      <  2:1468455/150‑1 (MQ=255)
                                                                                  ccAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACTACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCt                                                                     <  2:2711321/150‑1 (MQ=255)
                                                                                   cAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTg                                                                    <  1:5321917/150‑1 (MQ=255)
                                                                                    aGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTgg                                                                   >  2:3550250/1‑150 (MQ=255)
                                                                                      cTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTc                                                                 >  2:3819107/1‑150 (MQ=255)
                                                                                       tCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCg                                                                >  1:1038331/1‑150 (MQ=255)
                                                                                       tCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCg                                                                >  2:208106/1‑150 (MQ=255)
                                                                                         aGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGa                                                              >  1:2012877/1‑150 (MQ=255)
                                                                                           gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt                                                                                                                                    >  1:3325243/1‑78 (MQ=255)
                                                                                           gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt                                                                                                                                    >  1:511525/1‑78 (MQ=255)
                                                                                             cGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGa                                                          <  1:541083/150‑1 (MQ=255)
                                                                                             cGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGa                                                          >  1:4252772/1‑150 (MQ=255)
                                                                                              gggTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGaa                                                         <  1:4064343/150‑1 (MQ=255)
                                                                                              gggTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGaa                                                         <  2:167422/150‑1 (MQ=255)
                                                                                               ggTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAAt                                                        >  2:826028/1‑150 (MQ=255)
                                                                                                 tCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGa                                                      >  1:1527119/1‑150 (MQ=255)
                                                                                                     gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGACGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                  >  2:3604182/1‑150 (MQ=255)
                                                                                                     gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                  >  1:5495292/1‑150 (MQ=255)
                                                                                                     gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                  <  2:3330000/150‑1 (MQ=255)
                                                                                                     gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                  >  1:3668236/1‑150 (MQ=255)
                                                                                                        ggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGc                                               <  1:3514970/150‑1 (MQ=255)
                                                                                                        ggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGc                                               <  2:3649977/150‑1 (MQ=255)
                                                                                                           gAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGCGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATc                                            >  2:5492088/1‑150 (MQ=255)
                                                                                                              gACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGt                                         <  1:2787584/150‑1 (MQ=255)
                                                                                                                   aCATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCAt                                    <  1:826057/150‑1 (MQ=255)
                                                                                                                     atgatgAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCa                                  <  1:3791047/150‑1 (MQ=255)
                                                                                                                         tgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcgc                              <  2:1517597/150‑1 (MQ=255)
                                                                                                                         tgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcgc                              >  2:2249894/1‑150 (MQ=255)
                                                                                                                          gaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGcgcg                             >  1:3002222/1‑150 (MQ=255)
                                                                                                                                ttGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAGTCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCCTCAGCGCGTACCAg                       >  2:2895166/1‑150 (MQ=255)
                                                                                                                                ttGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAg                       >  2:1942692/1‑150 (MQ=255)
                                                                                                                                      aTACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATa                 <  2:4083516/150‑1 (MQ=255)
                                                                                                                                       tACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATaa                <  2:5180262/150‑1 (MQ=255)
                                                                                                                                       tACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAAGGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATaa                <  2:3444339/150‑1 (MQ=255)
                                                                                                                                        cccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAg               >  1:5519232/4‑150 (MQ=255)
                                                                                                                                         ccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg              >  1:4980491/3‑150 (MQ=255)
                                                                                                                                         ccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg              >  1:2572525/3‑150 (MQ=255)
                                                                                                                                         ccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg              >  1:649703/3‑150 (MQ=255)
                                                                                                                                         cGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg              <  1:146858/150‑1 (MQ=255)
                                                                                                                                          cggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGc             >  1:2235843/2‑150 (MQ=255)
                                                                                                                                            gCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCtt           <  1:3308037/150‑1 (MQ=255)
                                                                                                                                             cACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTg          <  1:3605087/150‑1 (MQ=255)
                                                                                                                                             cACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTg          >  1:5458440/1‑150 (MQ=255)
                                                                                                                                             cACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTg          >  1:816276/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)
                                                                                                                                                     |                                                                                                                                                     
AAGAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTTT  >  NC_002947/126972‑127270

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: