Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA43.8% 125.9 / 56.7 105L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (42/17);  new base A (0/46);  total (42/63)
Fisher's exact test for biased strand distribution p-value = 1.03e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 109 total aligned reads displayed.

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357742‑358038
                                                                                                                                                      |                                                                                                                                                        
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                         >  2:723870/1‑150 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                        >  4:904540/1‑150 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       <  2:368107/151‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       <  2:1324119/151‑1 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                      <  2:707722/150‑1 (MQ=255)
     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                   >  4:761616/1‑151 (MQ=255)
     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                   >  2:991939/1‑151 (MQ=255)
     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                   >  4:262137/1‑151 (MQ=255)
      cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                                  >  2:360030/1‑151 (MQ=255)
         aGCGCATCAGGCAATTCCAGCCGCAGCCTTGTTTCAGCGGCACCAGTGATTCATTCCCGCCAACACCCGCGCATTAATCCAACTCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                                <  4:217965/150‑1 (MQ=18)
         aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                               >  2:1447584/1‑151 (MQ=255)
           cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                             >  2:1422943/1‑151 (MQ=255)
              aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                                                          >  2:449078/1‑151 (MQ=255)
                 aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa                                                                                                                                       >  1:457302/1‑151 (MQ=255)
                     atttccccccGCAGCCGTGTGTCGGCGCCTCCCGTGTTTCTTTCCCCCCACCACCCGCGCTTTTACCCAACCCCTTGTGCAAAAATCCCTCCGCGGCGGCTGTCTGCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                                   <  4:1367532/143‑1 (MQ=255)
                     atttCCAGCCGAAGACCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTATCCACCGCCATGTGCAAAAATGCTTTCGGGGCGGCTGTCTGCCACCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                                   <  1:525304/149‑1 (MQ=18)
                     aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                     >  4:730292/1‑149 (MQ=255)
                         ccAGCCTCAGCCCTGTGTCAGGGGCTCCCGTGATTAATTCCCGCCAACACCCGCGCATTTACCCAACCCCTTGTGAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                               <  4:29075/151‑1 (MQ=255)
                         ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                               >  4:972568/1‑151 (MQ=255)
                          cAGCCGCAGCCCTGTGCCAGCGGCTACCGTTATTCTTTCCCGCAAACACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTTCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                              <  2:789862/151‑1 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                              >  2:1423815/1‑151 (MQ=255)
                           aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                             >  4:848652/1‑151 (MQ=255)
                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg                                                                                                                         >  4:1106068/1‑149 (MQ=255)
                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg                                                                                                                         >  2:212207/1‑149 (MQ=255)
                               ccggACGTGTGAGAGGGGCGACTGTGATCCTTTCCCGCCAACAACCGCGCATTAATCCAGCGCGATGTAAAAAAATGCCTTCCGGGCGGCTGTCTACCGGGTATTTGTCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                         <  2:1393002/148‑1 (MQ=16)
                               ccagccttgtgTCGGCGGCTACCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTCCGCGGCGGCTGCCGGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                         <  4:1233297/144‑1 (MQ=255)
                                 agcccTGTGTCAGCGGCTCCCGTGTTTCATTCCCGCCAACACCCGCGCATTTATCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                                       <  2:1421307/148‑1 (MQ=255)
                                  gACTTGTGTCAGCGGCTCCCGTATTCCATTCCCCCAACCACCCGCGCTTTTACCCACCCCCATGTGAAAAAATCCCTCCGCGGCGCCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                      <  4:1040781/151‑1 (MQ=255)
                                  gACCTGTGTCAGGGGCTACCGTGTTTCATTTCCCCCAACAACCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                      <  2:883025/151‑1 (MQ=255)
                                   actggtgtCACCGGCTCCCGTGATTCATTCCCCCAAACACCCCCGCTTTTATCCACCGCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCGGCAAGCTTTTTTCCCGCGACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                    <  3:20288/147‑1 (MQ=18)
                                    ccTGTGTCGGCGGCTCCCGTGTTTCATTTCCGCCAACAACCGCGCATTTACCCACCGCCATGGGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                                  <  2:1423092/151‑1 (MQ=255)
                                     cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaga                                                                                                                    >  4:122787/1‑143 (MQ=255)
                                      tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcggg                                                                                                                     >  4:1051733/1‑142 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaggaa                                                                                                                  >  4:1216314/1‑141 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGccaaaaaaa                                                                                                                  >  2:1080993/1‑142 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGACgagataa                                                                                                               >  3:474650/1‑148 (MQ=255)
                                         gtCAGCGGCTCCCGTGATTCATTCCCGCCAACAACCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                             <  2:982231/151‑1 (MQ=25)
                                         gccggcgcCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCTTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                             <  4:331378/147‑1 (MQ=255)
                                          tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaa                                                                                                             >  4:801666/1‑147 (MQ=255)
                                              cGGCTACCGTGATTCATTTCCCCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                                        <  2:1233179/151‑1 (MQ=255)
                                               ggCTACCGTGATTCTTTCCCGCCAACACCCGCGCTTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                                       <  2:1185234/151‑1 (MQ=25)
                                                 cTACCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATCCTTTCGCGGCGGCTGTCTGCCACTTATTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                                     <  1:190382/151‑1 (MQ=25)
                                                  tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAACttttaacc                                                                                                    >  2:1284696/1‑146 (MQ=255)
                                                   cccgtttttCTTTCCCCCCACCACCCGCGCATTTCCCCACCGCCTTGTAAAAAAATGCTTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                                     <  4:1422481/142‑1 (MQ=255)
                                                   cccGTGTTTCATTTCCGCCAACAACCGCGCTTTTACCAAACCCCATGTGCAAAAAGGCTTCCGCGGCGGCTGTCTCCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                                     <  2:1104562/148‑1 (MQ=255)
                                                   aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGaaaaaaacttaaacc                                                                                                    >  1:414803/1‑138 (MQ=255)
                                                    ccgtatttCATTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCTTCCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                                  <  4:1364152/145‑1 (MQ=255)
                                                      gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACttttaaaccccc                                                                                                >  2:628715/1‑142 (MQ=25)
                                                        gATTCATTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                                              <  2:912299/151‑1 (MQ=37)
                                                          ttCATTCCCGCCACCACCCGCGCTTTTACCCAACCCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                                             <  4:439113/150‑1 (MQ=255)
                                                            cATTTCCGCCAACACCCGCCCATTTACCCAACGCCATGTGCAAAAATGCCTCCGCGGCGGCTTTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                                           <  2:1215652/150‑1 (MQ=38)
                                                            cATTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                                          <  1:358644/151‑1 (MQ=37)
                                                             ctttCCCCCAAACACCCGCTCATTTACCCACCGCCATGGCAAAAAATGCCTCCGCGGCGGTGGCCTGCAGGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                                          <  3:69406/150‑1 (MQ=255)
                                                             aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAGATTTAAccccccctgtttt                                                                                          >  4:722925/1‑142 (MQ=25)
                                                               ttCCCCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                                                                                        <  2:498897/150‑1 (MQ=255)
                                                                ccccccAACACCCCCGCTTTTTCCCACCCCCATGTGAAAAATTGCCTTCGCGCCGGCTGCCTCCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                      <  2:1417728/147‑1 (MQ=255)
                                                                 cccccAACACCCGCGCATTTACCCAACGCAATGTGAAAAAATGCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                                     <  1:551792/148‑1 (MQ=255)
                                                                 ccGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCGCCATGTTTTTTgg                                                                                     >  2:1289536/1‑151 (MQ=25)
                                                                    ccaaacaccCGCGCATTTACCCAACCCCATGTGAAAAAAGCCCTTCGCGGCGGCTGTCTGCCGGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                                   <  2:896112/148‑1 (MQ=255)
                                                                      caaaaacccGCGAATTTACCCACCGCCATGTGAAAAAAGGCCTTCGGGGCGGCTGTCGGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                                  <  1:224229/144‑1 (MQ=255)
                                                                       aacaccCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                                  <  2:669551/103‑1 (MQ=255)
                                                                       aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                                  >  4:669548/1‑103 (MQ=255)
                                                                          cacccGCGCTTTTACCCACCGCCATGTCCAAAAAGCCCTCCGGGGCGGCTGCCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                              <  4:1336974/148‑1 (MQ=255)
                                                                               gctcATTTCCCCACCGCCTTGTAAAAAAATGCCTCCGGGGCGGTTCGCCCCCCCTTTTTTTCCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGAAATAGctggtg                                                                        <  4:830259/148‑1 (MQ=255)
                                                                                 gcttttaccaaccGCCATGTAAAAAAATGCTTTCGCGGCGGCTGTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                      <  2:1207005/139‑1 (MQ=255)
                                                                                  tatttccccAGCGCTATGTGCAAAAATGCCTTCGGGGCGGTTGTTTACAAGATATTTTTCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                      <  2:1427528/143‑1 (MQ=255)
                                                                                  cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAAATTTATCCCGTCATATTTTGTGGGaaaggcttccggg                                                                       >  2:77141/1‑139 (MQ=25)
                                                                                   aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGAAAAAAGTTTTAACCCCCCCTGTTGTGTgggggtcgcggccgggg                                                                      >  4:1460577/1‑135 (MQ=25)
                                                                                    tttaccaacccccATGTGAAAAAATGCCTGCGCGGGGGCTTTCGGCCAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                                    <  4:546657/139‑1 (MQ=255)
                                                                                    tttACCCACCCCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                                   <  1:91023/151‑1 (MQ=255)
                                                                                      tacccacccccATGTGCAAAAAGGCCTCCGCGGCGGTTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                                 <  1:514636/142‑1 (MQ=255)
                                                                                       acccAACGCCATGTAAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                                     <  3:436155/113‑1 (MQ=255)
                                                                                        ccaaccGCAATGTGAAAAAAGGCTTCCGGGGCGGCTTTCTGCCGGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                               <  3:419756/146‑1 (MQ=255)
                                                                                              cccATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                                        <  3:87145/150‑1 (MQ=255)
                                                                                               ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAAAAAAAGATTTTAGCCCGGCACGTTGTGTGGGAAACGCTGCTTGGTACCaaaaggaa                                                           >  2:232322/1‑143 (MQ=37)
                                                                                                    tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCCCCCCGTTGTGTGGGGATCGCTGCTGGTGACAATGTGGGAAgggg                                                    >  4:1458743/1‑148 (MQ=255)
                                                                                                     gCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGAAAAAAATTTTAAACCCCCCTTTTGTGTGGGGATCCCTGCCGGGCACAAAATGGGAAAGGAAgg                                                 >  4:1066061/1‑150 (MQ=255)
                                                                                                      cAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAGAAAAAAATTTTAAACCCCCATGTTGTGTGGGAATCGCTGCTGGGAACAATGTGggaagggaag                                                  >  2:412348/1‑149 (MQ=255)
                                                                                                        aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGGGAAAAAGCATTAAGCCCGCCCTGTTTTTTGGCGATCGCCGCGGGTGACAATGTGGGAAGGGagggggg                                              >  4:65584/1‑145 (MQ=255)
                                                                                                        aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAACTTTAACCCCCCATGGTTTGTGGGGAACGCCGCTTGGAAACATTTGGGAAGGGAAGGgcg                                              >  4:1139846/1‑151 (MQ=255)
                                                                                                             gctttCGGGGCGGCTGTCTCCCAGCTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa                                          <  2:1226729/147‑1 (MQ=255)
                                                                                                              ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGAAGAAAAAAACTTTAAAACGCCATGTTGTGTGGGGATAGCCGCTGGGAAAAAAGTGGGAAGGGGAggggggggga                                         >  2:132701/1‑141 (MQ=255)
                                                                                                                 ccgcgGCGGTTGTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                                     <  4:1423818/150‑1 (MQ=255)
                                                                                                                    ggggggcTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                                  <  4:1271465/147‑1 (MQ=255)
                                                                                                                      ggggcTTTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                <  4:1243140/149‑1 (MQ=255)
                                                                                                                          cTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAACCCCCCCTTTTTTGTGGGAATAGCTTCTGGTGACAATGTGGGAAAGGGAGGGGGGGGGAACAGCTCCaaaa                             >  2:780670/1‑150 (MQ=255)
                                                                                                                            gccTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                            <  2:1297782/147‑1 (MQ=255)
                                                                                                                            gTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGAAAAAACTTTAAACCCCGCAAGTTGTGTGGGGATTGCCGCTGGTGAAAAAGTGGGAAGGGGAGGGCGGCGGAACAGCTCCAGAAccc                          >  2:574983/1‑150 (MQ=255)
                                                                                                                                gCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg                        <  4:183887/149‑1 (MQ=255)
                                                                                                                                   agtttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt                   <  2:801669/146‑1 (MQ=255)
                                                                                                                                   aGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAGGAAAAACTTTTAAACCCCCATGTTGTGGGGGGAATGCCGGTGGTAACAATGTGGGAAGGGAAGGGCGGGGGAACAACTCCCGAAACgcgcccat                   >  2:467379/1‑147 (MQ=255)
                                                                                                                                      ttttttcCCGCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                 <  4:1248303/148‑1 (MQ=255)
                                                                                                                                       atttcgccgccgCAACAAAGCCGCCGTTCTCTCCAGTAGTTGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg               <  4:1313127/145‑1 (MQ=255)
                                                                                                                                       aTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGCAGAGAAAAAAACTTAAAACCCTCCCGTTGTGTGGGGGTCCCTGCTGGGGACAAAATGGGAAGGGGAAGGGGGCGGAGCACCTCCCGAACCCCGCCaaatta                >  4:171259/1‑146 (MQ=255)
                                                                                                                                         ttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAATTAAACCCCTCCTGTTTTGTGGGAATTGCTGGTGGTGAAAAAGTGGGAAGGGGAGGGGGGGGGAACAACTCCCCAACCGCGCCAAATGAggt             >  2:1463101/1‑149 (MQ=255)
                                                                                                                                         tttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt             <  1:178900/151‑1 (MQ=255)
                                                                                                                                           tcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTAGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc            <  4:306493/148‑1 (MQ=255)
                                                                                                                                            cccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg           <  2:1290294/149‑1 (MQ=255)
                                                                                                                                             ccgccgCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAAATTAAACCCCCCCTGTTGTGTGGGGATAGCTGCTGGTGACAATGTGGAAAGGGGAGGGGGGGGGAGCAGCCCCAGAACCGCGCCAATTTAGTTTGcc         >  4:1463095/1‑151 (MQ=255)
                                                                                                                                                    aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa  <  2:848661/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357742‑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: