Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC40.3% 101.6 / 31.9 67T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (38/2);  new base C (0/27);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 9.89e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655
                                                                                                                                                           |                                                                                                                                                      
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAacac                                                                                                                                                       >  2:776968/1‑151 (MQ=255)
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCATGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACTGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAacac                                                                                                                                                       >  1:1090681/1‑151 (MQ=255)
gCTAACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAaca                                                                                                                                                        <  2:1006778/150‑1 (MQ=255)
 cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACAcc                                                                                                                                                      >  1:151532/1‑151 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAacac                                                                                                                                                       <  2:495256/147‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACAcc                                                                                                                                                      <  1:1031257/147‑1 (MQ=255)
      cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGc                                                                                                                                                 >  1:1153905/1‑151 (MQ=255)
      cgacgATGACTTCGCCGCGTAGCATGGAAGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGCTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCg                                                                                                                                                  >  1:548233/4‑150 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                                >  1:65894/1‑151 (MQ=255)
         cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                               >  2:61593/1‑150 (MQ=255)
          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                             <  2:290625/151‑1 (MQ=255)
                       cgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                                 >  1:24913/1‑150 (MQ=255)
                            cgggggcGGCGCCCGAAAAATGGTGAGTGTTTCCCTGCCCGTACGGGAACAGGTCGCC‑GCGCT‑TCCGGGACCGGTGT‑TGTCTCCCGACCCTTA‑AAC‑CCTC‑TAAAAGGATAATATGACAACCCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  1:299177/145‑1 (MQ=255)
                              tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                          >  1:1063888/1‑150 (MQ=255)
                               ggACGGGAACCGCAAAAGGGGGAGTTATTCCCTGCGGGTAGTGAAACGGGCCGCC‑CCGCG‑TCCCGAACCGGTGG‑TTTATCTC‑AACCGTT‑AAACCCGT‑TATGAGAATAAATTGATAACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                        <  1:526956/151‑1 (MQ=255)
                                  cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                     >  2:823480/1‑151 (MQ=255)
                                   ggCGACCGCAAAATGGGGAGTTATTCCCTGGGCGTACGGGACCGGGCCCCC‑GCGGG‑TCCCGGACCGGTGG‑TGTATCTCGACCCGTA‑AAC‑CCGC‑TATCAGGATCATATGACAAACCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  1:371312/144‑1 (MQ=255)
                                   gccagccGAAAATTGGTGAGTTATTCCCGGCGCTTCCAGAAACGGGGCGCC‑GCGCG‑GCCCGAATCCTTGT‑TTTATCTCGACCCTTA‑AAC‑CCAC‑TATCAGGATAATATAACCAACCCCCTCGCCAATGCCCGCCGTTACATCGTCGaaaa                                                                                                                      <  2:1216423/144‑1 (MQ=255)
                                    gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    >  1:83896/1‑150 (MQ=255)
                                     cGACCGAAAAAGGGGGAGTTATCCCCCGCGCGTACTGAAAAAGTTCCCC‑GCGCG‑CACCGGACCCGGGG‑TGTTTCTCGACCCGTT‑AAC‑CCGC‑TATCAGGATCATATAATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  2:433555/150‑1 (MQ=255)
                                        cgaaaaaaTGGTGAGTTTTCCCCTCCGTTTAGTGAACCGGGTCGCC‑GCGCG‑TACCGGACCAGTGT‑TGTATCTC‑AACCCTT‑AAA‑ACTCTTACCGGCAAAATAGGATCACCCCCCCCGCGGATCGCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  2:65899/143‑1 (MQ=255)
                                        ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               >  1:605015/1‑151 (MQ=255)
                                        ccGAAAAATGGTGAGTGATCCCCTGCGCTTCCTGGAACAGGCGCCC‑GCGCG‑TCCCGAAACCGTGG‑TGTATCTCGACCCTTT‑AAC‑CCGC‑TATCAGGATCATATGATCACCCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               <  2:90377/151‑1 (MQ=255)
                                        ccGAAAAAGGGGGAGTGATCCCCGGCGCTTACGGAAAAAGGTCGCC‑GCGCG‑TCCGGGACCCGGGG‑TTTCTCCCG‑ACCGTT‑AAACCCGC‑TATAAGGATCATATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                  <  2:122678/148‑1 (MQ=255)
                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                              >  1:1058964/1‑150 (MQ=255)
                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                              >  1:796330/1‑150 (MQ=255)
                                               aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                        >  1:656792/1‑151 (MQ=255)
                                                gggggagattctCCCTGCGTGT‑CTGGGAACGGTCGCC‑‑CGCGCTCCCGGACACGTGG‑TTTTTCCCGGACCCTTAAAA‑CCTC‑TAGAGGAATAATATAAAAAACCCCCTGGCCGAGACCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                          <  1:23311/137‑1 (MQ=255)
                                                 ggtgggtgTTCCCCGGCGCGTACGGGAACGGGCCGCC‑GCGCT‑TCCGGGACCGGTGT‑TGTACCCCGACCCGTA‑AAC‑CCGC‑TACCAGGTTAATATTAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                       <  1:838978/150‑1 (MQ=255)
                                                      gtgaTTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  2:636055/1‑147 (MQ=255)
                                                         ttccccTGCGC‑TTCTGGAAAGGGCCCCCGGCGCG‑TCCGGGACCGGTGT‑TTTTTCCCAAACCGTA‑AAC‑CCTC‑TATAGGGAATAAATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                                 <  1:816275/145‑1 (MQ=255)
                                                           tCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                            >  2:1157013/1‑151 (MQ=255)
                                                                 cgggtatgggAACAGGGCGCC‑GCGTG‑TCCCGGAACCGGGG‑TGTATCTCGAACTTTT‑AAC‑CCGA‑TATCAGGAACAAAGGATCAGCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          <  2:197026/139‑1 (MQ=255)
                                                                   cgTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                    >  2:32300/1‑151 (MQ=255)
                                                                     tACGGGAACGGGCCGCC‑GCGCT‑TACCGGAACGGTGG‑TGTATCTCGACCCTTA‑AAA‑CCGC‑TATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   <  2:1035057/150‑1 (MQ=255)
                                                                         gaaaCGGGCCGCC‑GCGCG‑CCCCGGACCCGTGG‑TGTACCGCAACCCGTT‑AAC‑CCGC‑TAGAAGGATCAAAGGAACACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                                  <  2:732614/145‑1 (MQ=255)
                                                                           aaCGGGCCGCC‑GCGCG‑TCCCGGAACCGTGG‑TGTATCTCGACCCTTA‑AAC‑CCGC‑TATCAGAATAATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            <  1:897183/151‑1 (MQ=255)
                                                                               ggccgcc‑gcGCG‑TA‑CGGAACCGGGGTTTTACCCCGACCCGTA‑AAC‑CCGC‑TTTAAGAATAATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          <  1:332441/146‑1 (MQ=255)
                                                                                     c‑gcgcg‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  2:747696/1‑150 (MQ=255)
                                                                                            gtcccGGACCGGTGG‑TTTATCCCAACCCGTT‑AAC‑CCGC‑TATCAGAATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                                  <  2:832293/143‑1 (MQ=255)
                                                                                               ccGGACCCGGGT‑TGTATCTCGAACCGTT‑AAC‑CCGC‑TACAAGGAACATATAATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  1:556052/149‑1 (MQ=255)
                                                                                                    aCCGTGG‑GGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATTAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGACGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                     >  2:182531/1‑151 (MQ=255)
                                                                                                     ccGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    >  1:782336/1‑151 (MQ=255)
                                                                                                       gtgg‑ttTATCTCAAACCGTA‑AAA‑CCTC‑TATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                   <  2:1032090/150‑1 (MQ=255)
                                                                                                       gtgg‑tgTACCTCGACCCTTA‑AAC‑CCGC‑TATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  2:832152/141‑1 (MQ=255)
                                                                                                         gg‑tgTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  >  2:629793/1‑149 (MQ=255)
                                                                                                         gg‑tgTATCCCGACCCGTA‑AAA‑CCGT‑TATCAGGATCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                                <  2:548602/151‑1 (MQ=255)
                                                                                                            tttATCTCGAACCATA‑AAC‑CCGC‑TATAAGAATAAAATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                               <  1:840895/148‑1 (MQ=255)
                                                                                                                tctcGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          >  2:308157/1‑151 (MQ=255)
                                                                                                                 ctcGACCCTTA‑AAC‑CCGC‑TATCAGAATCATATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                           <  2:135359/149‑1 (MQ=255)
                                                                                                                 ctcGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                         >  1:641747/1‑151 (MQ=255)
                                                                                                                    gAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAGGATAACCTGGCGCAGGCGCTGCATg                                      >  2:779667/1‑151 (MQ=255)
                                                                                                                     aaCCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      >  2:770529/1‑150 (MQ=255)
                                                                                                                           t‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                               >  1:115787/1‑151 (MQ=255)
                                                                                                                              aa‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                              >  1:1121956/1‑150 (MQ=255)
                                                                                                                               aaccctttataAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                                  <  1:1003933/140‑1 (MQ=255)
                                                                                                                                 ccGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                           >  2:5030/1‑151 (MQ=255)
                                                                                                                                    gttaccAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                          <  1:506568/145‑1 (MQ=255)
                                                                                                                                    c‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                        >  2:671020/1‑151 (MQ=255)
                                                                                                                                        tCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     >  2:825658/1‑151 (MQ=255)
                                                                                                                                         cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGAAGGCGCTGCATGACAACCGCGACCTGCTCg                    >  2:663037/1‑151 (MQ=255)
                                                                                                                                          aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                   >  1:819004/1‑151 (MQ=255)
                                                                                                                                            gatattattaacAACCCCCCCCCCGACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                    <  1:508612/137‑1 (MQ=255)
                                                                                                                                            gATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATCACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 >  1:393270/1‑151 (MQ=255)
                                                                                                                                                 tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg            >  1:259400/1‑151 (MQ=255)
                                                                                                                                                  atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  1:1081852/1‑151 (MQ=255)
                                                                                                                                                   tGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 <  2:966706/144‑1 (MQ=255)
                                                                                                                                                      tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  1:516238/1‑151 (MQ=255)
                                                                                                                                                        aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  1:467471/1‑151 (MQ=255)
                                                                                                                                                        aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  2:470363/1‑151 (MQ=255)
                                                                                                                                                        aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  2:400032/1‑151 (MQ=255)
                                                                                                                                                           acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:792602/1‑151 (MQ=255)
                                                                                                                                                           acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:1072186/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                      
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCG‑TACCGGAACCGTGG‑TGTATCTCGAACCGTT‑AAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655

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