Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA53.6% ‑0.5 / 58.1 56L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (17/9);  new base A (0/30);  total (17/39)
Fisher's exact test for biased strand distribution p-value = 3.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-03

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357742‑358027
                                                                                                                                                           |                                                                                                                                            
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACCACCGCGCATTTATCCAACGCCATGTGCCAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                            >  1:571782/1‑151 (MQ=255)
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                            >  2:1644016/1‑151 (MQ=255)
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                            >  1:439492/1‑151 (MQ=255)
     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                       >  2:517849/1‑151 (MQ=255)
     gCGTAGCGCACCACGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAACCCCCCATGTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCGGTCTGGCTGCTATTGTTACTCCGCGACAAAAccg                                                                                                                                       >  1:1690266/1‑151 (MQ=17)
        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                    >  2:1279743/1‑151 (MQ=255)
                 aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa                                                                                                                           >  2:1161824/1‑151 (MQ=255)
                  ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                                                          >  1:1341210/1‑151 (MQ=255)
                   gaatttCCCGCCGCACCCCTGGGGCGGCGGCTCCCGTGTTTATTTTCCCCCAACACCCGCGTTTTTATCCAACGCCTTGTGAAAAAATGCCTCCGGGGCGGCTTTCTCCCAGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTcc                                                                                                                            <  2:1207122/144‑1 (MQ=255)
                   cattccccACCCGCAGACTTTTGTCAGCGGCTCCGGTATTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCATTGTCAAAAAATCCCTCCGCGGCGCCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagt                                                                                                                         <  1:1192672/145‑1 (MQ=255)
                      aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                      >  1:1649749/1‑151 (MQ=255)
                        cccccccGCGGCCCGGTGCCAGCGCCCCCCGTTTTTCTTTCCCCCCAACACCCGCGTTTTTACCCACCCCCATGTGAAAAAATGCCTCCGGGGCGGCTTTCCCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    <  2:1457282/146‑1 (MQ=255)
                          cacccGCAGCCCTGTGCCAGGGGCTCCGGTTATTCTTTCCCCCAAACACCCGCCCTTTTACCCACCCCCATGGGAAAAAATCCTTCCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                 <  2:1475505/148‑1 (MQ=255)
                            gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                              >  1:1013436/1‑151 (MQ=255)
                                  gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCCACAAAACCGCCGTTCTCTCCAGTAGGGGCggggggaga                                                                                                        >  2:1085092/1‑146 (MQ=255)
                                   ccgtttgtCGGCGCCTCCCTTGTTTCTTTCCCCCCAACAACCGCACTTTAACCCACCGCCATGTAAAAAAATGCTTTCGGGGGGCCTGTCTGCAATCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                      <  2:94825/146‑1 (MQ=255)
                                   aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCagag                                                                                                      >  2:1382828/1‑151 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgggaaaa                                                                                                   >  1:238201/1‑141 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaggaaa                                                                                                   >  1:898773/1‑141 (MQ=255)
                                        tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaagaaaa                                                                                                  >  2:1052313/1‑140 (MQ=255)
                                         tgtcggcgcctccctGATTCCATTCCCGCCACCACCCGCGATTTTACCCACCCCCTTGTGAAAAAAGCCCTTCGCGGCGGTTGTCTCCCGGCTTTTTGCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                  <  2:916136/136‑1 (MQ=255)
                                          gcAGCGGCTCCGGTGTTTCTTTCCCCCAACCACCCGCGCTTTTACCCACCGCCATGTAAAAAATTGCCTCCGCGGGGGCGGTCGCCCATTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                <  2:934236/149‑1 (MQ=255)
                                          ccggCGGCTACCGGTATTCTTTCCCCCCAACACCCGCGCTTTTACCAACCGCAAGGTAAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                               <  2:857561/148‑1 (MQ=255)
                                              cGGCTCCGGTAATTCTTTCCCCCAAACACCCGCGCATTTACCCAACCCCAGGTGAAAAAATGCCTCCGGGGGGGCGGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                             <  2:189769/149‑1 (MQ=255)
                                                 ctcccGTGTTTCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                        <  2:1003489/148‑1 (MQ=255)
                                                  tccggtgtTTCTTTCCCCCCAACACCCCCCATTTTACCCACCGCAAGTGGAAAAATGGCCTCGGGGGGGGCTGTCGCCCAGTTTTTTTCCCCCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                       <  2:595868/147‑1 (MQ=255)
                                                  caCCGTGATTCATTCCCCCCAACACCCGCCTTTTTCCCCACCCCAGTGTGAAAAAATCCCTTCGCGGCGGTTTTCTCCCGGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                       <  1:489991/150‑1 (MQ=255)
                                                    ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAGGGAAAAAATTTTAAAcc                                                                                       >  2:17377/1‑150 (MQ=25)
                                                    ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaaagtttaaacc                                                                                      >  2:879901/1‑137 (MQ=25)
                                                        taTTCTTTCCCGCCAACAACCGCCCTTTTACCCACCGCCATTTGCAAAAATCCCTCCGCGGCGCTTGTCGGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                                 <  1:241940/150‑1 (MQ=255)
                                                          ttCTTTCCCCCCAACACCCGCCCTTTTACCCACCGCAATGTGAAAAAAGCCCTTCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                               <  2:1482349/151‑1 (MQ=255)
                                                             ctttccccccAACACCCGCCATTTTACCCAACGCAATGTGAAAAAATCCCTCCGCGGCGGTTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1487730/143‑1 (MQ=255)
                                                             cATTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1197807/151‑1 (MQ=255)
                                                             cATTCCCCCCAACACCCCCCCTTTTACCCAACCCCATGTGAAAAATTGCCTTCGGGGGGGCTGTCGGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1425841/151‑1 (MQ=255)
                                                                  cccccaccaCCCGCGTATTAACCCAACGCCAGGGCAAAAATTCCCTCCGCGGGGGTTGTCGGCCAGTTTTTTCCCCGCAGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  1:24692/148‑1 (MQ=255)
                                                                  ccccaaacaccCGCGCTTTTCCCCACCCCCGTGTACAAAAATCCCTCCGCGGGGGCTGTCGGCCAGTTATTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  2:1401319/146‑1 (MQ=255)
                                                                            acCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGGCAGAGAGAAAAGTCTCAAGCCCCTCAGTTTTTGTGGCacccccg                                                                >  1:71896/1‑145 (MQ=21)
                                                                                  cATTTTTCCACCCCCTTGTCCAAAATTCCCTCCGCGGCGGCTGTCTGCCAGCTATTTTTCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                         <  2:1217167/151‑1 (MQ=255)
                                                                                  cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGGAAAAAACTTTAACCCCCCCCCGTTGTTGGGAAGCgctgcgcggg                                                         >  1:1380470/1‑146 (MQ=25)
                                                                                    tttACCCACCGCAATGTGAAAAAATGCTTCCGGGGGGCCTGTCGGCCGTCTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  1:505403/151‑1 (MQ=255)
                                                                                    tttACCCAACCCCATGTAAAAAAATCCCTTCGGGGGGGCTGTCGCCAGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  2:1205507/151‑1 (MQ=255)
                                                                                    ttaacccACCGCCGGGTGAAAAAATCCTTTCGCGGGGGCTGCCCCACAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  2:639640/146‑1 (MQ=255)
                                                                                        ctcaaccgctgggtgAAAAAATGCCTGCGGGGCGGCTTTCGGCCAGTTATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                    <  1:626670/139‑1 (MQ=255)
                                                                                        cccACCGCCATGTGAAAAAATCCTTTCGGGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                   <  2:1063106/150‑1 (MQ=255)
                                                                                            aCGCCAGGGGAAAAAAGCCCTTCGGGGGGGCTTTCGCCAAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                               <  2:1299161/151‑1 (MQ=255)
                                                                                             cccccttgtaaaaaaaTCCCTTCGGGGGGGCTGTCTGCCATTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                               <  2:1654513/140‑1 (MQ=255)
                                                                                             cGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGAAAAAATTTTTAACCCCCCATGTTTTTTGGGCAACGCTGCTGGtgaaaaagaggga                                              >  2:887553/1‑139 (MQ=25)
                                                                                             acccccatttaaaaaaaTCCTTCGGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                <  2:1564697/140‑1 (MQ=255)
                                                                                                   tgtgAAAAAATCCCTCCGGGGGGGCTGCCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  2:365101/151‑1 (MQ=255)
                                                                                                   tgaaaaaaaaTCCCTCCGCGGCGCCTTTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  1:1161254/146‑1 (MQ=255)
                                                                                                                              ttGTCTGCCGCTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:684132/150‑1 (MQ=255)
                                                                                                                              cttccTGCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:637447/147‑1 (MQ=255)
                                                                                                                              cggccggCCAGTTTTTTTCCCGCCCAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:517635/145‑1 (MQ=255)
                                                                                                                              cTGTCTGCCGGTTTTTTTCCCGCCACAAAAAGACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:1555412/151‑1 (MQ=255)
                                                                                                                              cTGTCGGCCAGCTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  2:1352039/151‑1 (MQ=255)
                                                                                                                              cTGTCAGCCAGTTGTTTCCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:246352/151‑1 (MQ=255)
                                                                                                                                           tttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg    <  2:393153/150‑1 (MQ=255)
                                                                                                                                             tttcCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:439645/151‑1 (MQ=255)
                                                                                                                                              ttccccgccgAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:1389164/147‑1 (MQ=255)
                                                                                                                                                           |                                                                                                                                            
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357742‑358027

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