Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,788,870 C→A 7.3% P257Q (CCG→CAG)  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,8700CA 7.3% 321.1 / 10.6 96P257Q (CCG→CAG) yibQdivergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base C (47/42);  new base A (5/2);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

GCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGC  >  NC_000913/3788724‑3789018
                                                                                                                                                  |                                                                                                                                                    
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:821028/1‑150 (MQ=255)
          ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttc                                                                                                                                         >  1:1327682/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                                                                                                                       >  1:379490/1‑147 (MQ=255)
                               cAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGc                                                                                                                       <  2:379490/147‑1 (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:1103894/150‑1 (MQ=255)
                                             taACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGAt                                                                                                      >  2:1037698/1‑150 (MQ=255)
                                                ccTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGTACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGa                                                                                                   <  1:1597715/150‑1 (MQ=255)
                                                ccTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCGGTGCATCAGGCACAGCGttt                                                                                                                                                >  1:312444/1‑103 (MQ=255)
                                                ccTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCGGTGCATCAGGCACAGCGttt                                                                                                                                                <  2:312444/105‑3 (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                                                                      <  1:214711/150‑1 (MQ=255)
                                                                             ccAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAg                                                                      <  2:627819/150‑1 (MQ=255)
                                                                             ccAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAg                                                                      <  2:92538/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:766006/1‑150 (MQ=255)
                                                                                                                           gAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCa                        <  1:1292566/150‑1 (MQ=255)
                                                                                                                           gAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCa                        <  2:1327682/150‑1 (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)
                                                                                                                                     gtgCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTgg              >  1:1017845/1‑150 (MQ=255)
                                                                                                                                       gCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCa            >  2:196575/1‑150 (MQ=255)
                                                                                                                                         cTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAg          <  1:268317/150‑1 (MQ=255)
                                                                                                                                             tGCACCGCGTAATCCGTTACGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTg      >  2:1498568/1‑150 (MQ=255)
                                                                                                                                              gCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTg               >  2:607669/1‑140 (MQ=255)
                                                                                                                                              gCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTg               <  1:607669/140‑1 (MQ=255)
                                                                                                                                              gCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTgg     <  1:438741/150‑1 (MQ=255)
                                                                                                                                              gCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTgg     >  1:175166/1‑150 (MQ=255)
                                                                                                                                                 ccGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGc  >  2:420743/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGC  >  NC_000913/3788724‑3789018

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: