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

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGC  >  NC_000913/1376356‑1376646
                                                                                                                                                         |                                                                                                                                             
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAa‑cac                                                                                                                                              >  2:410061/1‑151 (MQ=255)
  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAa‑cac                                                                                                                                              <  1:941966/149‑1 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTc                                                                                                                                          <  2:1040130/151‑1 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGcc                                                                                                                                       <  1:55301/151‑1 (MQ=255)
           aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATg                                                                                                                                   <  1:1162336/151‑1 (MQ=255)
                  cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCccgccg                                                                                                                            <  2:660740/151‑1 (MQ=255)
                  cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCccgccg                                                                                                                            >  2:1214173/1‑151 (MQ=255)
                   gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGt                                                                                                                           >  1:1029555/1‑151 (MQ=255)
                     cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGtt                                                                                                                          >  2:352782/1‑150 (MQ=255)
                               ggACGGCGCCCAAAAAATGGGGGGTAATCCCCTGCGCGTACTGGAACAGTCCGCCGCGCGTCCCGGACCG‑GTGGTGTATCTC‑AACCCTTAAACCCGT‑TATCAGTATAAAATGACCAC‑CCCCCCCCCCGATGCCCGCCGTTACAtcgtcg                                                                                                                 <  1:841855/149‑1 (MQ=255)
                               ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                >  2:1337518/1‑150 (MQ=255)
                                            aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                  >  2:301978/1‑151 (MQ=255)
                                            aaaaGGGGGAGTGATCCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGACCC‑GTGGTGTATCTCGAACCGTAAAA‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                   <  2:214254/150‑1 (MQ=255)
                                             aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                  >  1:1256511/1‑150 (MQ=255)
                                              aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                >  2:922008/1‑151 (MQ=255)
                                                 ggggagttaTCCCCGGCGGGTTCTGAAACGGGCCCCCGCGCGTCCCGGACCC‑GTGTTGTATCCCGACCCGTTAAC‑CCCC‑TATCAGGATCATATGACCAC‑CCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                              <  1:987282/147‑1 (MQ=255)
                                                    gagttaTTCCCGGCGCTTTCGGAAACAGGTCGCCGCGCGTACGGGACCG‑GGGGTGTATCTCGACCCGTTAAC‑CCGC‑TATGATTATTATATAATCAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                           <  1:1259536/150‑1 (MQ=255)
                                                              cTGCGGTTACTGGAACAGGCCCCCGCGCG‑ACCCGAACCGGGGGTTTACCTCGAACCGTTAAC‑CCCC‑TACCAGAATCAAATGACCAC‑CCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                 <  1:175820/150‑1 (MQ=255)
                                                              cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                >  2:1073180/1‑151 (MQ=255)
                                                                gcgcGTACGGGAACAGGGCGCCGCCCGTACGGGACCC‑GTGGTTTATCTCAACCCTTTAAC‑CCGT‑TATCAGGATAATATTTACAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                <  2:734457/149‑1 (MQ=255)
                                                                 cgcgTACGGGAACAGGCCCCCGCGCGTACCGGAACC‑GTGGTTTATCTCGAACCGTTAAA‑CCGT‑TATCAGGATCATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                               <  2:307217/149‑1 (MQ=255)
                                                                   cgTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                            >  1:18223/1‑150 (MQ=255)
                                                                     tatggcaaagggcccccgcGCGTGCGGCAGCG‑GTGGTTTATCCCAACCCTTTAAC‑CCGC‑TATAAGGTTCATAGGATCAC‑CACCCCCGCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                          <  1:1148730/135‑1 (MQ=255)
                                                                     tagggaaaCAGGTCCCCGCGCGTCCCGCACCG‑GGGTTTTATAACGAACCTTAAAC‑CCGC‑TATAAGAATCTAATGATCAC‑CACCCTCGCTGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                               <  2:831379/139‑1 (MQ=255)
                                                                      aCTGGAACAGGTCGCCGCGCGTACCGGACCC‑GTGGTGTATCCCGAACCGTTAAC‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                        <  2:881059/151‑1 (MQ=255)
                                                                           aaCAGGCCCCGGCGCGTCCCGGAACG‑GGGTGGTTCCCCGACCCGTAAAC‑CCGT‑TATCAGGATAATATGACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                    <  1:262475/150‑1 (MQ=255)
                                                                               ggTCGCCGCGCGTACCGGACCC‑GTGGTGTATCCCGACCCGTAAAC‑CCGT‑TAAAAGGATAATATGACAAC‑CCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                               <  1:770268/151‑1 (MQ=255)
                                                                                         cgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                       >  2:15524/1‑149 (MQ=255)
                                                                                            cccgaaccc‑GTGGTGTATCTCAAACCTTTAAC‑CCGT‑TATCAGGATCATATAATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                  <  2:583372/144‑1 (MQ=255)
                                                                                              cggaagccttggttTACCTCG‑ACCGTTAAA‑CCCCTTTTCAGATAAATATGATCAC‑CCCCCCCGCCTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                              <  2:148122/129‑1 (MQ=255)
                                                                                               ggAACC‑GGGGTTTATCCCGACCCTTTAAC‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                  <  2:963792/148‑1 (MQ=255)
                                                                                                 aaCG‑GTGGTGTATCTCGAACCGTTAAC‑CCGC‑TATCAGGATCATATGATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                             <  2:598007/151‑1 (MQ=255)
                                                                                                  aCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                            >  2:733757/1‑151 (MQ=255)
                                                                                                   cc‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                           >  2:566902/1‑151 (MQ=255)
                                                                                                      gtggtgTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                         >  1:88960/1‑151 (MQ=255)
                                                                                                      ggggtgTCTCCCAAGCCGTTAAC‑CCGC‑TATAAGAAAAATATAACCAC‑CCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                 <  1:135362/141‑1 (MQ=255)
                                                                                                       gtggtatacctccaCCGTTAAACCCAT‑TATCAGTATAATAGGATCAA‑CCCCCCCTCCTTGTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                              <  1:591820/133‑1 (MQ=255)
                                                                                                          tgtATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                     >  2:778968/1‑151 (MQ=255)
                                                                                                            tATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                   >  2:989522/1‑151 (MQ=255)
                                                                                                             aTCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                  >  1:1164317/1‑151 (MQ=255)
                                                                                                                tcGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                 >  2:162605/1‑149 (MQ=255)
                                                                                                                 ccgaacGTTAAACCCGA‑TATCGGGAT‑AAATGAACAACCCCCCCCGCCGTGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                               <  1:873306/146‑1 (MQ=255)
                                                                                                                  caaCCCTTAAA‑CCCCTTATAGGAATAATAGGAACAC‑CCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                  <  1:1055213/146‑1 (MQ=255)
                                                                                                                   aaCCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                             >  1:1158442/1‑150 (MQ=255)
                                                                                                                    aCCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                           >  1:643022/1‑151 (MQ=255)
                                                                                                                    aCCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                           >  1:146476/1‑151 (MQ=255)
                                                                                                                      cGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                          >  1:718930/1‑150 (MQ=255)
                                                                                                                      cGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                          >  1:150663/1‑150 (MQ=255)
                                                                                                                           aaccccc‑TACCAGAATCATATGAACAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                        <  1:73878/142‑1 (MQ=255)
                                                                                                                            aacccgttacaaGGATAATAGAACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                      <  1:298915/139‑1 (MQ=255)
                                                                                                                            a‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                    >  2:141759/1‑150 (MQ=255)
                                                                                                                             agccgttacaaGGAACATATGACAAA‑CCCCCTCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                     <  2:345025/140‑1 (MQ=255)
                                                                                                                              ccGC‑TATCAGTATAATATGACCAA‑CCCCCTCGCCTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                     <  1:4245/148‑1 (MQ=255)
                                                                                                                               cGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                 >  1:651618/1‑151 (MQ=255)
                                                                                                                                gC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                >  2:499497/1‑151 (MQ=255)
                                                                                                                                       aGGATCAAATGACAAC‑CCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg               <  2:902315/146‑1 (MQ=255)
                                                                                                                                        ggATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt         >  2:474474/1‑151 (MQ=255)
                                                                                                                                        ggATCATATGACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc            <  2:706704/148‑1 (MQ=255)
                                                                                                                                        ggATCATATAATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt             <  1:1311091/147‑1 (MQ=255)
                                                                                                                                          aTCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc        >  1:970803/1‑150 (MQ=255)
                                                                                                                                          aTCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc        >  2:937456/1‑150 (MQ=255)
                                                                                                                                            cATAGGATCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg     <  2:542471/151‑1 (MQ=255)
                                                                                                                                             atatGATCAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc            <  2:731656/143‑1 (MQ=255)
                                                                                                                                             atatGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt    >  2:682712/1‑151 (MQ=255)
                                                                                                                                             atatGACCAC‑CCCCCCCCCCGACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc            <  1:210828/143‑1 (MQ=255)
                                                                                                                                              tatGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg   >  1:307116/1‑151 (MQ=255)
                                                                                                                                              tatGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg   >  1:594021/1‑151 (MQ=255)
                                                                                                                                              tatGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg     >  2:941619/1‑149 (MQ=255)
                                                                                                                                               atGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc  >  1:212818/1‑151 (MQ=255)
                                                                                                                                                tGATCAA‑CCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg   <  2:836667/149‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                             
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGC  >  NC_000913/1376356‑1376646

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