Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA49.2% 27.4 / 42.5 63L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base A (0/31);  total (18/45)
Fisher's exact test for biased strand distribution p-value = 2.62e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357740‑358038
                                                                                                                                                          |                                                                                                                                                       
cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                         >  2:894672/1‑151 (MQ=255)
   aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                      <  2:1080937/151‑1 (MQ=255)
      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTTCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                   <  2:1093957/151‑1 (MQ=255)
      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                   <  2:725644/151‑1 (MQ=255)
        cGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                                 <  2:497664/151‑1 (MQ=255)
            gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                                                                             >  2:831014/1‑151 (MQ=255)
               cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                          >  2:1059713/1‑151 (MQ=255)
                aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                                                         >  2:1248067/1‑151 (MQ=255)
                 tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                                        >  1:726310/1‑151 (MQ=255)
                  cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGCGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCCATTTTTCCGCCGCAACCAAACCGCCGTTCTCTcc                                                                                                                                       >  2:747198/1‑151 (MQ=255)
                   aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                       >  1:579912/1‑150 (MQ=255)
                       aaTCCCAGCCGCAGCCTTGTGTCAGCGCCTCCCTTGTTTCATTCCCGCCAACACCCCCGCTTTTATCAAACGCCATGGGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagta                                                                                                                                   <  2:694596/150‑1 (MQ=255)
                           ccAGCCGCAGCCCTGTGTCGGCGGCTCCCGTTTTTCTTTCCCGCAAACACCCGCGCTTTTACCCACCGCCATGTAAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                                <  1:1303261/149‑1 (MQ=255)
                           ccACCCCCAGCCTTGTGTCAGCGCCTCCCTTGATCCATTCCCCCCAACCCCCGCCCTTTTACCCACCCCCATGTAAAAAAATGCTTCCGCGGCGCTTTTCTGCAAGTTTTTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCagta                                                                                                                                   <  2:534347/146‑1 (MQ=255)
                            cccccGCAGCCTTGTTTCCGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGTTTTTACCCACCCCCATGCACAAAAATCCTGCCGCGGCGGTTGTCTGCAAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                                   <  2:479296/142‑1 (MQ=255)
                            cAGCTGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAGGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                             >  2:782552/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                             >  1:892113/1‑151 (MQ=255)
                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                            >  1:283682/1‑151 (MQ=255)
                               ccGCAGCCTTGTGTCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                          <  1:433027/151‑1 (MQ=255)
                                cGCAGCCTTGTGTCGGCGGCCCCCGTAATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTCAAAAAGTCCCTTCGCGCGGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                                <  2:772590/144‑1 (MQ=255)
                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                         >  2:691697/1‑150 (MQ=255)
                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                           >  1:415606/1‑149 (MQ=255)
                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                            >  2:1226613/1‑148 (MQ=255)
                                  cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                            >  2:1063172/1‑146 (MQ=255)
                                    ggccTGTGTCGGCGCTTGCCGTGATTCTTTCCGGCAAACACCCGCGCATTTCTCCAACGCAATGTCAAAAAAGCCTTCCGGGGCGGTTGTTGGCCGGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                     <  1:966288/149‑1 (MQ=255)
                                      cttgtgTCAGCGGCTCCCGTGATTCTTTCCCCCAAACACCCGCGCTTTTTCCCAACTCCATTCAAAAAATTGCTTCCGCGGCGCCTCCCTCCCCGTTTTTTCCCCCCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                   <  1:895282/148‑1 (MQ=255)
                                        tgtgTCAGCGCCTCCCGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                   <  1:916213/148‑1 (MQ=255)
                                                cgcctcccgtttttCATGCCCCCCAACACCCGCCCTTTTCCCCACCGCCATGAGAAAAAATGCCTCCGCGGGGGCTGTCTGCCGGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                                                                             <  2:356490/133‑1 (MQ=255)
                                                 ggCTCCCGTGTTTCTTTCCCCCCAAAACCCGCGCTTTTATCCAACGCCTTGTGAAAAAATGCTTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                                        <  1:1197288/149‑1 (MQ=255)
                                                            ttCATTCCCCCCAACACCCCCGCATTTCCCCACCGCCATGTGCAAAAATGCTTTCGCGCCGCCTTTCTCCCAGCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                                                                              <  1:417763/148‑1 (MQ=255)
                                                              cttgcccGCCAACACCCGCGCATTCCCCCACCGCCGGGGGCAAAAAGGCCGCGGGGGGGGCGGGCGGCGGGTTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                                 <  2:78909/138‑1 (MQ=255)
                                                              cttccccGCCAACACCCGCGCTTTAATCCACCGCCATGTGAAAAATGCCCTCGGCGGCGCTTGTCTCCCAGCTATTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                                                                              <  2:1103724/141‑1 (MQ=255)
                                                                ttccccccaccaCCCGCGCTTTTACCAACCCCATTGTAAAAAAATGCCTTCGCGCGGCCTGTCTGCCGTTTATTTGCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                           <  2:824138/141‑1 (MQ=255)
                                                                ttccccccAACACCCCCGCTTTTACCCACCCCATTGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                                                        <  1:109085/144‑1 (MQ=255)
                                                                 ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAACTTTTAGCCCCTCCtgtt                                                                                          >  1:37287/1‑147 (MQ=25)
                                                                  tCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGAGAGGAAAAAACTTTaaaaccctaccttttgt                                                                                       >  2:601505/1‑135 (MQ=21)
                                                                  ccccccAACACCCGCTCTTTCACCAAACGCCATGTCCAAAAATCCCTCCGGGGCGCTTTCCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                                         <  2:39567/143‑1 (MQ=255)
                                                                      ccaccaccccctcATTTTCCCACCCCCGTGTCCAAAAAGCCCTGCGCGGCGGCTGTCGGCCAGCTATTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                                 <  2:164420/139‑1 (MQ=255)
                                                                      ccAACACCCGCGTATTTACCCACCGCCATGTGAAAAATTGCCTTCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                                   <  1:862678/149‑1 (MQ=255)
                                                                        aacaccCGCGCTTTTATCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                               <  1:54560/151‑1 (MQ=38)
                                                                         acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAACTTTAACCCCTCATGTTTTGTGGGaa                                                                                 >  1:1272931/1‑148 (MQ=25)
                                                                          cacccccgattttccccAACGCCAGGTGAAAAAATGCTTTCGGGGGGGTGGTCTGCCAGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                              <  2:891505/135‑1 (MQ=255)
                                                                                   aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACACAACCGCCGTTCTCTCCAGTAGTGGCGGCGCACGCGAAATATCATTTAAACCGTCTGGTGGTGTGGGCGTCACTGCTgg                                                                       >  2:770576/1‑148 (MQ=21)
                                                                                    ttcccccAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCGGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                                    <  2:2457/145‑1 (MQ=255)
                                                                                            acccccttgacaaaaaaTGCCTCGGCGGCGCCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                           <  2:443190/139‑1 (MQ=255)
                                                                                               ccATGTAAAAAAACCCCTTCGCGGCGGTTGTCTGCCATCTTTTTCCCCGCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                                <  2:1022275/142‑1 (MQ=255)
                                                                                                ccttgtgAAAAAATGCCTTGGCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                                        <  1:1144117/147‑1 (MQ=255)
                                                                                                     acaaaaagTGCTTTCGCGGTGTCTGCTTCCCAGTTTTTTTCCCGCCAAAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                                        <  2:550665/142‑1 (MQ=255)
                                                                                                        aaaaaTGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                                <  1:1215201/150‑1 (MQ=255)
                                                                                                               ccTGCGGGGCGGCTGTCTCCCAGTTTTTTCCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                        <  1:225901/151‑1 (MQ=255)
                                                                                                                  tccgcgGGGCCTGCCGGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                                         <  2:176129/146‑1 (MQ=255)
                                                                                                                        ggcggtTTTCTCCCGGTTTTTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                        <  2:156785/144‑1 (MQ=255)
                                                                                                                        ggcggtTGTCTGCCAGCTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                                         <  1:73981/143‑1 (MQ=255)
                                                                                                                                    ccactttttttcCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                     <  1:266082/144‑1 (MQ=255)
                                                                                                                                    ccATCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                           <  1:1165358/145‑1 (MQ=255)
                                                                                                                                     ccgtcttttttCCCCCGAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                            <  2:532071/138‑1 (MQ=255)
                                                                                                                                       gttATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa                    <  1:579935/147‑1 (MQ=255)
                                                                                                                                       gCTTTTTCTCCCCGGCAACAAAAGCGGCCTTCTCTCCAGGAGTGGCGGGGCCAGAGGAATGGGTTGAAGGCCGTCATGTTGGGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                         <  2:889798/144‑1 (MQ=255)
                                                                                                                                           ttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                 <  1:987768/148‑1 (MQ=255)
                                                                                                                                           ttttcCCGCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                 <  2:360540/148‑1 (MQ=255)
                                                                                                                                           ttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa                <  1:1171027/149‑1 (MQ=255)
                                                                                                                                            ttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTTAACCCCCCCTGTTGTGTGGGGATCGCTGCTGGGGACAAATGGGGAAGGGGAGGGGGGGGGAACAACTCCCGAAACCCGCCAAtt                  >  1:300546/1‑146 (MQ=255)
                                                                                                                                            tttcCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt              <  1:701640/150‑1 (MQ=255)
                                                                                                                                               cccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg           <  2:1310865/149‑1 (MQ=255)
                                                                                                                                                        aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa  <  2:1238899/151‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                       
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357740‑358038

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