Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC50.0% 150.0 / 27.4 94T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (42/5);  new base C (0/47);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 3.20e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376356‑1376652
                                                                                                                                                          |                                                                                                                                                   
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAacac                                                                                                                                                    >  2:787137/1‑151 (MQ=255)
  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAacac                                                                                                                                                    <  2:1573761/149‑1 (MQ=255)
  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACAccc                                                                                                                                                  >  2:1304270/1‑151 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACAccc                                                                                                                                                  <  1:1670152/149‑1 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCt                                                                                                                                                 <  2:1662228/150‑1 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTc                                                                                                                                                >  1:295947/1‑151 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTc                                                                                                                                                >  2:1071948/1‑151 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTc                                                                                                                                                >  2:378570/1‑151 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACA‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACCCCCTCGcc                                                                                                                                             <  1:52909/151‑1 (MQ=255)
        ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             <  1:1149098/150‑1 (MQ=255)
          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACCCCCTCGCCGAt                                                                                                                                          <  1:481999/151‑1 (MQ=255)
               ccgcgccggGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAGCAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                                     <  1:1078243/148‑1 (MQ=255)
                    ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                >  1:1429214/1‑151 (MQ=255)
                    ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                >  2:1184814/1‑151 (MQ=255)
                        gTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAt                                                                                                                            >  1:1597797/1‑151 (MQ=255)
                            ccgggACGGAGACCGAAAAATGGTGAGTGCTTCCCTGCGCGTACGGAAAAAGGTCGCCGCGCGTCCCGGACCC‑GTGT‑TGTATCTCGAACCGTTAAC‑CCGT‑T‑ATCAGGATAATATGCCCACCCCCCTCGCCGATGCCCGCCGTTACAtcgt                                                                                                                         <  1:1609499/147‑1 (MQ=255)
                            ccgatcgggCGCCCACAAAATGGTGAGTGATCCCCGGCGCGTACGGGAACAGGCCCCCGCGCGTCCCGGACCG‑GTGT‑TTTATCTCAACCCGTAAAC‑CCGT‑T‑ATCAGAATAATATAATAACCCCCCCCGCCGATGCCCGCCGTTACAtcgt                                                                                                                         <  1:53292/143‑1 (MQ=255)
                            ccTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                        >  1:650641/1‑151 (MQ=255)
                             cTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                       >  2:801114/1‑151 (MQ=255)
                               gggcGGCGCCCGTAAATTGGTGAGTGTTTCCCTGCGCGTAGTAGAACGGCTCGCCGCGCTCCCGGGAACG‑GGGG‑TGTTTCTC‑ACCCTTTAAACCCGC‑T‑GTCGGGGTTATTAGACCAACCCCCCCCCCGATGCCCGCCGTTACAtcgtcg                                                                                                                       <  1:4438/146‑1 (MQ=255)
                               ggccGGCGCCCGAAAAATGGTGAGTGTTCCCCTGCGCGTACGGGAACAGGTCGCCGCGGGTCCCGGACCC‑GTGG‑TGTATCTC‑AACCTTTAAACCCGC‑T‑ACCGGGATCATATAATAACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGaa                                                                                                                     <  2:1122675/148‑1 (MQ=255)
                                 aCGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                    >  1:764692/1‑150 (MQ=255)
                                   ggCGCCCGAAAAATGGTGGGTGATTCCCGGCGCGTACGGGAACAGGCCGCCGCGCGTCCCGAACCC‑GTGG‑TGTTTCTCGAACCTTTAAC‑CCTC‑T‑ATAAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                    <  1:629778/148‑1 (MQ=255)
                                   ggCGCCCGAAAAATGGTGAGTGTTCCCCTCCGCTTATGGAAACGGCTCGCCGCGCGTCCCGAAACC‑GTGG‑TGTATCCCGACCCTTAAAC‑CCCC‑T‑ACCGGAATAATATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                  <  1:701612/150‑1 (MQ=255)
                                   ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgt                                                                                                                         >  2:53297/1‑143 (MQ=255)
                                     cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               >  1:951196/1‑151 (MQ=255)
                                        ccGCAAAAGGGTGAGTGTTCCCCTGCGCGTCCAGAAACAGGCCGCCGCGCGTCCCGGACCC‑GTGT‑TGTATCCCGAACCTTAAAA‑CCGC‑T‑ATCAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                            <  1:1576115/151‑1 (MQ=255)
                                        ccGAAAAATGGTGAGTGATCCCCGGCGCGTACGGGAACAGGTCGCCGCGCGTACCGAACCG‑GGGG‑TGTTTCTCAACCCGTTAAC‑CCGT‑T‑ACCAGGATCATATAACCAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              <  1:473398/149‑1 (MQ=255)
                                            aaaaggtggagtggTCCCCCGCGCTTGCTGAAAGGGGTCGCCGCGCGT‑CCGGAACCGGTGG‑TGTATCGCGACCCATGAAA‑CCGC‑T‑ATCAGGTTCATAAATCCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                        <  2:1657630/143‑1 (MQ=255)
                                            aaaaGGGTGAGTGTTTCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCC‑GTGG‑TGTATCTCAACCCGTAAAC‑CCGC‑T‑ATAAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                          <  2:1444392/149‑1 (MQ=255)
                                              aagggTGAGTGTTTCCCTGCGCGTACGGGAACAGGCCCCCGCGCGTCCCGGACCC‑GTGG‑TGTATCTCGACCCTTAAAC‑CCGC‑T‑ACCAGGATAAAATGATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                      <  1:218143/148‑1 (MQ=255)
                                                  gtgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  2:371304/1‑151 (MQ=255)
                                                         tttcctggggctgtcTGGAACGGGCCGCCGGGCGTCCGGGACCT‑GGGG‑TGT‑TCCCAAACCCTTTAA‑CCCCTT‑ATGGGGAAAATATAACCAACCCCCCCGCCTATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                             <  1:1052277/135‑1 (MQ=255)
                                                          ttCCCGGCGCGTACTGAAACGGGTCGCCGCGCGTACCGGACCC‑GTGG‑TTTATCTCGACCCTTTAAC‑CCGT‑T‑ATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                           <  2:585762/150‑1 (MQ=255)
                                                             ccggcgcTTCCTGGACCAGGTCCCCGCGCTTCCGGGAGCT‑GTGG‑TTTATCACGACCCTTTAAC‑CCTC‑T‑ATCAGGATCATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                            <  1:220576/143‑1 (MQ=255)
                                                             ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                       >  2:1038420/1‑151 (MQ=255)
                                                              cgggccgtatgaaaccAGGTCGCCACGCGTCCCGAACCT‑GTGG‑ATTATCCCACCCCGTAAAA‑CCGA‑T‑CACAGGATAATATTACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                       <  2:1280095/135‑1 (MQ=255)
                                                                 cgcgagcgggacaaGGCCGCGGCGCGTGCCGAAACT‑GTGT‑TTTATCTCGACCCGTTAAA‑CCGC‑T‑ATAAGGATAATATGAAAACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                     <  1:49518/136‑1 (MQ=255)
                                                                 cgcgTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                   >  2:1244545/1‑151 (MQ=255)
                                                                     tACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                >  2:672100/1‑150 (MQ=255)
                                                                           aaCGGGCCGCCGCGCGTCCCGGACCC‑GGGG‑TGTATCCCAACCCGTAAAC‑CCCC‑T‑ATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                          <  1:546219/150‑1 (MQ=255)
                                                                            acgggTCCCGGCGCGT‑CCGGAACC‑CTGGTTTTATCCCGACCCATAAAC‑CCGT‑T‑ATGAGAAAAATTTAATCAACCCCCCCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                           <  1:1523520/129‑1 (MQ=255)
                                                                            acgggTCCCCGCGCTTCCCGAAACG‑GTGG‑TTTATCTC‑AACCCTTAAA‑CCCCTT‑ATCAGGATAATATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                         <  2:1587456/147‑1 (MQ=255)
                                                                            aCAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                        >  2:1306116/1‑151 (MQ=255)
                                                                               ggTCGCCGCGCGTCCCGGACCC‑GGGG‑TTTATCCCGACCCTTAAAC‑CCGT‑T‑ATAAAGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                          <  2:722462/146‑1 (MQ=255)
                                                                               ggTCCCCGCGCGTACCGAACCC‑GTGT‑TTTATCTCAACCCGTAAAC‑CCGC‑T‑ACCGGGATCATATAATCAACACCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                      <  1:1453493/150‑1 (MQ=255)
                                                                                  cgccgcGCGT‑CCGGAACCGGTGG‑TGTATCCCGACCCGTAAAC‑CCGC‑T‑ATAAGGAAAATATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                       <  1:198928/146‑1 (MQ=255)
                                                                                   gccgcGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                 >  2:1421128/1‑151 (MQ=255)
                                                                                    ccgcGCGTCCCGAAACC‑GTGG‑TTTATCTCGACCCTTAAAA‑CCGC‑T‑ATAAGAATCATATGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                    <  1:552165/147‑1 (MQ=255)
                                                                                    ccgcCCGTCCGGGAACC‑GGGT‑TTTATCCCGACCCGTTAAC‑CCGC‑T‑ATAAGAATAATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                <  1:602975/151‑1 (MQ=255)
                                                                                     cgcgcgTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                >  2:1140358/1‑150 (MQ=255)
                                                                                       cgcgTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                >  2:288480/1‑148 (MQ=255)
                                                                                        gcgTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  1:977530/1‑150 (MQ=255)
                                                                                        gcgTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  1:731011/1‑150 (MQ=255)
                                                                                         cgTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  2:837544/1‑149 (MQ=255)
                                                                                          gtcccGAAACC‑GTGG‑TTTATCTCAACCCGTTAAC‑CCGC‑T‑ATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                           <  1:846654/147‑1 (MQ=255)
                                                                                          gTACCGGACCC‑GGGG‑TGTATCCCGACCCGTTAAC‑CCGC‑T‑ATAAGGATCATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                          <  1:748793/151‑1 (MQ=255)
                                                                                             ccGGACCG‑GGGG‑TGTATCTC‑AACCCTTAAACCCGC‑T‑ATCAGGATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                              <  1:1573012/144‑1 (MQ=255)
                                                                                              cGGACCC‑GGGG‑TTTATCTCGACCCTTTAAC‑CCGT‑T‑ATAAGAATAATATAACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                          <  1:923958/147‑1 (MQ=255)
                                                                                               ggACCC‑GGGG‑TTTATCTCGACCCG‑TAAA‑CCCC‑TAATAAGGAAAATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                     <  1:1601982/151‑1 (MQ=255)
                                                                                                ggaccggtgg‑ttTATCCCCAGCAGTAAAC‑CCGC‑T‑ATCAGGATTTTAGGACCACCCCCCCCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                           <  2:466239/139‑1 (MQ=255)
                                                                                                 aaCC‑GTGT‑TGTCTCCCGACCCTTAAAC‑CCTC‑T‑ACCAGGTTAATATAATAACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  1:1036823/151‑1 (MQ=255)
                                                                                                  aCC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                  >  1:1123285/1‑151 (MQ=255)
                                                                                                    c‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGACGAAGATAACCTGgcg                                                >  2:1621338/1‑151 (MQ=255)
                                                                                                    c‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                >  2:527988/1‑151 (MQ=255)
                                                                                                        gg‑tgTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               >  1:1209873/1‑149 (MQ=255)
                                                                                                           tgtATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           >  1:778280/1‑151 (MQ=255)
                                                                                                           tgtATCCCAACCCGTTAAC‑CCGC‑T‑ACCGGAATAATATAACAAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 <  1:1244447/145‑1 (MQ=255)
                                                                                                            ttttctcAACCCGTAAAC‑CCTC‑T‑ATCAGTATCATATGATAACCCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           <  2:483721/147‑1 (MQ=255)
                                                                                                                ctcGAACCTTAAAC‑CCCC‑T‑ATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                         <  2:1597894/148‑1 (MQ=255)
                                                                                                                ctcGAACCGTTAAC‑CCGC‑T‑ATAAGAATCATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                          <  1:286608/147‑1 (MQ=255)
                                                                                                                    aaCCTTTAAC‑CCGC‑T‑ATAAGAATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                             <  1:377814/140‑1 (MQ=255)
                                                                                                                    aaCCTTAAAC‑CCCC‑T‑AACAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGa                                  <  2:51677/151‑1 (MQ=255)
                                                                                                                     aCCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                 >  1:1477468/1‑151 (MQ=255)
                                                                                                                      cctttAAC‑CCGC‑T‑TGCAGAATCATATGACCAACCCCCTCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                    <  2:1625932/144‑1 (MQ=255)
                                                                                                                         ttAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             >  1:48707/1‑151 (MQ=255)
                                                                                                                            aacccGC‑T‑ATCAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           <  1:766207/148‑1 (MQ=255)
                                                                                                                                        cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                 >  2:976450/1‑151 (MQ=255)
                                                                                                                                        cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                 >  2:794807/1‑151 (MQ=255)
                                                                                                                                        cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                 >  2:439023/1‑151 (MQ=255)
                                                                                                                                         aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                >  2:580136/1‑151 (MQ=255)
                                                                                                                                         aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                >  2:616889/1‑151 (MQ=255)
                                                                                                                                          ggttaatatGATCACCCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                      <  2:1190636/139‑1 (MQ=255)
                                                                                                                                          ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt               >  2:322421/1‑151 (MQ=255)
                                                                                                                                             tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt            >  2:284588/1‑151 (MQ=255)
                                                                                                                                             tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt            >  2:1073738/1‑151 (MQ=255)
                                                                                                                                              cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg           >  1:532461/1‑151 (MQ=255)
                                                                                                                                              cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg           >  1:1666397/1‑151 (MQ=255)
                                                                                                                                               atatGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc              <  1:1184725/147‑1 (MQ=255)
                                                                                                                                               atatGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc              <  1:1087429/147‑1 (MQ=255)
                                                                                                                                                 ataaccAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg         <  1:1616494/145‑1 (MQ=255)
                                                                                                                                                 atGAACAACCCCCCCCCCATTCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa             <  2:1263800/146‑1 (MQ=255)
                                                                                                                                                  tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc        >  2:1485847/1‑150 (MQ=255)
                                                                                                                                                   aaTCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat      <  1:496473/150‑1 (MQ=255)
                                                                                                                                                       aaCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc              <  2:24848/139‑1 (MQ=255)
                                                                                                                                                       aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  1:1636754/1‑151 (MQ=255)
                                                                                                                                                        acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  1:1029810/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                   
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376356‑1376652

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