Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC54.1% 47.4 / 91.5 61K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (26/2);  new base C (0/33);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 3.10e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-13

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357730‑358026
                                                                                                                                                         |                                                                                                                                                     
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                       <  1:265689/151‑1 (MQ=255)
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                       <  1:346958/151‑1 (MQ=255)
 cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                     >  1:911970/1‑151 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                   >  1:750579/1‑151 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                   >  2:903222/1‑151 (MQ=255)
           tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                            >  1:648115/1‑150 (MQ=255)
                   gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                   <  2:112854/151‑1 (MQ=255)
                   gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                   >  2:1112945/1‑151 (MQ=255)
                    tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                  >  2:332735/1‑151 (MQ=255)
                        gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                              >  2:329810/1‑151 (MQ=255)
                            cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                          >  2:1059497/1‑151 (MQ=255)
                              ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                          >  2:194929/1‑149 (MQ=255)
                                 aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                       >  2:1234288/1‑149 (MQ=255)
                                 aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                       >  2:6552/1‑149 (MQ=255)
                                 aaTCCCAGCCCCGGACCTGTGTCGGCGGTTCCCGTATTTCATTCCCGCCACCAACCGCGCTTTACCCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAAccgccg                                                                                                                                  <  1:833313/138‑1 (MQ=14)
                                  aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    >  1:933705/1‑151 (MQ=255)
                                  aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                     >  2:214672/1‑150 (MQ=255)
                                   ttCCAGCCGCAGACCTGTGTCGGCGGCTCCCTTGTTTCTTTCCCCCCAACACCCGCCCTTTTACCCACCCCCATGTGCAAAATTCCCTCCGCGGCGGCTGTCTCCCGTTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                      <  2:1194886/148‑1 (MQ=255)
                                   ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   >  1:696447/1‑151 (MQ=255)
                                   tcccggCCGCGGCCTTGTGCCAGCGGCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTCAAAAAATGCCTCCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    <  2:1069745/145‑1 (MQ=255)
                                    tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                  >  2:916307/1‑151 (MQ=255)
                                          cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                            >  1:961136/1‑150 (MQ=255)
                                             agcctgtggccggCGGCTACCGTTTTTCTTTCCCCCCACCACCCGCCCATTTACCAACCCCCAGGGCAAAAAATCCCTCCGGGGGGCTTGTCGCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  1:986350/137‑1 (MQ=255)
                                             agcccttttccAGCGCCTCCGGTGTTTCTTTCCCGCCACCACCCGCGCTTTTACCCACCCCCATGGGAAAAAATGCCTTGGGGGGGGTTGTCTGCCGTTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  2:193202/139‑1 (MQ=255)
                                             agcccTGTGTCGGCGGCTACCGTGTTTCATTCCCCCCAACAACCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                         <  2:402876/148‑1 (MQ=255)
                                             aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                             >  1:570104/1‑147 (MQ=255)
                                                ccTGTGTAGGCGGCTACCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                    <  2:175879/151‑1 (MQ=21)
                                                     gtCGGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGTATTTACCCACCGCCATGGGAAAAAATCCTTGCGGGGCGGCTGTCGGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                 <  2:1160199/149‑1 (MQ=255)
                                                     gccgcgggcacccGTAATTCATTCCCCCCAACACCCGCCATTTTACCCAACCCCTTGTAAAAAAATCCTTCCGCGGGGCCTGCCTCCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                   <  1:567842/136‑1 (MQ=255)
                                                           gtctcccgtttttCTTTCCCCCCAAAACCCGCCCTTTTCCCAAACCCCTTGTGAAAAATTCCCTGCGGGGGGGCTGTCTCCAGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                               <  1:589440/134‑1 (MQ=255)
                                                           ggCTACGGTGATTCATTCCCGCCAACACCCGCTCATTTACCCAACCCCATGTGAAAAAATGCCTTCGGGCCGGTTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                           <  2:761641/149‑1 (MQ=255)
                                                           ggCTACCGTGATTCTTTCCCCCCAACACCCGCGCTTTTATCCAACGCCATGTACAAAAATGCCTGCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                         <  2:1205044/151‑1 (MQ=255)
                                                              tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAATaaaatt                                                                                           >  2:481212/1‑141 (MQ=255)
                                                               aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaaaatttaaccc                                                                                     >  2:195228/1‑138 (MQ=25)
                                                                 ccgtttttctttCCCCCAACCACCCGCCCTTTTACCCAACGCCATGGAAAAAATGCCCTTCGGGGGGGTTTTCTGCCGGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                    <  1:181566/126‑1 (MQ=255)
                                                                    gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAAATTTAAGcc                                                                                      >  1:814665/1‑146 (MQ=25)
                                                                      ttatttcccGCAAAAAACCGCGAATTTACCCACCGCCATGTACAAAAATGCTTCCGGGGCGGTTTTTTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                     <  2:133435/137‑1 (MQ=255)
                                                                      ttCTTTCCCGCCAACACCCGCCTTTTTACCCACCCCCTTTTAAAAAAATGCCTCCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                               <  2:872545/150‑1 (MQ=255)
                                                                      tactttccccccAACACCCCCCCTTTTCCCCACCGCCATTTAAAAAAGTCCCTCCGCGGGGGTTGCCTGCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                         <  1:158964/130‑1 (MQ=255)
                                                                             cccccAAAACCCGCGCTTTTACCCACCGCATTGTAAAAAAATGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                        <  2:1006117/147‑1 (MQ=255)
                                                                                 caacaaCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGGAAAAAATTTTAACCCGTCACGGTTTGTgggg                                                                     >  1:533778/1‑147 (MQ=21)
                                                                                     cacccGCCCTTTTACCCACCCCCATGTAAAAAATTGCCTCCGGGGCGCTGGTCTCCCATTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                  <  1:573724/146‑1 (MQ=255)
                                                                                        ccgcgtttttccccAACCCCTTGTGAAAAAATCCCTCCGGGGCGGCTTTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                             <  1:234971/139‑1 (MQ=255)
                                                                                               tttaccaaccGCCGTGTGAAAAAATCCCTGCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                         <  1:632244/139‑1 (MQ=255)
                                                                                               tttACCCAACGCCTTGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                        <  2:809712/149‑1 (MQ=255)
                                                                                                ttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACTTTAAATCCGCCCTGTTTTGTGGCGATCg                                                                 >  1:808891/1‑139 (MQ=37)
                                                                                                      acccccttgtgAAAAAATCCCTCCGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                <  2:481975/143‑1 (MQ=255)
                                                                                                      acccccATGGGAAAAAAGGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                <  2:163306/146‑1 (MQ=255)
                                                                                                      accGCAATGTGAAAAAAGGCCTTCGCGGCGGCTGTCTGCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                   <  1:1054077/145‑1 (MQ=255)
                                                                                                      aaCCCCATGTACAAAAATGCCTTCGCGGCGGTTGTCGGCAAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                               <  1:239940/151‑1 (MQ=255)
                                                                                                        ccctaTGTGAAAAAATGCCTCCGCGGGGCCTGTCTCCAAGTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                               <  1:884427/145‑1 (MQ=255)
                                                                                                          ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAAGTTTAAACCCTCCTGTTTTTTGGGGATCGCCGCTGGGGACAAAGTgg                                               >  1:194294/1‑147 (MQ=19)
                                                                                                          ccATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCGGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                               <  2:278012/147‑1 (MQ=255)
                                                                                                          ccATGTAAAAAAATGCCTCGGGGGCGGCTTTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                   <  2:1037713/143‑1 (MQ=255)
                                                                                                           cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTAAACCCCCCCTGTTTTGTGGGGATTGGCGCTGGGGACAAAGTGGgaa                                            >  2:967400/1‑149 (MQ=21)
                                                                                                             tgtgAAAAATTCCCTCCGCGGCGGCTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                         <  2:199620/150‑1 (MQ=255)
                                                                                                                   aaaatTGCCTGCGCGGCGGTCGTCGTCCAGTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                      <  2:350674/148‑1 (MQ=255)
                                                                                                                     aaaTGCCTTCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                 <  2:922042/151‑1 (MQ=255)
                                                                                                                     aaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAATTTAAAACCCCCCTGTTTTGTGGGGAACGCTGCTGGGAACAATGTGGAAAGGGAAgggggg                                 >  2:629566/1‑146 (MQ=255)
                                                                                                                           cctcgggggcggcTGTCGCCCTTTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                            <  1:55694/144‑1 (MQ=255)
                                                                                                                           ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAAATTTTAACCCCCCCTGTTTTGTTGGGATCGCTGCTGGGGAAAAAGTGGGAAGGGGAggggggggg                              >  2:956126/1‑141 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  2:957230/147‑1 (MQ=255)
                                                                                                                                                   tATTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg     <  1:812165/150‑1 (MQ=255)
                                                                                                                                                      tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  2:183117/150‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357730‑358026

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