Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA38.7% 181.3 / 27.5 94E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (11/46);  new base A (36/0);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 2.80e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCA  >  NC_000913/3508512‑3508810
                                                                                                                                                          |                                                                                                                                                    
ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAAt                                                                                                                                                      <  1:1886951/151‑1 (MQ=255)
ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAAt                                                                                                                                                      <  2:548217/151‑1 (MQ=255)
  aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTc                                                                                                                                                    <  2:598748/151‑1 (MQ=255)
         ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCacacac                                                                                                                                              <  1:292180/150‑1 (MQ=255)
         ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCacacac                                                                                                                                              <  2:1967408/150‑1 (MQ=255)
          tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACg                                                                                                                                             <  1:53381/150‑1 (MQ=255)
           ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATAACACACggg                                                                                                                                           >  1:1581994/1‑150 (MQ=255)
           ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACgg                                                                                                                                            <  2:1367730/150‑1 (MQ=255)
            tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGGG‑CG‑G‑‑AAAAAATtaca                                                                                                                                                 >  1:1495876/1‑144 (MQ=255)
             aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTa                                                                                                                                          <  1:1042692/150‑1 (MQ=255)
             aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTa                                                                                                                                          <  1:497979/150‑1 (MQ=255)
              aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGG‑CG‑G‑‑GAAAAATCACACACGGTa                                                                                                                                          >  1:1804964/1‑149 (MQ=255)
              aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑GGAAAAAACTcacacaggcat                                                                                                                                          >  1:1144001/1‑146 (MQ=255)
                  cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑AAAAAATCACACAACGTAGccccc                                                                                                                                    >  2:429373/1‑151 (MQ=255)
                   ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATCAAAACCGGCAGcccccc                                                                                                                                   >  2:1030033/1‑150 (MQ=255)
                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGG‑CGCG‑‑AAAAATT‑ACACACCGTAGccccc                                                                                                                                    >  2:1145369/1‑145 (MQ=255)
                          tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑GG‑G‑‑GGAAAATAACACAAAGGAGCCCCcac                                                                                                                                  >  1:149635/1‑145 (MQ=255)
                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATCACACACGGTAGCCCCcac                                                                                                                                  >  2:139602/1‑144 (MQ=255)
                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATCACACACGGGAGccccccc                                                                                                                                  >  2:314506/1‑142 (MQ=255)
                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATTACACACCGTGGCCCCcac                                                                                                                                  >  1:1561835/1‑144 (MQ=255)
                             atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGG‑CG‑G‑‑GAAAAATTAAACACCGTAGCCCCCACCTCAGaaa                                                                                                                          >  2:931789/1‑150 (MQ=255)
                               atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                                                                                                       <  2:1712578/151‑1 (MQ=255)
                                tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCCG‑G‑‑AAAAATT‑ACACACCGGAGCCCCCACATCCGaaa                                                                                                                          >  1:1511031/1‑147 (MQ=255)
                                  cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAAGCGCC‑CG‑G‑‑ACAAAATTACAGACGGGAAACCCCCCATCCAAAAACgtgg                                                                                                                    >  2:1135846/1‑151 (MQ=255)
                                   aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATTACACACCGGAGCCCCCCCCACACAAAACgtgg                                                                                                                    >  1:1287007/1‑150 (MQ=255)
                                    gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtggtg                                                                                                                  <  1:1829265/151‑1 (MQ=255)
                                    gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtggtg                                                                                                                  <  2:896925/151‑1 (MQ=255)
                                     aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCG‑CG‑G‑‑GAAAAATCAAACACCGGAGCCCCCCCACCAGaaaaagggg                                                                                                                    >  2:309178/1‑142 (MQ=255)
                                      tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATCACACCCGGGAGGCCCCCCATCAAAAAGCggggtg                                                                                                                  >  2:1296505/1‑149 (MQ=255)
                                       gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                                <  1:1595256/150‑1 (MQ=255)
                                          cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                            <  1:721667/151‑1 (MQ=255)
                                           aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCG‑GG‑G‑‑GAAAAATTACAACCGGGAGCCCCCCCACCAGaaaaagggggg                                                                                                                  >  1:1614311/1‑136 (MQ=255)
                                             gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                                                                         <  2:519452/151‑1 (MQ=255)
                                             gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAAATTACACACCGGAGACCCCACCTCACAAAACGGGGGGGGCCCAcc                                                                                                          >  2:1451902/1‑150 (MQ=255)
                                              ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                        <  2:1378009/151‑1 (MQ=255)
                                                cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATCACAAACGGGAACCCCCCCATCAAAAAACGGGGGGGGCACAcc                                                                                                          >  1:501010/1‑147 (MQ=255)
                                                   gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                                   <  2:2075622/151‑1 (MQ=255)
                                                    ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGG‑CG‑G‑‑GAAAAATCACACACGGTAGCCCCCCCATCCGAAAACGGGGGGGGCACACCCAGCTGc                                                                                                   >  2:912202/1‑150 (MQ=255)
                                                     gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATTAAACACCGTAGGCCCCCCACCAGAAAAAGGGGTgcggcaacccggccgcg                                                                                                  >  2:1829959/1‑136 (MQ=255)
                                                       gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATCAAACACCGTAAGCCCCCCAACACAAAAGgtgggggggac                                                                                                             >  1:53853/1‑130 (MQ=255)
                                                          aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                            <  2:20909/151‑1 (MQ=255)
                                                            ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATCACAAACCGGAGGCCCCCCCACAGAAAACGAGGGGCGCACACCCAGCcgcg                                                                                                  >  1:1799595/1‑143 (MQ=255)
                                                              gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                        <  1:2000506/151‑1 (MQ=255)
                                                                   aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                   <  1:1215882/151‑1 (MQ=255)
                                                                     cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                                  <  2:915760/150‑1 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATCACACACCGGAACCCCCCCACCAAAAAAAGGGGGGGGCCCACCCAGCCGCGTGGGGAGCATCCCtt                                                                                   >  1:2070219/1‑148 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATCACACACCGTAGCCCCCCCATCAGAAAGCGTGGGGCGCCCACCCAGCTGCGTAGCGATCATCCCTTGc                                                                                 >  2:1706441/1‑150 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                                                                <  1:1963556/151‑1 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCG‑CG‑G‑‑GAAAAATAACAAAAGGGAGGCCCCCCAACAGAAAAAGGGGGGGGCACACCGAGCCGCGGAGCGGGCATccc                                                                                     >  1:1985647/1‑146 (MQ=255)
                                                                       aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGGCCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                               <  2:645595/151‑1 (MQ=255)
                                                                        aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATTACACACGGGAGCCCCCACATCAGAAAGGGTGGGGTGCCCACCCCGCCGCCTGGCGTTCCTCCCCTTCTg                                                                               >  1:1314388/1‑150 (MQ=255)
                                                                            cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                            <  2:90141/149‑1 (MQ=255)
                                                                             aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                         <  2:167370/151‑1 (MQ=255)
                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATCACAAACGGGAGCCCCCCCATCAGAAAGCGTGGGGGGCACACCGAGCCGCGTAGCGATCATCCCTTGCTGATaaa                                                                          >  2:1355179/1‑148 (MQ=255)
                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                        <  1:1314682/150‑1 (MQ=255)
                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                       <  1:1661270/151‑1 (MQ=255)
                                                                                 ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                                     <  1:1087100/151‑1 (MQ=255)
                                                                                   tgtCTCCTGCAAACCGGAGGGAGGGAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATTAAACAAGGGAGCCCCCCCATCAAAAAACGGGGGGCGCACACCCAGCCGCGGAGCGATCATCCCTTGCTGATaaaac                                                                        >  2:954572/1‑144 (MQ=255)
                                                                                   tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGTGGTGGGCACCCCGAGCTGCGTAGCGGTCATCCCCTGCTGGTAAAACAacac                                                                   >  2:1021740/1‑151 (MQ=255)
                                                                                    gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                                  <  2:337474/151‑1 (MQ=255)
                                                                                      ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAATTCAAACACCGGAGCCCCCCCATCAAAAAACGGGGGGCGCACACCCAGCCGCGGAGCGATCATCCCTTTCTGATAAACCAGCACTgg                                                                >  2:1059155/1‑151 (MQ=255)
                                                                                       tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGGGCG‑CG‑G‑‑GAAAAATTAAACACGGGAGCCCCCCCATCAGAAAACGGGGGGGGCCCACCCAGCTGCGGATCGACCATCCCTTGTTGGTAAAACAGCACTgg                                                                >  2:913823/1‑150 (MQ=255)
                                                                                        ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                                <  2:537706/149‑1 (MQ=255)
                                                                                        ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                              <  1:52382/151‑1 (MQ=255)
                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac                                                             <  1:868030/149‑1 (MQ=255)
                                                                                               aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATTACAAACGGTAGCCCCCACATCAAAAAACGGGGGGCGCACACCGAGCTGCCGAGCGATCATCCCCTGCTTGTAAAACAGCACCGGcaccac                                                          >  2:693886/1‑148 (MQ=255)
                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAATTTAAACACGGGAGCCCCCCCCACACAAAACGTGGGGCGCCCACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAGCACTGGCCCCACTTcgc                                                     >  2:1190583/1‑151 (MQ=255)
                                                                                                  ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATTACAAACGGGAGCCCCCACATCAGAAAACGGGGGGGGCCCACCGAGCTGCGGAGCGATCATCCCCTTCTGATAAAACAACACTGGcaccacct                                                        >  2:1197852/1‑145 (MQ=255)
                                                                                                   gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATTAAACACGGTAGCCCCCACAACCGAAAAGGGGGGGCGCACACCGAGCCGCGTAGCGATCATCCCCTTCTGATAAAACAACACCGGCCCCCCCTCCCgg                                                   >  2:579514/1‑151 (MQ=255)
                                                                                                   gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAAATAAACACCGTAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGGTCATCCCCTCCTGATAAAACAGCACTGGCCCCACCTCCCgg                                                   >  1:622118/1‑151 (MQ=255)
                                                                                                    agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑GAAAAATCACAAACCGGAGCCCCCCCATCAGAAAAGGTGGTGCGCACACCGCGCCGCGTAGCGATCCTCCCCTGCTGGTAAAACAACACTGGCACCACCTCCCggg                                                  >  1:1059348/1‑151 (MQ=255)
                                                                                                     gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                 <  2:494398/151‑1 (MQ=255)
                                                                                                      ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTa                                                <  2:198458/151‑1 (MQ=255)
                                                                                                        aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑CG‑G‑‑GAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCACCTCGCggg                                                  >  2:1584031/1‑147 (MQ=255)
                                                                                                          gaaagaaaAGAGCGACCATTTTGTGGCGGGGGGGCG‑CG‑G‑‑GAAAAATTAAACACAGTAGCCCCCCCCTCAGAAAACGGGGTGGGCACACCCAGCTGCGGTGCGATCATCCCTTGCTGGTAAAACACCACCAGGACCACTTCGCcggg                                                 >  2:859141/1‑145 (MQ=255)
                                                                                                               aaaagaGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtg                                       <  2:1561090/151‑1 (MQ=255)
                                                                                                                 aagaGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                     <  2:415559/151‑1 (MQ=255)
                                                                                                                         ccATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                             <  1:1504333/151‑1 (MQ=255)
                                                                                                                           aTTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                           <  2:1599918/151‑1 (MQ=255)
                                                                                                                           aTTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACGGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                           <  1:1347/151‑1 (MQ=255)
                                                                                                                                  ggcggtggcgCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGCATGGTGGCTGTGACCAATGGTCGGATCg                     >  1:873613/1‑150 (MQ=255)
                                                                                                                                   gcggtggcgCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                   <  2:431600/151‑1 (MQ=255)
                                                                                                                                      gtggcgCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc                >  1:218489/1‑151 (MQ=255)
                                                                                                                                        ggcgCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt              <  1:1438342/151‑1 (MQ=255)
                                                                                                                                         gcgcg‑cG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg             <  1:152647/151‑1 (MQ=255)
                                                                                                                                          cgcg‑cg‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg             >  1:1838415/1‑150 (MQ=255)
                                                                                                                                             g‑cg‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc         <  2:105432/151‑1 (MQ=255)
                                                                                                                                             g‑cg‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc         >  2:1837957/1‑151 (MQ=255)
                                                                                                                                             g‑cg‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc         <  2:95087/151‑1 (MQ=255)
                                                                                                                                               cg‑G‑‑AAAAATTCACACACGGTAGACCCCACATCAGAAAGCGTGGTGCGCACACCGAGATGCGTAGAGATAAGCCCTTGCGGAAAAAACAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg        <  1:2020660/151‑1 (MQ=255)
                                                                                                                                                g‑g‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc       >  1:145849/1‑151 (MQ=255)
                                                                                                                                                g‑g‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc       <  1:1145335/151‑1 (MQ=255)
                                                                                                                                                g‑g‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc       <  1:1060448/151‑1 (MQ=255)
                                                                                                                                                       aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc   >  1:152396/1‑151 (MQ=255)
                                                                                                                                                       aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc   >  2:363649/1‑151 (MQ=255)
                                                                                                                                                        aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCa  <  2:1605484/151‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                    
CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑CG‑G‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCA  >  NC_000913/3508512‑3508810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑