Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC45.1% 171.0 / 21.4 103T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (42/14);  new base C (0/46);  total (42/61)
Fisher's exact test for biased strand distribution p-value = 8.88e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648
                                                                                                                                                            |                                                                                                                                               
cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAa‑ca                                                                                                                                                 <  1:886433/151‑1 (MQ=255)
    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CAcc                                                                                                                                               <  1:35025/149‑1 (MQ=255)
    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CAcc                                                                                                                                               <  1:1419993/149‑1 (MQ=255)
     aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTc                                                                                                                                            >  1:8086/1‑151 (MQ=255)
       cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCg                                                                                                                                           <  4:63239/150‑1 (MQ=255)
       cgccgACGACCTCGCCGCGCAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGc                                                                                                                                          >  4:796688/1‑151 (MQ=255)
        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGcc                                                                                                                                         <  2:122441/151‑1 (MQ=255)
        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGcc                                                                                                                                         <  1:222564/151‑1 (MQ=255)
          cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCg                                                                                                                                        >  1:866986/1‑150 (MQ=255)
          cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCg                                                                                                                                        >  4:1432106/1‑150 (MQ=255)
          cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGa                                                                                                                                       <  1:123015/151‑1 (MQ=255)
           gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGAt                                                                                                                                      >  1:778423/1‑151 (MQ=255)
           gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGAt                                                                                                                                      >  1:26883/1‑151 (MQ=255)
            aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CTCCCTCGCCGa                                                                                                                                       <  4:628276/149‑1 (MQ=255)
            aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATg                                                                                                                                     <  1:882336/151‑1 (MQ=255)
             tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGc                                                                                                                                    >  4:708439/1‑151 (MQ=255)
                ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCccg                                                                                                                                 <  4:115248/151‑1 (MQ=255)
                  tCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCccgcc                                                                                                                               >  4:1331195/1‑151 (MQ=255)
                   cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCccgccg                                                                                                                              >  3:153831/1‑151 (MQ=255)
                     ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGGTCATATGATCAA‑CACCCTCGCCGATGCCCGCCGtt                                                                                                                            >  1:1262101/1‑151 (MQ=255)
                     ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGtt                                                                                                                            >  4:1413760/1‑151 (MQ=255)
                         gTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACAt                                                                                                                        >  2:43962/1‑151 (MQ=255)
                             ccgggccGGCGCCCGCAAAAGGGTGAGTGTTCCCCTGCGCGTACGGG‑AACAGG‑TCCCCGCGGGTACCGGAACG‑GGGGTGTATCTCGACCGTTTAAC‑CCGT‑TATCAGGATAATATGATCAC‑CCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                    <  1:135150/145‑1 (MQ=255)
                             ccgggACGGCGCCCGAAAAATGGTGAGTTTTCCCCTGGGCTTACGGA‑AACGGG‑TCGCCGCGCGTACGGGACCG‑GTGGTGTATCCCGAACCTTTAAC‑CCGC‑TATCAGGATAATATAATCAA‑CCCCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                      <  3:80993/146‑1 (MQ=255)
                               tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                   >  1:1448415/1‑150 (MQ=255)
                                    ggCGCCCGAAAAATGGTGAGTGTTTCCCGGCGGGTACGGG‑AACAGG‑TCGCCGCGCTTACCGGACCG‑GTGGTGTATCTCGAACCGTTAAC‑CCCC‑TATCAGAATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  1:672941/150‑1 (MQ=255)
                                    ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             >  1:799116/1‑151 (MQ=255)
                                    ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             >  1:1126305/1‑151 (MQ=255)
                                       gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                          >  1:1290112/1‑151 (MQ=255)
                                         ccGCAAATTGGTGGGTGTTTCCCGGCGCGACCGGA‑AACGGG‑CCGCCGCGCGTCCGGCAGCC‑GTGTTGTATCTCGACCGGTAAAC‑CCGT‑TTCCGGGATAATATAACCAC‑CCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             <  1:144382/146‑1 (MQ=255)
                                         ccGAAAAGTGGTGGTTGATCCCCCGCGCGTGCTGG‑AACGGC‑CCGCCGCGC‑CCCCGGAACCTGTGTTTTCTCTCGAACCGTTAAC‑CCCC‑TACAAGTATCATATGACCAC‑CCCCCCCGCCGACGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                         <  4:1008162/150‑1 (MQ=255)
                                         ccGAAAAGTGGGGAGTGTTCCCCGGCGCGTATGGG‑AACAGG‑CCGCCGCGCGT‑CCGGAACCGGTGGTTTATCCCGACCCTTTAAC‑CCGC‑TATAAGAATAATATAACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                          <  4:1450770/149‑1 (MQ=255)
                                             aaaggggTGAGTGTTCCCCTGCGCGTATGGG‑AACGGG‑CCCCCGCGGGTCCCGAACCG‑GTGTTTTTTCTCG‑ACCCTTAAC‑CCCCTTATAAGGATAATATGAACAC‑CCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        <  4:1054550/142‑1 (MQ=255)
                                             aaaaTGGTGAGTGTTTCCCTGCGCGTACTGG‑AACGGG‑TCGCCGCGCGTCCCGGAACC‑GTGGTGTATCTCGACCCGTTAAC‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                     <  1:379518/150‑1 (MQ=255)
                                             aaaaTGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                    >  1:1076199/1‑151 (MQ=255)
                                             aaaaGGGTGAGTGATTCCCGGCGCGTACGGG‑AACAGT‑TCCCCGCGCGTACCGGACCC‑GTGTTTTATCTCGACCCGTTAAA‑CCGC‑TACCAGGATCATATAATCAC‑CCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                         <  1:702327/146‑1 (MQ=255)
                                                  ggTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  1:660753/1‑150 (MQ=255)
                                                   gtgagtgtTCCCCTGCGC‑TTCCGG‑AACCGGCTCGCCGCGCGTACGGGACCG‑GGGGTGTCTCCCGACCCTTAAAC‑CCGT‑TATCAGGTTAATATGAACAC‑CACCCCCCCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                <  4:14769/149‑1 (MQ=255)
                                                     gactgaTCCCCTGCGAGT‑CTGGAAACCGG‑TCCCGGCGCGTCCGGGCCCG‑GTGGTGTTTCTCCACCCATAAAC‑CCCC‑TCAAAGGACAAAAACATAAC‑CACCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                            <  4:1159373/148‑1 (MQ=255)
                                                       gtgaTCCCGGGCGCGTACGGA‑AACAGG‑CCGCCGCGGGTACCGGAGCG‑TGGGTTTATCTCG‑ACCCTTAAA‑CCACTTACGAGGATAATATGAAAAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                           <  4:1413378/150‑1 (MQ=255)
                                                         gaTTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTACCTCGACCCGTAAAC‑CCGC‑TATCAGGATAATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         <  4:180573/150‑1 (MQ=255)
                                                           ttCCCTGCGCGTACGGG‑AACAGG‑TCGCCCCGCGTCCCGAACCG‑GTGGTGTATCCCAACCCGTTAAC‑CCGC‑TAGAAGGTTCATATGATCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       <  1:1266576/150‑1 (MQ=255)
                                                           tccccTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGTTGTATCTCGACCCGTAAAC‑CCGC‑TATCAGGATCATATGAACAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       <  1:458672/148‑1 (MQ=255)
                                                           tcccaggcgtGTATTGG‑AACAGT‑TCGCCGCGGGTGCCGAAACC‑GGGGTTTCTCTGTAACCTTTAAA‑CCGC‑TATAAGGACCATATTATAAA‑CCCCCTCCCTGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       <  1:605076/144‑1 (MQ=255)
                                                             cccGGCGCTTACTGG‑AACAGG‑CCGCCGCGCGTA‑CGGAACCGGTGGTTTATCTCGACCCGTTAAC‑CCGC‑TATCGGAACAATATGTTAAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                     <  3:2704/150‑1 (MQ=255)
                                                              ccggcgcTTACTGA‑AACAGT‑TCGCCGGGCGTACCGAACCG‑GTGGTTTATCCCGACCCTTTAAC‑CCGT‑TATAAGGATCATATGATCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                     <  4:1213938/146‑1 (MQ=255)
                                                              ccTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                   >  1:765033/1‑151 (MQ=255)
                                                               cTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  4:856564/1‑151 (MQ=255)
                                                                tGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  4:1392283/1‑150 (MQ=255)
                                                                tGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  4:626543/1‑150 (MQ=255)
                                                                       aCTGA‑AACAGT‑TCGCCGCGCGTACCGGACCC‑GTGGTTTATCTCGACCCGTAAAC‑CCGC‑TATCAGAATCATATGACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  <  1:1124202/143‑1 (MQ=255)
                                                                         acggaaacggttCGCCGCGCGTCCCGGACCG‑GTGGTGTACCTCGACCCGTAAAC‑CCGC‑TATCAGGATCATATGATCAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              <  1:697395/136‑1 (MQ=255)
                                                                          gg‑AACAGG‑TCGCCGCGCGTCCCGGACCC‑GGGGTGTATCTCGACCCGTAAAC‑CCGC‑TATCAGGATAATATGATAAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                        <  1:1098440/150‑1 (MQ=255)
                                                                          gg‑AACAGG‑TCGCCGCGCGTCCCGGAACG‑GTGGTGTATCTCGACCCTTTAAC‑CCGT‑TATGAGAATCATATGACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                   <  2:58505/139‑1 (MQ=255)
                                                                             aaCAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  1:1471459/150‑1 (MQ=255)
                                                                               cAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                    >  4:441091/1‑150 (MQ=255)
                                                                                 gg‑TCGCCGCGCGTACGGGAACG‑GTGGTGTTTCTCGACCCGTTAAC‑CCGC‑TATCAGGATCATATGATCAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                    <  1:357882/148‑1 (MQ=255)
                                                                                 gg‑TCGCCGCGCGTACCGGACCC‑GTGGTGTATCTCGACCCGTTAAA‑CCGC‑TATCAGGATCATATAACAAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                  <  1:148675/150‑1 (MQ=255)
                                                                                 gg‑TCGCCGCCCGTCCCGGAACC‑GTGGTTTATCGC‑AACCCTTAAACCCGT‑TACCAGGATCATATGACCAC‑CCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                     <  2:39895/147‑1 (MQ=255)
                                                                                  g‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                 >  1:563085/1‑150 (MQ=255)
                                                                                       ccgcGGGTCCGGGACCC‑GTGGTGTATCTC‑AACCCTTAAACCCGC‑TATCAGGATAATATGATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTACCACATGAATATCga                                                            <  1:126235/151‑1 (MQ=255)
                                                                                       ccgcGCGTCCCGGACCG‑GTGGTGTATCTCGAACCGTTAAC‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               <  1:1272265/148‑1 (MQ=255)
                                                                                       ccgcGCGACCCGAACCC‑GGGTTGTATCACGACCCTTAAAA‑CCGC‑TACCAGAATTATATGACCAC‑CCCCCGCCCCTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   <  4:828990/144‑1 (MQ=255)
                                                                                         gcgcgTCCGGGAACC‑GTGGTTTATCTCGACCCTTAAAC‑CCGC‑TATCGGAATCATATGATCAA‑CCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                  <  4:275373/143‑1 (MQ=255)
                                                                                            gcgtccGGAACCGGTGGTTTATCTCGAACTTTAAAC‑CCGC‑TATCAGAATAATATGACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                           <  1:1108033/144‑1 (MQ=255)
                                                                                            cgTCCCGGAACC‑GTGGTTTATCTCGACCCGTAAAA‑CCGC‑TATCAGGATAATATGACCAC‑CACCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            <  4:497289/146‑1 (MQ=255)
                                                                                                ggtttaCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                    >  4:469271/6‑150 (MQ=255)
                                                                                                ccGGACCC‑GTGTTGTATCTCGACCCGTTAAC‑CCGC‑TATAAGTATCATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                       <  1:1003220/147‑1 (MQ=255)
                                                                                                   ggaccgggggtgTATCTC‑AACCCTTAAACCCGC‑TATCAGGATCATATGATAAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                     <  1:725465/139‑1 (MQ=255)
                                                                                                       cagtggtgTATCCCGACCCCTAAAA‑CCGA‑TATCAGGCTTAT‑TTATCAACCCCCCCCTCCCACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                               <  4:532149/147‑1 (MQ=255)
                                                                                                           ggtgTATCTCGAACCGTTAAC‑CCGC‑TATCAGAATAATATGATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           <  1:916661/149‑1 (MQ=255)
                                                                                                           ggtgTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  2:148334/1‑149 (MQ=255)
                                                                                                             tgtATCTCGACCCGTAAAC‑CCGT‑TACCAGGATATTATACTAAA‑CCCCCCCGCGGAGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                            <  1:1296927/146‑1 (MQ=255)
                                                                                                             tgtATCTCGAACCTTTAAC‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                        <  2:11252/150‑1 (MQ=255)
                                                                                                              gtATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                      >  4:632113/1‑151 (MQ=255)
                                                                                                                aTCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    >  1:1436626/1‑151 (MQ=255)
                                                                                                                     gAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               >  4:661666/1‑151 (MQ=255)
                                                                                                                     gAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               >  4:515949/1‑151 (MQ=255)
                                                                                                                      acccGTAAAC‑CCGC‑TATCAGGATCATATAATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                   <  2:26767/144‑1 (MQ=255)
                                                                                                                       aCCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                             >  1:1366110/1‑151 (MQ=255)
                                                                                                                        cctttAAC‑CCGT‑TACCGGGATAATATTACAAA‑CCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                            <  1:432812/148‑1 (MQ=255)
                                                                                                                        ccGTAAAC‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                             <  1:535254/150‑1 (MQ=255)
                                                                                                                         cGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                            >  1:1441824/1‑150 (MQ=255)
                                                                                                                           ttAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                         >  4:575770/1‑151 (MQ=255)
                                                                                                                            tAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                        >  3:28123/1‑151 (MQ=255)
                                                                                                                             taaacccgcttaaaGGATCATATAACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                             <  4:1308944/134‑1 (MQ=255)
                                                                                                                             aaa‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                       >  1:1290768/1‑151 (MQ=255)
                                                                                                                              aacccGC‑TATCAGAATAATATAACCAC‑CCCCCTCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                        <  1:417759/147‑1 (MQ=255)
                                                                                                                              aacccGC‑TAACAGAATCATAGGATCAC‑CCCCCTCCCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                          <  1:1068730/145‑1 (MQ=255)
                                                                                                                                acccgttagatggATTATATAACCAA‑CCCCCCCTGCTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                          <  1:1012290/136‑1 (MQ=255)
                                                                                                                                 ccGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                    >  1:1376821/1‑151 (MQ=255)
                                                                                                                                    c‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 >  2:10488/1‑151 (MQ=255)
                                                                                                                                       acgagaataatatGACCAA‑CCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt               <  4:170124/142‑1 (MQ=255)
                                                                                                                                        tCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc              >  1:318048/1‑151 (MQ=255)
                                                                                                                                         cAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGcgcg             >  1:1458071/1‑151 (MQ=255)
                                                                                                                                           ggttaatatGATCAC‑CACCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc              <  1:904180/143‑1 (MQ=255)
                                                                                                                                           ggATCATATGATCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt               <  1:1380436/147‑1 (MQ=255)
                                                                                                                                                   tGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc    >  2:21304/1‑150 (MQ=255)
                                                                                                                                                     aTCAC‑CCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt      <  1:1363756/146‑1 (MQ=255)
                                                                                                                                                     aTCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  4:85354/1‑150 (MQ=255)
                                                                                                                                                            |                                                                                                                                               
CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGG‑TCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648

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