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→GATyhfWphosphopentomutase‑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

CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCA  >  NC_000913/3508512‑3508810
                                                                                                                                                          |                                                                                                                                                    
ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAAt                                                                                                                                                      <  1:1886951/151‑1 (MQ=255)
ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAAt                                                                                                                                                      <  2:548217/151‑1 (MQ=255)
  aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTc                                                                                                                                                    <  2:598748/151‑1 (MQ=255)
         ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCacacac                                                                                                                                              <  1:292180/150‑1 (MQ=255)
         ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCacacac                                                                                                                                              <  2:1967408/150‑1 (MQ=255)
          tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACg                                                                                                                                             <  1:53381/150‑1 (MQ=255)
           ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATAACACACggg                                                                                                                                           >  1:1581994/1‑150 (MQ=255)
           ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACgg                                                                                                                                            <  2:1367730/150‑1 (MQ=255)
            tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGGGCGG‑‑AAAAAATtaca                                                                                                                                                 >  1:1495876/1‑144 (MQ=255)
             aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTa                                                                                                                                          <  1:1042692/150‑1 (MQ=255)
             aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTa                                                                                                                                          <  1:497979/150‑1 (MQ=255)
              aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGG‑‑GAAAAATCACACACGGTa                                                                                                                                          >  1:1804964/1‑149 (MQ=255)
              aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAAACTcacacaggcat                                                                                                                                          >  1:1144001/1‑146 (MQ=255)
                  cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑AAAAAATCACACAACGTAGccccc                                                                                                                                    >  2:429373/1‑151 (MQ=255)
                   ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATCAAAACCGGCAGcccccc                                                                                                                                   >  2:1030033/1‑150 (MQ=255)
                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGCGCG‑‑AAAAATTACACACCGTAGccccc                                                                                                                                    >  2:1145369/1‑145 (MQ=255)
                          tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGG‑‑GGAAAATAACACAAAGGAGCCCCcac                                                                                                                                  >  1:149635/1‑145 (MQ=255)
                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATCACACACGGTAGCCCCcac                                                                                                                                  >  2:139602/1‑144 (MQ=255)
                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATCACACACGGGAGccccccc                                                                                                                                  >  2:314506/1‑142 (MQ=255)
                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATTACACACCGTGGCCCCcac                                                                                                                                  >  1:1561835/1‑144 (MQ=255)
                             atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGG‑‑GAAAAATTAAACACCGTAGCCCCCACCTCAGaaa                                                                                                                          >  2:931789/1‑150 (MQ=255)
                               atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                                                                                                       <  2:1712578/151‑1 (MQ=255)
                                tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCCGG‑‑AAAAATTACACACCGGAGCCCCCACATCCGaaa                                                                                                                          >  1:1511031/1‑147 (MQ=255)
                                  cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAAGCGCCCGG‑‑ACAAAATTACAGACGGGAAACCCCCCATCCAAAAACgtgg                                                                                                                    >  2:1135846/1‑151 (MQ=255)
                                   aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATTACACACCGGAGCCCCCCCCACACAAAACgtgg                                                                                                                    >  1:1287007/1‑150 (MQ=255)
                                    gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtggtg                                                                                                                  <  1:1829265/151‑1 (MQ=255)
                                    gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtggtg                                                                                                                  <  2:896925/151‑1 (MQ=255)
                                     aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGCGG‑‑GAAAAATCAAACACCGGAGCCCCCCCACCAGaaaaagggg                                                                                                                    >  2:309178/1‑142 (MQ=255)
                                      tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATCACACCCGGGAGGCCCCCCATCAAAAAGCggggtg                                                                                                                  >  2:1296505/1‑149 (MQ=255)
                                       gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                                <  1:1595256/150‑1 (MQ=255)
                                          cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                            <  1:721667/151‑1 (MQ=255)
                                           aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGGGG‑‑GAAAAATTACAACCGGGAGCCCCCCCACCAGaaaaagggggg                                                                                                                  >  1:1614311/1‑136 (MQ=255)
                                             gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                                                                         <  2:519452/151‑1 (MQ=255)
                                             gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAAATTACACACCGGAGACCCCACCTCACAAAACGGGGGGGGCCCAcc                                                                                                          >  2:1451902/1‑150 (MQ=255)
                                              ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                        <  2:1378009/151‑1 (MQ=255)
                                                cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATCACAAACGGGAACCCCCCCATCAAAAAACGGGGGGGGCACAcc                                                                                                          >  1:501010/1‑147 (MQ=255)
                                                   gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                                   <  2:2075622/151‑1 (MQ=255)
                                                    ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGG‑‑GAAAAATCACACACGGTAGCCCCCCCATCCGAAAACGGGGGGGGCACACCCAGCTGc                                                                                                   >  2:912202/1‑150 (MQ=255)
                                                     gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATTAAACACCGTAGGCCCCCCACCAGAAAAAGGGGTgcggcaacccggccgcg                                                                                                  >  2:1829959/1‑136 (MQ=255)
                                                       gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATCAAACACCGTAAGCCCCCCAACACAAAAGgtgggggggac                                                                                                             >  1:53853/1‑130 (MQ=255)
                                                          aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                            <  2:20909/151‑1 (MQ=255)
                                                            ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATCACAAACCGGAGGCCCCCCCACAGAAAACGAGGGGCGCACACCCAGCcgcg                                                                                                  >  1:1799595/1‑143 (MQ=255)
                                                              gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                        <  1:2000506/151‑1 (MQ=255)
                                                                   aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                   <  1:1215882/151‑1 (MQ=255)
                                                                     cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                                  <  2:915760/150‑1 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATCACACACCGGAACCCCCCCACCAAAAAAAGGGGGGGGCCCACCCAGCCGCGTGGGGAGCATCCCtt                                                                                   >  1:2070219/1‑148 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATCACACACCGTAGCCCCCCCATCAGAAAGCGTGGGGCGCCCACCCAGCTGCGTAGCGATCATCCCTTGc                                                                                 >  2:1706441/1‑150 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                                                                <  1:1963556/151‑1 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGG‑‑GAAAAATAACAAAAGGGAGGCCCCCCAACAGAAAAAGGGGGGGGCACACCGAGCCGCGGAGCGGGCATccc                                                                                     >  1:1985647/1‑146 (MQ=255)
                                                                       aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGGCCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                               <  2:645595/151‑1 (MQ=255)
                                                                        aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATTACACACGGGAGCCCCCACATCAGAAAGGGTGGGGTGCCCACCCCGCCGCCTGGCGTTCCTCCCCTTCTg                                                                               >  1:1314388/1‑150 (MQ=255)
                                                                            cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                            <  2:90141/149‑1 (MQ=255)
                                                                             aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                         <  2:167370/151‑1 (MQ=255)
                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATCACAAACGGGAGCCCCCCCATCAGAAAGCGTGGGGGGCACACCGAGCCGCGTAGCGATCATCCCTTGCTGATaaa                                                                          >  2:1355179/1‑148 (MQ=255)
                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                        <  1:1314682/150‑1 (MQ=255)
                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                       <  1:1661270/151‑1 (MQ=255)
                                                                                 ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                                     <  1:1087100/151‑1 (MQ=255)
                                                                                   tgtCTCCTGCAAACCGGAGGGAGGGAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATTAAACAAGGGAGCCCCCCCATCAAAAAACGGGGGGCGCACACCCAGCCGCGGAGCGATCATCCCTTGCTGATaaaac                                                                        >  2:954572/1‑144 (MQ=255)
                                                                                   tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGTGGTGGGCACCCCGAGCTGCGTAGCGGTCATCCCCTGCTGGTAAAACAacac                                                                   >  2:1021740/1‑151 (MQ=255)
                                                                                    gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                                  <  2:337474/151‑1 (MQ=255)
                                                                                      ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAATTCAAACACCGGAGCCCCCCCATCAAAAAACGGGGGGCGCACACCCAGCCGCGGAGCGATCATCCCTTTCTGATAAACCAGCACTgg                                                                >  2:1059155/1‑151 (MQ=255)
                                                                                       tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGGGCGCGG‑‑GAAAAATTAAACACGGGAGCCCCCCCATCAGAAAACGGGGGGGGCCCACCCAGCTGCGGATCGACCATCCCTTGTTGGTAAAACAGCACTgg                                                                >  2:913823/1‑150 (MQ=255)
                                                                                        ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                                <  2:537706/149‑1 (MQ=255)
                                                                                        ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                              <  1:52382/151‑1 (MQ=255)
                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac                                                             <  1:868030/149‑1 (MQ=255)
                                                                                               aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATTACAAACGGTAGCCCCCACATCAAAAAACGGGGGGCGCACACCGAGCTGCCGAGCGATCATCCCCTGCTTGTAAAACAGCACCGGcaccac                                                          >  2:693886/1‑148 (MQ=255)
                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAATTTAAACACGGGAGCCCCCCCCACACAAAACGTGGGGCGCCCACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAGCACTGGCCCCACTTcgc                                                     >  2:1190583/1‑151 (MQ=255)
                                                                                                  ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATTACAAACGGGAGCCCCCACATCAGAAAACGGGGGGGGCCCACCGAGCTGCGGAGCGATCATCCCCTTCTGATAAAACAACACTGGcaccacct                                                        >  2:1197852/1‑145 (MQ=255)
                                                                                                   gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATTAAACACGGTAGCCCCCACAACCGAAAAGGGGGGGCGCACACCGAGCCGCGTAGCGATCATCCCCTTCTGATAAAACAACACCGGCCCCCCCTCCCgg                                                   >  2:579514/1‑151 (MQ=255)
                                                                                                   gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAAATAAACACCGTAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGGTCATCCCCTCCTGATAAAACAGCACTGGCCCCACCTCCCgg                                                   >  1:622118/1‑151 (MQ=255)
                                                                                                    agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑GAAAAATCACAAACCGGAGCCCCCCCATCAGAAAAGGTGGTGCGCACACCGCGCCGCGTAGCGATCCTCCCCTGCTGGTAAAACAACACTGGCACCACCTCCCggg                                                  >  1:1059348/1‑151 (MQ=255)
                                                                                                     gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                 <  2:494398/151‑1 (MQ=255)
                                                                                                      ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTa                                                <  2:198458/151‑1 (MQ=255)
                                                                                                        aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGG‑‑GAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCACCTCGCggg                                                  >  2:1584031/1‑147 (MQ=255)
                                                                                                          gaaagaaaAGAGCGACCATTTTGTGGCGGGGGGGCGCGG‑‑GAAAAATTAAACACAGTAGCCCCCCCCTCAGAAAACGGGGTGGGCACACCCAGCTGCGGTGCGATCATCCCTTGCTGGTAAAACACCACCAGGACCACTTCGCcggg                                                 >  2:859141/1‑145 (MQ=255)
                                                                                                               aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtg                                       <  2:1561090/151‑1 (MQ=255)
                                                                                                                 aagaGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                     <  2:415559/151‑1 (MQ=255)
                                                                                                                         ccATTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                             <  1:1504333/151‑1 (MQ=255)
                                                                                                                           aTTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                           <  2:1599918/151‑1 (MQ=255)
                                                                                                                           aTTTTGTGGCGGTGGCGCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACGGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                           <  1:1347/151‑1 (MQ=255)
                                                                                                                                  ggcggtggcgCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGCATGGTGGCTGTGACCAATGGTCGGATCg                     >  1:873613/1‑150 (MQ=255)
                                                                                                                                   gcggtggcgCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                   <  2:431600/151‑1 (MQ=255)
                                                                                                                                      gtggcgCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc                >  1:218489/1‑151 (MQ=255)
                                                                                                                                        ggcgCGCGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt              <  1:1438342/151‑1 (MQ=255)
                                                                                                                                         gcgcgcGG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg             <  1:152647/151‑1 (MQ=255)
                                                                                                                                          cgcgcgG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg             >  1:1838415/1‑150 (MQ=255)
                                                                                                                                             gcgG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc         <  2:105432/151‑1 (MQ=255)
                                                                                                                                             gcgG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc         >  2:1837957/1‑151 (MQ=255)
                                                                                                                                             gcgG‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc         <  2:95087/151‑1 (MQ=255)
                                                                                                                                               cgG‑‑AAAAATTCACACACGGTAGACCCCACATCAGAAAGCGTGGTGCGCACACCGAGATGCGTAGAGATAAGCCCTTGCGGAAAAAACAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg        <  1:2020660/151‑1 (MQ=255)
                                                                                                                                                gg‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc       >  1:145849/1‑151 (MQ=255)
                                                                                                                                                gg‑‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc       <  1:1145335/151‑1 (MQ=255)
                                                                                                                                                gg‑‑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)
                                                                                                                                                          |                                                                                                                                                    
CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGG‑‑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: