Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,563,554 G→A 6.6% intergenic (+36/+170) glpD → / ← yzgL aerobic glycerol 3‑phosphate dehydrogenase/putative uncharacterized protein YzgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,563,5540GA 6.6% 283.2 / 6.8 91intergenic (+36/+170)glpD/yzgLaerobic glycerol 3‑phosphate dehydrogenase/putative uncharacterized protein YzgL
Reads supporting (aligned to +/- strand):  ref base G (35/49);  new base A (3/3);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTT  >  NC_000913/3563405‑3563697
                                                                                                                                                     |                                                                                                                                               
gacgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                 <  1:17793/150‑1 (MQ=255)
    cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCTGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTg                                                                                                                                             <  2:1609543/150‑1 (MQ=255)
    cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTGACCTg                                                                                                                                             <  2:1508232/150‑1 (MQ=255)
      cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTgag                                                                                                                                           >  1:1149147/1‑150 (MQ=255)
        gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTTGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGAc                                                                                                                                         <  2:889363/150‑1 (MQ=255)
          ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtg                                                                                                                                       >  1:1057620/1‑150 (MQ=255)
             gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                    <  2:1158330/150‑1 (MQ=255)
                gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                    >  2:1587721/1‑147 (MQ=255)
                gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                    >  2:461375/1‑147 (MQ=255)
                gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                    <  1:1587721/147‑1 (MQ=255)
                gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                    <  1:461375/147‑1 (MQ=255)
                 cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACaa                                                                                                                                <  1:570050/150‑1 (MQ=255)
                 cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACaa                                                                                                                                <  2:100329/150‑1 (MQ=255)
                 cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACaa                                                                                                                                <  1:1719246/150‑1 (MQ=255)
                   caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAac                                                                                                                              <  2:1447718/150‑1 (MQ=255)
                    aaaaCAAGGCATGTTGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAaca                                                                                                                             <  2:1643574/150‑1 (MQ=255)
                    aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAaca                                                                                                                             >  1:1441318/1‑150 (MQ=255)
                     aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAacac                                                                                                                            >  2:1403134/1‑150 (MQ=255)
                        cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaa                                                                                                                         <  2:447855/150‑1 (MQ=255)
                          aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAAc                                                                                                                       <  1:978266/150‑1 (MQ=255)
                             cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGc                                                                                                                    <  1:144031/150‑1 (MQ=255)
                                     tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                            <  2:1441318/150‑1 (MQ=255)
                                     tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                            <  2:26095/150‑1 (MQ=255)
                                     tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                            <  1:141066/150‑1 (MQ=255)
                                           cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGc                                                                                                      <  2:1613520/150‑1 (MQ=255)
                                            ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCt                                                                                                     <  2:1057620/150‑1 (MQ=255)
                                            ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTGACCTgat                                                                                                                                           <  2:557591/112‑2 (MQ=255)
                                            ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTGACCTgat                                                                                                                                           >  1:557591/1‑111 (MQ=255)
                                             gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCtt                                                                                                    <  1:1766446/150‑1 (MQ=255)
                                              aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTTATACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTa                                                                                                   <  2:1391010/150‑1 (MQ=255)
                                              aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGTCtgatg                                                                                                                                    >  1:104254/1‑117 (MQ=255)
                                              aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGTCtgatg                                                                                                                                    <  2:104254/117‑1 (MQ=255)
                                                caacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATc                                                                                                 <  1:1754634/150‑1 (MQ=255)
                                                 aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                >  1:1993916/1‑150 (MQ=255)
                                                 aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                >  2:1687677/1‑150 (MQ=255)
                                                 aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                >  2:1680815/1‑150 (MQ=255)
                                                   caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                              <  1:188944/150‑1 (MQ=255)
                                                    aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                                             <  2:1849635/150‑1 (MQ=255)
                                                      tctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCt                                                                                           >  1:639451/1‑150 (MQ=255)
                                                              gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGtt                                                                                   >  2:692269/1‑150 (MQ=255)
                                                              gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGtt                                                                                   <  1:1269845/150‑1 (MQ=255)
                                                                  cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATg                                                                               <  2:1498967/150‑1 (MQ=255)
                                                                   aGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGc                                                                              >  2:1340665/1‑150 (MQ=255)
                                                                   aGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGc                                                                              <  2:1149147/150‑1 (MQ=255)
                                                                   aGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGc                                                                              <  1:1127226/150‑1 (MQ=255)
                                                                      ggCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAat                                                                           >  2:228096/1‑150 (MQ=255)
                                                                      ggCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAat                                                                           >  2:914474/1‑150 (MQ=255)
                                                                       gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCaaat                                                                          >  1:1208511/1‑148 (MQ=255)
                                                                          ggtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAatatat                                                                       <  2:1314034/150‑1 (MQ=255)
                                                                              gAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGaa                                                                   <  2:1001693/150‑1 (MQ=255)
                                                                                  ataCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTg                                                               <  1:1403134/150‑1 (MQ=255)
                                                                                     cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGTCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                            <  2:175147/150‑1 (MQ=255)
                                                                                     cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                            >  1:535789/1‑150 (MQ=255)
                                                                                      gcagcaGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATg                                                           <  2:639451/150‑1 (MQ=255)
                                                                                       cagcagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACaaa                                                                                                                               <  2:833368/81‑1 (MQ=255)
                                                                                       cagcagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACaaa                                                                                                                               >  1:833368/1‑81 (MQ=255)
                                                                                        agcagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGt                                                         >  2:497653/1‑150 (MQ=255)
                                                                                           agagGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCtt                                                      >  1:34346/1‑150 (MQ=255)
                                                                                              ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACTTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                   <  2:652581/150‑1 (MQ=255)
                                                                                              ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTg                                                                     >  1:1650754/1‑132 (MQ=255)
                                                                                              ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTg                                                                     <  2:1650754/132‑1 (MQ=255)
                                                                                              ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                   >  2:367743/1‑150 (MQ=255)
                                                                                              ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                   <  1:744918/150‑1 (MQ=255)
                                                                                              ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                   <  1:1549230/150‑1 (MQ=255)
                                                                                              ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                   <  1:1541899/150‑1 (MQ=255)
                                                                                                ttATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCATGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAgg                                                 <  1:845606/150‑1 (MQ=255)
                                                                                                 tATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGc                                                >  1:608800/1‑150 (MQ=255)
                                                                                                  aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGcc                                               <  1:497653/150‑1 (MQ=255)
                                                                                                  aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGcc                                               >  1:426741/1‑150 (MQ=255)
                                                                                                       tGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATa                                          >  2:803833/1‑150 (MQ=255)
                                                                                                          cgtcgtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGa                                       <  1:1687677/150‑1 (MQ=255)
                                                                                                          cgtcgtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGa                                       <  1:1680815/150‑1 (MQ=255)
                                                                                                           gtcgtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGCCTGCTGACAACCACACAACTGCCTGCTGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGAc                                      >  2:456622/1‑150 (MQ=25)
                                                                                                              gtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGc                                                              >  1:6972/1‑123 (MQ=255)
                                                                                                              gtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGc                                                              <  2:6972/123‑1 (MQ=255)
                                                                                                              gtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGtt                                   >  2:765383/1‑150 (MQ=255)
                                                                                                              gtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGtt                                   >  1:72903/1‑150 (MQ=255)
                                                                                                                aaaTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCa                                 <  2:506456/150‑1 (MQ=255)
                                                                                                                  aTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACg                               <  1:1655815/150‑1 (MQ=255)
                                                                                                                    tAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTc                             <  1:786813/150‑1 (MQ=255)
                                                                                                                     aaCGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCg                            <  1:1769629/150‑1 (MQ=255)
                                                                                                                        gTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCAt                         >  1:1657085/1‑150 (MQ=255)
                                                                                                                         tAAGTTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGcc                                               <  2:1260067/127‑1 (MQ=255)
                                                                                                                         tAAGTTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGcc                                               >  1:1260067/1‑127 (MQ=255)
                                                                                                                            ggTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCgg                     <  2:1350109/150‑1 (MQ=255)
                                                                                                                              gtctCAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCa                   <  2:989774/147‑1 (MQ=21)
                                                                                                                              gtctCAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCa                   >  1:989774/4‑150 (MQ=21)
                                                                                                                                  cAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGaa               >  1:207670/1‑150 (MQ=255)
                                                                                                                                     gTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACaa            >  1:1026389/1‑150 (MQ=255)
                                                                                                                                       cAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCATACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAg          >  2:1782575/1‑150 (MQ=255)
                                                                                                                                             tCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACt    >  2:1225962/1‑150 (MQ=255)
                                                                                                                                               aaTCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACttt  >  1:599266/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTT  >  NC_000913/3563405‑3563697

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: ‑