Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5340CT28.3% 109.6 / 13.3 53R326K (AGA→AAA) lacIDNA‑binding transcriptional repressor LacI
Reads supporting (aligned to +/- strand):  ref base C (4/34);  new base T (15/0);  total (19/34)
Fisher's exact test for biased strand distribution p-value = 6.20e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACC  >  NC_000913/366304‑366782
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                
aTAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGGGGC‑A‑GGGGGGGTTTTTTTTTCACCCGTGAGA‑Cgggg                                                                                                                                                                                                                                             >  1:731956/1‑250 (MQ=255)
      gTTTCCTTTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGTTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGG‑CGGGCAACAg                                                                                                                                                                                                                                        <  2:365347/250‑1 (MQ=255)
      gTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAg                                                                                                                                                                                                                                        <  1:708657/250‑1 (MQ=255)
                    aTTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑Ac                                                                                                                                                                                                                        <  2:266953/251‑1 (MQ=255)
                      tGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑Acc                                                                                                                                                                                                                       <  1:459882/250‑1 (MQ=255)
                             cGCTCACAATGCCCCACAACATCCGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATCAGTGGGCTAACTCCCATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAGCATGTCGTGCCAGCTGCATTAATGAATCGGCAAACTCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTGCTTTTTTCCAGTGAGA‑CGGGAAAACGCTGATTGCCCTTC‑GTCGCCTGGcc                                                                                                                                                                                                               <  2:56883/251‑1 (MQ=255)
                              gCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGccc                                                                                                                                                                                                              <  1:356468/251‑1 (MQ=255)
                               cTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGt                                                                                                                                                                                                                                                       <  1:33789/211‑1 (MQ=255)
                                                      aGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGCGCGGCGGGTTTTGTATT‑TGGCGCC‑A‑GGGGGGGTTTTTTTTTCACCCTTGTGA‑CGGGGAACAGCTTGTTGCCCCTC‑ACCCCCCCGCCCTGCGAGAG‑G‑TGCCGCGAGCgggcc                                                                                                                                                                                    >  2:426561/1‑248 (MQ=255)
                                                       gCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATT‑GGGGGCCCG‑GGGGGGTTTTTTTTTTC‑CCCGTGAGAAGGGGCACCA‑CTGGTTGCCCTTC‑CCCCCCTGGCCCTGAGAGGG‑T‑TGCAGCAAACgggccc                                                                                                                                                                                   >  1:169566/1‑247 (MQ=255)
                                                         cGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGCGCC‑A‑GGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGGGAAAGCTGATTGGCCTTT‑ACCGCCCGGCCCCGAGAGAG‑T‑TGCAGCAAGCgggccccc                                                                                                                                                                                 >  1:325125/1‑245 (MQ=255)
                                                                aTAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGGGCC‑A‑GGGGGGGTTTTTTTTTCACCCGTGAGA‑CGGGCAACAACTGATTGTCCTTC‑ACCGCCTGGCCCCGAGAGAG‑T‑TGCAGCAAGCGGGCCCCGC‑Tgggtt                                                                                                                                                                         >  1:462561/1‑248 (MQ=255)
                                                                aTAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGttt                                                                                                                                                                         <  2:120836/251‑1 (MQ=255)
                                                                       gtAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTAGCCTT‑TGGCGCC‑A‑GGGGGGGTTTTTTTTTCCCCAATGGGA‑ACGGCAAAAACTGATTTCCCTTCCACC‑CCCGGGCCGGAGAGAG‑G‑TGCACGAAGCAGGCCACGC‑TGGGTTGtccccg                                                                                                                                                                  >  1:299824/1‑249 (MQ=255)
                                                                          aaGCCTGGGGTGCCTAATGATTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCcagcag                                                                                                                                                               <  2:521222/251‑1 (MQ=255)
                                                                                     gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTATT‑TGGCGCC‑C‑GGGGGGTTTTTTTTTTCACCAGGGAAAAGGGGAAACA‑CTGATTGCCCTTC‑CCCCCCTGGCCCCGAGAGAG‑G‑TGCAGCAAACGGTCCCCCC‑TGGGTTTCCCCCGCCGGCG‑GAAATccc                                                                                                                                                   >  1:719215/1‑250 (MQ=255)
                                                                                     gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGTGTATT‑TGGGGGC‑A‑GGGGGGGTTTTTTTTTTACCAGGGAGA‑CGGGCAAAAGCCGATTGCCCTTC‑ACCCCCTGGGCCTGGGAGAG‑GGTGC‑GCAAACGGGCCAAGC‑TGGTTTGCCCCCGCAgggggaaaaccc                                                                                                                                                    >  1:297695/1‑241 (MQ=255)
                                                                                         aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGTG‑AAAATCCTGtt                                                                                                                                                <  1:516633/250‑1 (MQ=255)
                                                                                         aaTGAGTGAGCTAACTCACATTAATTGCGTCGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCATCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAAGCCTGttt                                                                                                                                               <  2:546975/251‑1 (MQ=255)
                                                                                            gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATT‑TGGGGCC‑C‑GGGGGGGTTTTTTTTTCACCCGTGAGA‑CGGGCAACAGCTGATTGCCCCTC‑ACCCCCTGGCCCTTAGAGGG‑T‑TGCCGCAAGCGGGCCACCC‑CGGGTTTCCCCAGCAGGCGGAAAATC‑TGttttgt                                                                                                                                            >  1:499776/1‑248 (MQ=255)
                                                                                                 agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTACTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATT‑TGGGGCC‑AAGGGGGTTTTTTTTTTTCCCCA‑TGAGA‑CGGGGAACAGCTGATTGCCCCTC‑CCCCCCCGGCCCTGAGAGAG‑G‑TGCAGAAAACGGGCCACGC‑CGGTTTGCCCCCGCAGGCG‑GAAAACCTtttttgtggggg                                                                                                                                       >  1:94882/1‑243 (MQ=255)
                                                                                                  gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATggtgg‑t                                                                                                                                     <  1:467709/251‑1 (MQ=255)
                                                                                                   cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTgaga                                                                                                                                                                                                         <  1:195179/187‑1 (MQ=255)
                                                                                                      aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAAc                                                                                                                                 <  1:665367/251‑1 (MQ=255)
                                                                                                         cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATT‑TGGGGCC‑C‑GGGGGGGTTTTCTTTTTACCCGGGAGA‑CGGGCAACAGCTGGTTGGCCTTC‑ACCACCCGGCCCCGAGAGAG‑G‑TGCAGCAAACGGTCCCCGCCGGGTTTGCCCC‑GCAGGCG‑GAAAACCCGTTTGATGGGGG‑TTTAcggg                                                                                                                              >  1:201890/1‑250 (MQ=255)
                                                                                                           caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAAcggcgg                                                                                                                            <  2:128224/251‑1 (MQ=255)
                                                                                                              tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATT‑TGGCGCC‑A‑GGGGGGGTTTTTTTTTTACCAGTGAGA‑CGGGCAAAAACTGATTTCCCTTC‑ACCGCCTTGCCCTGGGGGAGGT‑TG‑AGCAAGCGGTCCACGCCGGGTTTGCCCC‑GCCGGCG‑GAAATCCTGttttttggggggtaaccggggggg                                                                                                                          >  1:393769/1‑230 (MQ=255)
                                                                                                               aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGat                                                                                                                         <  1:346422/250‑1 (MQ=255)
                                                                                                                aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGGCC‑‑‑CGGTGGTTTTTTTTTTTACCCGTGAGA‑CGGGCAAAAGCTGATTGTCCTTT‑ACCGCCCGGCCCTTTGGGAG‑G‑TGTAGCAAGGGGGCCAAGC‑CGGGTTGCCCCAGAAGGAG‑AAAAATCTGTTTTAtggtggggtacggcggggaat                                                                                                                        >  1:459766/1‑236 (MQ=255)
                                                                                                                         cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGc                                                                                                              <  1:485605/251‑1 (MQ=255)
                                                                                                                          gcTCTCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAATGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACAGGAGCt                                                                                                             <  1:104722/251‑1 (MQ=255)
                                                                                                                          gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGCGCC‑A‑GGGGGGGTTTTTTTTTTTCCCGTGGGA‑CGGGCAACAGCTGATTTCCCCTC‑ACCCCCTGGCCCTGAGGGAG‑G‑TGCAGCAAACGGGCCCCCC‑TGGTTTGGCCCAGCAGGCG‑GAAAACCTTTTTGATGGTGG‑GTTACGGCgggggataaaaaacgat                                                                                                             >  1:534853/1‑237 (MQ=255)
                                                                                                                          gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCt                                                                                                             <  2:169566/251‑1 (MQ=255)
                                                                                                                             cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGGGCC‑C‑GGGGGGGTTTTTTTTTCACCAGGGAGA‑CGGGCAACAGCCGATTGCCCTTC‑ACCGCCCGGCCCTGAGAGAG‑A‑TGCAGCAAGCGGTCCACGC‑TGGGTTGCCCCAGCAGGCG‑AAAAAACTGTTTGTTGGTGG‑TTAACGGCggggttttaaatgagatgtg                                                                                                          >  1:588742/1‑234 (MQ=255)
                                                                                                                             cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGt                                                                                                           <  2:357609/250‑1 (MQ=255)
                                                                                                                                 gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATT‑TGGGGCC‑C‑GGGGGGGTTTTTTTTTTACCAGTGAGA‑AGGGGCACAGCTGATTGCCCTTT‑CCCCCCTTGCCCTGAGAGAG‑G‑TGCAGCAAGCGGGCCACGC‑CGGGTTGCCCCCGCAGGGG‑GAAATCCTGTTTTATGGTGGGTTAA‑GGCggggggaaaaaatatgaggtcttt                                                                                                      >  1:318247/1‑230 (MQ=255)
                                                                                                                                 gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTcn                                                                                                      <  2:335676/251‑2 (MQ=255)
                                                                                                                                   ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGt                                                                                                    <  1:117527/251‑1 (MQ=255)
                                                                                                                                             gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTAt                                                                                          <  2:689517/251‑1 (MQ=255)
                                                                                                                                               cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATc                                                                                         <  1:660572/250‑1 (MQ=255)
                                                                                                                                                            gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGata                                                                           <  1:496902/251‑1 (MQ=255)
                                                                                                                                                                     gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTGGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCACAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCAcc                                                                   <  2:307676/250‑1 (MQ=255)
                                                                                                                                                                                         gcgcgGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCggn                                             <  2:325125/251‑2 (MQ=255)
                                                                                                                                                                                               ggAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCACCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTTCCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGgc                                       <  2:201890/251‑1 (MQ=255)
                                                                                                                                                                                                     gCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCAt                                 <  2:534853/251‑1 (MQ=255)
                                                                                                                                                                                                       ggTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGGGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGATTGATGGTGG‑TTAACGGCGGGGTATAACAATAGATGTCTTCCGTATCCTCGGATCCCACTACCCGGGTATCCTCACCCAAC‑CCCAGCCCGGGCTCGGGAAAGGCGCGCATTg                               >  2:105655/1‑251 (MQ=255)
                                                                                                                                                                                                            gCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGccc                          <  2:462561/251‑1 (MQ=255)
                                                                                                                                                                                                                    gggCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑C‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGccc                   >  1:216362/1‑250 (MQ=255)
                                                                                                                                                                                                                     ggCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCAt                  <  2:499776/251‑1 (MQ=255)
                                                                                                                                                                                                                            a‑GGGTGGTTTTTCTTTTCACCATTGAGA‑CGGGCACCAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTgcgcaca                         <  2:283412/238‑4 (MQ=255)
                                                                                                                                                                                                                            a‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCCCACC‑AACGCGCCGCCCGGACTCGGTAATGGCGCGCATTGCGcccc                         >  1:283412/1‑237 (MQ=255)
                                                                                                                                                                                                                                   gtTTTTCTTTTCACCAGTGAGA‑CGGGTAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGGCCAGCGCCATCTGATCGTTGGc      <  2:475660/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCTCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCACAGCGCCACCTGATCGTTGGc      <  1:512307/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tttcttttcACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCCTCTGATCGTTTGCAAc   >  1:599053/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tttcttttcACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAAAAAGAGCTGTCTTCGGGATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCCTCTGATCGTTGGCAAc   >  2:240295/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tcttttcACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAAcc  <  1:268972/250‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                
ATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCC‑A‑GGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAG‑T‑TGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCG‑AAAATCCTGTTTGATGGTGG‑TTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACC‑AACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACC  >  NC_000913/366304‑366782

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