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

CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCG  >  NC_000913/1376361‑1376654
                                                                                                                                                           |                                                                                                                                                     
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTc                                                                                                                                                  <  1:1027982/150‑1 (MQ=255)
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCg                                                                                                                                                 <  2:1369189/151‑1 (MQ=255)
 cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGc                                                                                                                                                <  1:535295/151‑1 (MQ=255)
     gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGAt                                                                                                                                            >  1:1119449/1‑151 (MQ=255)
             cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCccgccg                                                                                                                                    >  2:65526/1‑151 (MQ=255)
               ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGtt                                                                                                                                  >  1:131879/1‑151 (MQ=255)
                       ccTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                          >  1:374747/1‑151 (MQ=255)
                         tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                         >  2:680953/1‑150 (MQ=255)
                                cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                 >  1:1220077/1‑151 (MQ=255)
                                   ccacaaaggggTGAGTCATTCCCGCCGCGTGCGGC‑AACGGGC‑CGCCGCGCG‑CCCCGAACC‑GGGGTTTTATCTCGAACCATAAAC‑CCG‑TT‑A‑TCAGCA‑TTATA‑TAACAAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                  <  1:40656/138‑1 (MQ=255)
                                       aaaaTGGGGAGTGCTCCCCTGCGTTTTCGGA‑AACGGGC‑CCCCG‑GGGGCCCCGAACC‑CTGGTTTTATCTCGAACGGTTAAC‑CCG‑CT‑A‑CCGGCA‑TCATA‑TGAAAAA‑CCCCCTCGCCTCTGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                           <  1:1358247/150‑1 (MQ=255)
                                        aaaTGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                          >  1:5715/1‑150 (MQ=255)
                                           tGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                      >  1:531589/1‑151 (MQ=255)
                                                agtgtTTCCGTGCGTTTTCGGA‑AACGGGC‑CGCCGCGGTTCCCGGACCC‑GTGG‑TGTATCAGGAACCGTTAAC‑CCG‑CT‑A‑ACAGGA‑TCATA‑TGACCAC‑CCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                                  <  1:1272019/150‑1 (MQ=255)
                                                agtgaTTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                    >  2:424804/1‑148 (MQ=255)
                                                 gttatccccTGCGCTTCCGGA‑AACGGGC‑CGCCGCGCGTCCCGGACCC‑GTGG‑TGTATCCCAACCCTTTAAC‑CCG‑CT‑A‑TAGGTA‑TCATA‑TGACAAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                                  <  2:878039/143‑1 (MQ=255)
                                                 gtgttccccTGCGC‑TTCTGGAAACGGGC‑CCCCGCGCGTCCCGAACCG‑GGGG‑TTTATCTCAAACCTTAAAC‑CCG‑CT‑A‑TAGGGA‑TAATA‑TGACCAC‑CCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                        <  2:1341872/137‑1 (MQ=255)
                                                       cccggggcgTAC‑GG‑AACCGGTCCGCCGCGCGCCCGGGACCG‑GTGT‑TGTACCCCAACCCGTAAAC‑CCG‑CT‑A‑TCAGAA‑TAATA‑TGACAAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                          <  2:1344882/145‑1 (MQ=255)
                                                       cccTGCGCGTCCTGA‑AAAGGGT‑CCCCGGGCGCA‑CGGAACCGGGGG‑TGTACCCCGACCCTTAAAC‑CCG‑CT‑A‑TAAGGA‑TAAAA‑TAACAAC‑CCCCCCCGCGGAGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                          <  2:371114/151‑1 (MQ=255)
                                                       cccTGCGCGTACTGA‑AACGGTT‑CGCCGCGCGTCCGGGAACG‑GTGG‑TTTATCTCGACCCGTTAAC‑CCG‑CT‑A‑TGAGGA‑TAATA‑TGATCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                             <  1:1244616/148‑1 (MQ=255)
                                                       cccGGCGCGTACGGA‑AACAGGC‑CGCCGCGCGTCCCGGACCG‑GTGT‑TGTCTCCCCACCCGTAAAA‑CCG‑CT‑A‑TCAGAA‑TTATA‑TGATCAA‑CCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                             <  2:1261498/148‑1 (MQ=255)
                                                         cTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                        >  1:1216828/1‑151 (MQ=255)
                                                          tGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                        >  2:362748/1‑150 (MQ=255)
                                                           gcgcGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                      >  2:729844/1‑151 (MQ=255)
                                                           gcgcGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                      >  1:574484/1‑151 (MQ=255)
                                                               gTACTGA‑ACCGGGC‑CCCCCCGCGT‑CCGGAACCGGTTG‑TTTATC‑CCAACCCTTAAA‑CCC‑CTAA‑TAAGTA‑TAATA‑TAAACAC‑CCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                   <  2:1269527/150‑1 (MQ=255)
                                                                 aCTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                                >  2:46712/1‑151 (MQ=255)
                                                                     gaaaCCGGT‑CGCCGCGTGTGCGGGACGG‑GGGT‑TGTTTCCCG‑ACCGTTAAACCCG‑TT‑A‑T‑AGGATTCATA‑TGACAAA‑CCCCCGCCCGGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                                <  2:1309513/146‑1 (MQ=255)
                                                                        acggtt‑CGCCGCGGGTCCGGGAACC‑GTGC‑TGTATCTC‑AACCCTTAAACCCG‑CT‑A‑TCAGTA‑TAATA‑TGACCAC‑CCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                            <  1:755374/144‑1 (MQ=255)
                                                                        aacagtcccccgcGTGTCCCGGACCC‑GGGG‑GGTCTCTCGACCCTTTAAA‑CCG‑TT‑A‑CAAGAT‑TAATA‑TGACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                   <  1:1122635/134‑1 (MQ=255)
                                                                         acggtcccccccgcgTCCGGGACCT‑GTGG‑TGTTCCCCGAACCTTAAAA‑CCC‑CT‑A‑TAGGTT‑TCAAA‑TGACAAC‑CCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                           <  1:1105179/139‑1 (MQ=255)
                                                                          cgggcccccgcGGGTACCGGACCC‑GGGG‑TGTATCCCGCCCCTTTAAC‑CAG‑CT‑A‑T‑AAGA‑TAAAAGTAACCAC‑CACCCGGGCGAAACCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                         <  1:49178/144‑1 (MQ=255)
                                                                            ggcccccgccgcccccGGACCC‑GTGG‑TTTCTCTCGACCCTTAAAC‑CCG‑CT‑A‑TCAGAA‑TAATA‑TAAACAC‑CCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                          <  2:245918/134‑1 (MQ=255)
                                                                               cgcccggcgccccGGACCG‑GTGG‑TGTATCTCG‑ACCTTTAAA‑CCGTCT‑A‑TAAGAA‑TCATA‑TAACAAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                       <  2:264370/137‑1 (MQ=255)
                                                                                 ccgcGCGTCCCGGACCC‑GTGG‑TGTTTCCCGACCCTTAAAC‑CCG‑CT‑A‑TCAGTA‑TCATA‑TACACAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                      <  1:927306/147‑1 (MQ=255)
                                                                                  cgcgcgTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                  >  2:360374/1‑150 (MQ=255)
                                                                                   gcgcgTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                  >  2:1004444/1‑149 (MQ=255)
                                                                                    cgcgTCCGGGACCC‑GGGG‑TGTACCTCAACCCGTAAAC‑CCC‑CT‑A‑AAAGGA‑TCATA‑TAACCAA‑CCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                                 <  2:278397/149‑1 (MQ=255)
                                                                                      ggtcccGGGACCGGTGG‑TGTACCCCAAACCGTAAAC‑CCT‑CT‑A‑TCAGGA‑TAATA‑TGACCAC‑CCCCCCCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                                <  2:566977/145‑1 (MQ=255)
                                                                                      cgTCCCGGAACC‑GTGG‑TGTCTCCCAACCCGTAAAC‑CCG‑CT‑A‑CAAGAA‑TAATA‑TAACCAC‑CCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                                 <  2:1340057/147‑1 (MQ=255)
                                                                                           ccgAACCGGTGG‑TTTATCCCGAA‑CTTTAAA‑CCC‑GT‑ATTAAGAA‑TAATA‑TAACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                         <  2:47087/149‑1 (MQ=255)
                                                                                           cGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                        >  1:474435/1‑151 (MQ=255)
                                                                                           cGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                        >  1:436545/1‑151 (MQ=255)
                                                                                                 ctgtgt‑tgTATCTC‑AACCTTTAAACCCG‑TT‑A‑TCAGAA‑TAATA‑TAACCAA‑CCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                    <  2:1243189/148‑1 (MQ=255)
                                                                                                 c‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                  >  1:1278325/1‑151 (MQ=255)
                                                                                                    tgg‑tgTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                 >  2:1063546/1‑150 (MQ=255)
                                                                                                    tgg‑tgTACCTC‑AACCCTTAAACCCG‑TT‑A‑TAAGGA‑TAATA‑TAACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                      <  2:1022003/145‑1 (MQ=255)
                                                                                                            tctcGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                         >  1:759637/1‑151 (MQ=255)
                                                                                                            tctcGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                         >  1:1274907/1‑151 (MQ=255)
                                                                                                               ccaaCCCTTAAACCCG‑TT‑A‑TCAGAA‑TCATA‑‑AACCAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                             <  2:787845/143‑1 (MQ=255)
                                                                                                                 acccGTTAAC‑CCG‑CT‑A‑TCGGGA‑TCATA‑TAACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                     <  1:435022/148‑1 (MQ=255)
                                                                                                                  aCCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                      >  1:567815/1‑148 (MQ=255)
                                                                                                                  aCCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                   >  2:1241854/1‑151 (MQ=255)
                                                                                                                   cctttAAC‑CCT‑CT‑A‑TCAGAA‑TAATA‑TAATAAC‑CCCCCGCCCGGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                     <  2:951239/145‑1 (MQ=255)
                                                                                                                    cGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                  >  2:1241197/1‑150 (MQ=255)
                                                                                                                         aa‑CCGCCT‑A‑TAAGGA‑TCATA‑TAATAAC‑CACCCGCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                             <  1:382750/135‑1 (MQ=255)
                                                                                                                           acccccT‑A‑TCAGGA‑TCATG‑TGACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                               <  1:1186858/143‑1 (MQ=255)
                                                                                                                                   a‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGTGCAGGCGCTGCATGACAACCGCGACCTGc                      >  1:760316/1‑150 (MQ=255)
                                                                                                                                      cAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                   >  2:513261/1‑151 (MQ=255)
                                                                                                                                        ggA‑TCATA‑TGAACAA‑CCCCCCCGCCGCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                        <  1:856671/144‑1 (MQ=255)
                                                                                                                                          a‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                >  1:230109/1‑150 (MQ=255)
                                                                                                                                          a‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                >  1:167037/1‑150 (MQ=255)
                                                                                                                                            tCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt              >  2:382501/1‑151 (MQ=255)
                                                                                                                                            tCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt              >  1:1241098/1‑151 (MQ=255)
                                                                                                                                             cATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg             >  1:1340862/1‑151 (MQ=255)
                                                                                                                                               ta‑tGACCAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                   <  2:543839/143‑1 (MQ=255)
                                                                                                                                                  tGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc          >  2:1119291/1‑150 (MQ=255)
                                                                                                                                                    aTCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat        >  2:107953/1‑150 (MQ=255)
                                                                                                                                                      cAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt     >  1:1206697/1‑151 (MQ=255)
                                                                                                                                                        a‑cacCCTCGCCGATGCCCGCCGTTACCTCCTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg    >  2:519903/1‑150 (MQ=255)
                                                                                                                                                          cacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCg  >  1:1046696/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                     
CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGG‑AACAGGT‑CGCCGCGCGTACCGGAACC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CT‑A‑TCAGGA‑TCATA‑TGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCG  >  NC_000913/1376361‑1376654

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