Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA46.2% 58.5 / 52.3 78L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (32/10);  new base A (0/36);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 1.90e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039
                                                                                                                                                       |                                                                                                                                                        
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                         >  2:2026022/1‑150 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       <  1:1243860/151‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       <  1:1848082/151‑1 (MQ=255)
  tAAGCGTAGCGCTTCAGGCAATTCCAGCCGCGGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACACCCGCGCATTCATCCAACGCCATGTGCAAAAATGCATTCGCGTCGGCTGTCTGCCAGCTATTTTCCGCCGcaacaa                                                                                                                                                        <  2:1237731/148‑1 (MQ=14)
  tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       >  2:768295/1‑150 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                      <  2:1258916/150‑1 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                     >  1:1792281/1‑151 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                    >  2:352727/1‑151 (MQ=255)
       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                                 >  1:727587/1‑151 (MQ=255)
       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgac                                                                                                                                                 >  1:1616099/1‑149 (MQ=255)
          gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                                                                              >  1:1635819/1‑151 (MQ=255)
            gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                            >  1:668219/1‑151 (MQ=255)
             cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                           >  1:1114480/1‑151 (MQ=255)
             cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                           >  1:1683744/1‑151 (MQ=255)
                cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                        >  2:1068227/1‑151 (MQ=255)
                      aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                                  >  1:281082/1‑151 (MQ=255)
                      aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                                   >  2:1875817/1‑150 (MQ=255)
                       tccccgCCGCAGCCCTGTGTCGGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCAACGCCTTGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagta                                                                                                                                    <  1:1319676/143‑1 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGAGGCTAACGTGATTCATTTCCGAAAACAACCGCGCATTTATCCAACGCCATGTGCAAAACTGCCTTCGCGGACGCTGTCTGCCAGCTATTTTTCAGCAGCAACAACACCGCCGTTCACTCCAGTAGTGGc                                                                                                                               >  2:2139920/1‑150 (MQ=255)
                                     cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                         >  2:1700009/1‑143 (MQ=255)
                                      tgtgCCGGCGCCTCCCGTGATTATTTCCCCCCAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAATGCCTCCGCGGCGGTGGTCGGCCGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                                       <  2:1704582/145‑1 (MQ=255)
                                           cggcgcCTCCCGTCTTTCATTCCCCCCAACACCCCCGCATTTACCCACCCCGATGGGAAAAAGGGCCTCGGCGCCGGCTCTCTCCCGTTTATTTCCCGGCCGCACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                            <  1:1490716/147‑1 (MQ=255)
                                            aGCGGCTACCGTGAATCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGAGGCTGTCTGCCAGCTCTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggggaagca                                                                                                                 >  2:1358711/1‑136 (MQ=25)
                                             gcgcCTCCCGTGTTTCATTCCCCCCAACACCCGCCCTTTTACCCACCGCCATGTGAAAAAAGGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                                        <  1:36197/151‑1 (MQ=255)
                                             gCGGCTCCCGTGTTTTTTTCCCCCAAACACCCGCGCTTTTACCCAACGCCATGTAAAAAAATCCCTTCGGGGGGGCTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                                          <  2:2131483/149‑1 (MQ=255)
                                               ggCTCCCTTGATTCATTCCCGCCAACACCCGCGTTTTTACCCACCCCATTGTGAAAAAATGCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                                          <  1:420232/147‑1 (MQ=255)
                                               ggCTCCCGTATTTCTTTCCCGCCAACACCCGCGCTTTTATCCACCCCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                                       <  1:640800/150‑1 (MQ=255)
                                                gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggtgccggggaata                                                                                                               >  2:1180879/1‑131 (MQ=255)
                                                 atACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTAcgagggaaagatct                                                                                                         >  2:1299634/2‑137 (MQ=25)
                                                    cggTGATTCATCCCCCCCAACACCCGCGAATTAATCCACCGCATTTAGAAAAATTGCCTCCGCGGCGGCTGTCCCCCAGCTTTTTTCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                                  <  1:1910652/148‑1 (MQ=255)
                                                    ccGTGTTTCATTCCCGCCAACAACCGCCATTTTACCCAACCCCTTGTAAAAATTTGCCTCCGGGGCGCCTGTTTCCCAGCTTTTTCCCCCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                                  <  1:2177338/149‑1 (MQ=255)
                                                    ccGTGTTTCTTTCCCCCAAAAACCCGCCCTTTTACCCACCGCCATGTCAAAAAATCCCTTCGCGGCGGTTTTCTCCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                                     <  1:972508/147‑1 (MQ=255)
                                                      gtttttCTTTCCCCCCACCACCCGCGCTTTTATCCACCCCCTTGTGAAAAAAGGCCTTGGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                                <  1:482270/146‑1 (MQ=255)
                                                       tgttcctttccccccAACACCCGCCCTTTCCCCCACCGCCTTGTGAAAAAGCGCCTGCGCGGCGGTTTTCGCCCAGTTTTTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                                <  2:1097151/136‑1 (MQ=255)
                                                        gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTaatacc                                                                                                  >  1:146954/1‑142 (MQ=25)
                                                              tttCCCCCAACACCCGCCCTTTTCACCACCCCCATGTGAAAAAATCCCTTCGGGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                                         <  2:560526/149‑1 (MQ=255)
                                                              ttccccccAACACCCCCCCTTTTACCCACCGCCATGCGAAAAAATCCCTTGGCGGCGGCTGTCGGCCAGTTATTTTGCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                                <  2:1653482/136‑1 (MQ=255)
                                                               ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACTTTTAACCCCCCCTGTtttgt                                                                                      >  1:267943/1‑151 (MQ=25)
                                                                 cccccaccaCCCGCGCATTTATCCACCGCCTTGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                                                      <  1:464253/146‑1 (MQ=255)
                                                                      aacaccCGCCCTTTTACCCACCGCCTTGTGAAAAAAGCCCTTCGCGGCGGTTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                                  <  1:1032545/148‑1 (MQ=255)
                                                                         aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACTTTAAACCCGTCATGTTTTGTGGCaa                                                                                 >  1:382756/1‑146 (MQ=255)
                                                                             cgcgCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAGGGAAAAGACTTTAAACCCTCATGGTGTGTggggttag                                                                              >  2:333337/1‑140 (MQ=25)
                                                                              gccaaTTTACCCAACCCCTTGTCAAAAATCCCCTTCGGGGCGGTTTTCCGCCGTCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                         <  1:1508304/146‑1 (MQ=255)
                                                                                   tttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGAACGCTGCTgg                                                                       >  1:787669/1‑147 (MQ=21)
                                                                                   tttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAAACCCCCATGTTTTGTGGGGAAAGCTGCTgg                                                                       >  1:465384/1‑147 (MQ=25)
                                                                                   tttATCCAACCCCATGTAAAAAAATCCCTTCGCGGCGGCTGTCTCCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                                       <  1:342028/147‑1 (MQ=255)
                                                                                   tttACCCACCGCCATGTCAAAAAATCCCTCCGGGGCGGTTGTCTGCCAGTTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                        <  1:904170/146‑1 (MQ=255)
                                                                                   tttACCCACCCCCATGTGAAAAAATGCTTCCGCGGCGGCTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                                     <  1:1969844/149‑1 (MQ=255)
                                                                                   tttACCCAACGCCATGTGAAAAATTCCCTCGGGGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                         <  1:962394/145‑1 (MQ=255)
                                                                                        ccAACCCCATGTGCAAAAATCCCTCCGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                              <  1:1429417/151‑1 (MQ=255)
                                                                                         caccccaagggcaaaaaaTCCCTCCGGGGGGGCTTTCTCCCAGGTTTTTTCCCCCACCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                                  <  2:721553/133‑1 (MQ=255)
                                                                                         cAACCCCATGTGAAAAAAGTCCTTCCGGGCGTGTTGCCCCCATGTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                                  <  2:2177442/146‑1 (MQ=255)
                                                                                           aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTAAACCCCGCCTGTTTTGTGGGGAACGCTGCTGGGAaaaaaaggg                                                            >  1:887493/1‑145 (MQ=25)
                                                                                            cccctggtaaaaaaaTGCCTCCGGGGGGGTTGCCTGCAAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                                <  2:468096/135‑1 (MQ=255)
                                                                                            cGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAAAATTATAACCCCCCTCTTGTGTGGGGAACGCTGCTGGTGACaaaatgg                                                            >  1:1242600/1‑144 (MQ=21)
                                                                                              ccttgtaaaaaatTGCCTGCGGGCGGGCTGCCGGCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                             <  2:626541/138‑1 (MQ=255)
                                                                                               cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTATCCCCTCCTGTTGTGTGGGGAACGCTGCTGGTGACAAtgtg                                                             >  1:1613105/1‑145 (MQ=37)
                                                                                                    gaaaaaaaGCCTTCCGGGGCGGCTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCGTGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                                   <  2:2076678/149‑1 (MQ=255)
                                                                                                       aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACCTTTAAGCCCCCCTGTTGTTTTGGGATCCCTGCTGGTGAAAAATTGGGAAGGGGAggggg                                               >  1:810743/1‑147 (MQ=255)
                                                                                                            tcccTTCGCGGGGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                                            <  1:1498079/148‑1 (MQ=255)
                                                                                                              ccTCCGCGGCGCCTTTCGCCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                         <  2:434068/151‑1 (MQ=255)
                                                                                                               cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAGTTTTAACCCCCCCAGTGTTGTGGGGATAGCTGCTTGGGAAAAAGTGGGGAGGGAAGAgcgg                                              >  1:2049355/1‑145 (MQ=255)
                                                                                                                   cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCCCTTTTTTGTGGGGATACCTGCTGGGTACAAAGTGGGAAggggaaggggggggatc                                        >  1:307131/1‑134 (MQ=255)
                                                                                                                      ggcggcTGTCCCCCGGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                         <  2:1993011/144‑1 (MQ=255)
                                                                                                                       gcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACTGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  2:89198/151‑1 (MQ=255)
                                                                                                                            tttcgccgagtttttttcCCGCCAAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  2:1370889/133‑1 (MQ=255)
                                                                                                                               ctcccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGTGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                              <  2:506444/143‑1 (MQ=255)
                                                                                                                                  ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                          <  1:1550277/147‑1 (MQ=255)
                                                                                                                                  ccAGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                       <  1:932212/150‑1 (MQ=255)
                                                                                                                                  ccAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  1:1800706/140‑1 (MQ=255)
                                                                                                                                  ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGAAAAAAACTTAAACCCCCCCCGTTCTGTGGCGATCGCTCCGGGGGACCAAGTGGGAAGGGGAGGGCGGCGGAACAGCTccccaaacccgc                        >  2:814227/1‑139 (MQ=255)
                                                                                                                                         tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAACCCCCCCTTTTTTGTGGGGATCGCTGCTGGGAAAAAAGTGGAAAAGGAAGGGGGGGGGAACAGCCCCCCAACCGCGCCAATTTAg                >  1:135170/1‑150 (MQ=255)
                                                                                                                                         ttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                            <  2:1308814/138‑1 (MQ=255)
                                                                                                                                         ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa                 <  1:395976/149‑1 (MQ=255)
                                                                                                                                          tttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt               <  1:1516533/150‑1 (MQ=255)
                                                                                                                                             cccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc             <  2:748908/148‑1 (MQ=255)
                                                                                                                                              ccgccgAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGc           <  2:281047/150‑1 (MQ=255)
                                                                                                                                              cccccacaaaaaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc             <  1:1950260/142‑1 (MQ=255)
                                                                                                                                              ccacggcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                          <  2:1432384/130‑1 (MQ=255)
                                                                                                                                                  cgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAACCCGCCTGTTGTGTGGGAATTGCTGCTGGTAAAAAATTGGGAAGGGGGGGGGGGGGGATTCACTCCCGAACCCCGCTAATTTAGTTTGCCt         >  1:1069530/1‑148 (MQ=255)
                                                                                                                                                  cgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt         <  1:1426215/148‑1 (MQ=255)
                                                                                                                                                    caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGAAAAAAAAATTTAACCCCTCCTTGTTTTTGGCGCACGCCGCTGGTGAAAATGTGGGAAGGGGAGGGGGGGGGAACAGCCCCCGGAACCCGCCAAATGAGGTCGCCTccc      >  2:1498248/1‑147 (MQ=255)
                                                                                                                                                      acaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac  <  1:1056026/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                        
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039

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