Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,788,856 T→G 7.5% A252A (GCT→GCG yibQ → divergent polysaccharide deacetylase domain‑containing protein YibQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,788,8560TG 7.5% 275.7 / 11.0 92A252A (GCT→GCGyibQdivergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base T (45/40);  new base G (5/2);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTG  >  NC_000913/3788708‑3789005
                                                                                                                                                    |                                                                                                                                                     
cGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtg                                                                                                                                                      <  2:1284384/150‑1 (MQ=255)
 gaacattgccattgAGCATCCTCACCCTGCAACGGTACGCGTGATGCAACAGATGGTTAATAACCTGCCACCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                                                                                                                                                     >  2:313381/4‑150 (MQ=255)
 gACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                                                                                                                                                     >  2:321623/1‑150 (MQ=255)
 gACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                                                                                                                                                     >  1:425444/1‑150 (MQ=255)
   ccattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGc                                                                                                                                                   >  1:1552345/1‑150 (MQ=255)
    cattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGTTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGcc                                                                                                                                                  <  1:1491163/150‑1 (MQ=255)
    cattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGcc                                                                                                                                                  >  1:66196/1‑150 (MQ=255)
     attgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCt                                                                                                                                                 <  1:1546241/150‑1 (MQ=255)
        gccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACGCCGCCGAAAAACGCTGTGCCTGAt                                                                                                                                              >  2:994351/1‑150 (MQ=255)
                gCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgcg                                                                                                                                      >  2:38616/1‑150 (MQ=255)
                gCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgcg                                                                                                                                      <  2:1341007/150‑1 (MQ=255)
                 cATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGt                                                                                                                                     <  1:457714/150‑1 (MQ=255)
                   tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTaa                                                                                                                                   >  1:831758/1‑150 (MQ=255)
                   tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTaa                                                                                                                                   >  1:337883/1‑150 (MQ=255)
                      tCACCCTGCATCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtcc                                                                                                                                <  2:502064/150‑1 (MQ=255)
                      tCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCGGTGCATCAGGCACAGCGtttttc                                                                                                                                <  1:1138167/150‑6 (MQ=255)
                          ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttc                                                                                                                            >  1:1327682/1‑150 (MQ=255)
                          ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttc                                                                                                                            >  1:821028/1‑150 (MQ=255)
                           cTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttcc                                                                                                                           >  2:610338/1‑150 (MQ=255)
                           cTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttcc                                                                                                                           >  2:42601/1‑150 (MQ=255)
                           cTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttcc                                                                                                                           >  1:1627086/1‑150 (MQ=255)
                            tGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttccg                                                                                                                          <  1:321623/150‑1 (MQ=255)
                                aCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGc                                                                                                                      >  1:92538/1‑150 (MQ=255)
                                   gTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTg                                                                                                                   <  1:1117040/150‑1 (MQ=255)
                                     aCGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGaa                                                                                                                 <  1:69682/150‑1 (MQ=255)
                                        cgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCt                                                                                                              >  2:1187572/1‑150 (MQ=255)
                                            ttGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGc                                                                                                          <  2:832377/150‑1 (MQ=255)
                                             tGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATAGTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTTGCGTGAAGCTGTGCa                                                                                                         <  2:1627086/150‑1 (MQ=255)
                                             tGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCa                                                                                                         <  1:913215/150‑1 (MQ=255)
                                              gCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCaa                                                                                                        >  2:342343/1‑150 (MQ=255)
                                               cAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGc                                                                                                          <  2:379490/147‑1 (MQ=255)
                                               cAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGc                                                                                                          >  1:379490/1‑147 (MQ=255)
                                               cAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCaaa                                                                                                       >  1:261173/1‑150 (MQ=255)
                                                  cAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCg                                                                                                    <  2:337883/150‑1 (MQ=255)
                                                  cAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCg                                                                                                    >  2:438741/1‑150 (MQ=255)
                                                     aTGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCgaag                                                                                                 <  1:775975/150‑1 (MQ=255)
                                                         tttATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAAc                                                                                             >  1:1310423/1‑150 (MQ=255)
                                                         tttATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAAc                                                                                             <  2:345813/150‑1 (MQ=255)
                                                           tataACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCg                                                                                           >  1:32111/1‑150 (MQ=255)
                                                           tataACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCg                                                                                           <  2:1075179/150‑1 (MQ=255)
                                                             taACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGAt                                                                                         >  2:1037698/1‑150 (MQ=255)
                                                             taACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGAt                                                                                         <  2:1103894/150‑1 (MQ=255)
                                                                ccTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGTACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGa                                                                                      <  1:1597715/150‑1 (MQ=255)
                                                                ccTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCGGTGCATCAGGCACAGCGttt                                                                                                                                   <  2:312444/105‑3 (MQ=255)
                                                                ccTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCGGTGCATCAGGCACAGCGttt                                                                                                                                   >  1:312444/1‑103 (MQ=255)
                                                                        cAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGttt                                                                              >  1:911900/1‑150 (MQ=255)
                                                                          gACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTAt                                                                            <  2:1015038/150‑1 (MQ=255)
                                                                                 cGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATc                                                                     <  1:149282/150‑1 (MQ=255)
                                                                                  gCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCTGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCg                                                                    >  1:805389/1‑150 (MQ=255)
                                                                                  gCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACTCCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCg                                                                    <  1:1167356/150‑1 (MQ=255)
                                                                                      gTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGcttc                                                                <  1:1047952/150‑1 (MQ=255)
                                                                                       tGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGcttct                                                               <  2:1548806/150‑1 (MQ=255)
                                                                                        gAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGcttctt                                                              <  1:38616/150‑1 (MQ=255)
                                                                                          aaGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTg                                                            <  1:1071742/150‑1 (MQ=255)
                                                                                            gCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGaa                                                          <  2:1552345/150‑1 (MQ=255)
                                                                                             ccAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAg                                                         <  2:92538/150‑1 (MQ=255)
                                                                                             ccAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAg                                                         <  2:627819/150‑1 (MQ=255)
                                                                                             ccAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAg                                                         <  1:214711/150‑1 (MQ=255)
                                                                                                   gcTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTaa                                                   >  1:80030/1‑150 (MQ=255)
                                                                                                     ttGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGc                                                 <  1:1037698/150‑1 (MQ=255)
                                                                                                       gCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCATTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGa                                               >  1:212535/1‑150 (MQ=255)
                                                                                                        ccgaatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTACTCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGCACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGaa                                              >  1:352711/3‑150 (MQ=255)
                                                                                                        ccgaatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTCTATATGCTAATCGCTTCTTTGAAGTATTAAGCGaa                                              >  1:346760/3‑150 (MQ=255)
                                                                                                        cTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGaa                                              >  1:1527690/1‑150 (MQ=255)
                                                                                                         tgaatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAATCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGaaa                                             >  2:1009789/1‑150 (MQ=255)
                                                                                                          gaatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAg                                            >  1:1364824/1‑150 (MQ=255)
                                                                                                          gaatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAACGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAg                                            >  2:894315/1‑150 (MQ=255)
                                                                                                            atgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCTCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCa                                          <  2:1325464/150‑1 (MQ=255)
                                                                                                              gaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATc                                        >  2:1400370/1‑150 (MQ=255)
                                                                                                              gaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATc                                        <  2:420755/150‑1 (MQ=255)
                                                                                                               aaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCa                                       >  1:822971/1‑150 (MQ=255)
                                                                                                                aCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAg                                      >  1:780480/1‑150 (MQ=255)
                                                                                                                 ccGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGc                                     >  2:404698/1‑150 (MQ=255)
                                                                                                                  cGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGcc                                    <  1:1187572/150‑1 (MQ=255)
                                                                                                                   gCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCa                                   >  2:700806/1‑150 (MQ=255)
                                                                                                                      ggTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCagag                                >  1:1175375/1‑150 (MQ=255)
                                                                                                                      ggTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCagag                                >  2:442735/1‑150 (MQ=255)
                                                                                                                         tGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCAc                             >  1:1463869/1‑150 (MQ=255)
                                                                                                                               ttCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACGGAAGAAACCTATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTTTTAAGCGAAAGCATCAGCCAGAGTACGCTGAt                       <  2:821028/150‑1 (MQ=255)
                                                                                                                                 cTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCg                     <  2:780480/150‑1 (MQ=255)
                                                                                                                                   acacCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGtt                   >  2:912854/1‑150 (MQ=255)
                                                                                                                                      ccgccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTAc                >  1:160950/1‑150 (MQ=255)
                                                                                                                                       cgccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACt               <  2:1264144/150‑1 (MQ=255)
                                                                                                                                        gccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACtt              >  1:1436142/1‑150 (MQ=255)
                                                                                                                                         ccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTc             >  2:1292566/1‑150 (MQ=255)
                                                                                                                                         acggattaCGCGGTGCATCAGGCACAGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTc             >  1:799701/8‑150 (MQ=255)
                                                                                                                                          cggattaCGCGGTGCATCAGGCACAGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTcc            >  1:375055/7‑150 (MQ=255)
                                                                                                                                          cgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTcc            <  2:831758/150‑1 (MQ=255)
                                                                                                                                           ggattaCGCGGTGCATCAGGCACAGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCa           >  1:775387/6‑150 (MQ=255)
                                                                                                                                           gAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCa           <  2:1327682/150‑1 (MQ=255)
                                                                                                                                           gAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCa           <  1:1292566/150‑1 (MQ=255)
                                                                                                                                           gAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCa           >  2:766006/1‑150 (MQ=255)
                                                                                                                                             aaaaCGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGCAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGc         >  1:1089071/1‑150 (MQ=255)
                                                                                                                                             aaaaCGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGCACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGCGCACGCTGCTCGTTTACTTCCAGc         >  2:1453714/1‑150 (MQ=255)
                                                                                                                                              ttaCGCGGTGCATCAGGCACAGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCa        >  1:891731/3‑150 (MQ=255)
                                                                                                                                              aaaCGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCa        <  1:42601/150‑1 (MQ=255)
                                                                                                                                                    tgtgCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTg  <  1:1400370/150‑1 (MQ=255)
                                                                                                                                                    tgtgCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTg  <  2:160950/150‑1 (MQ=255)
                                                                                                                                                    tgtgCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTg  >  2:956884/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTG  >  NC_000913/3788708‑3789005

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