Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 367,573 G→A intergenic (‑63/+14) lacI ← / ← mhpR DNA‑binding transcriptional repressor LacI/DNA‑binding transcriptional activator MhpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913367,5730GA98.4% 181.8 / NA 61intergenic (‑63/+14)lacI/mhpRDNA‑binding transcriptional repressor LacI/DNA‑binding transcriptional activator MhpR
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (26/34);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCC  >  NC_000913/367329‑367814
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
gACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAACCGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTc                                                                                                                                                                                                                                             >  2:447276/1‑251 (MQ=255)
 aCTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCg                                                                                                                                                                                                                                            >  1:142548/1‑251 (MQ=255)
                 ttctggttccccacggtttgtaaTGTCATTCTGCTCCGCCATAGCCGTTTCCGCGTTGTCCCACGTTGTCGCATAAACGTGGCTGGCCTGATCCACCACGCGGGAACCAGTCCGATAAGCGACCCCGGCATACTCTGCGACATCGTATAACGTTACTGTTTTCTCATTACCCACCCTGGATTGTGTCGCGTCTGGGCGCTATCATGCCATCACC‑CGAAAGGCATTGCGCCACTCGGTGGTGGCATCGTaa                                                                                                                                                                                                                             <  2:286658/229‑1 (MQ=255)
                      gtgCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAatgca                                                                                                                                                                                                                        <  1:198262/250‑1 (MQ=255)
                        gCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATA‑GAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGcc                                                                                                                                                                                                                    >  1:403877/1‑251 (MQ=255)
                        gCCACGCGGTTGGGAATGTAATTCAGCTCCGCCACCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCCCCACCCTGAATTGACTCGCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAatgcatg                                                                                                                                                                                                                      <  2:557325/250‑1 (MQ=255)
                              cgGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAACGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTc                                                                                                                                                                                                               >  1:346785/1‑251 (MQ=255)
                               ggTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCg                                                                                                                                                                                                              >  1:470544/1‑251 (MQ=255)
                                    ggAATGTTATTCATCTCCGCCATCGTCGCTTCCCCTTTTTCCCGCTTTTGACCAGAAACGTGGCTGACCTGGTTCACCACTCGGGAAACGGTCTGATAAGAGACACCGTCATTCTCTGCGACATCGTATAACGTTACTTGTTTCACATTCACCACCCTGAATTTACTCTCTTCCTGGCGCTATCATGCCAT‑CCCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCttt                                                                                                                                                                                                         <  2:33856/251‑2 (MQ=255)
                                    ggAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTc                                                                                                                                                                                                         >  2:150882/1‑251 (MQ=255)
                                     gAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTcc                                                                                                                                                                                                        <  1:176229/251‑1 (MQ=255)
                                            aTTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATGGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGccacc                                                                                                                                                                                                 >  2:69234/1‑251 (MQ=255)
                                                  ccccGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCGGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACAGCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATa                                                                                                                                                                                           <  2:356850/249‑1 (MQ=255)
                                                         aTCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCg                                                                                                                                                                                    >  2:128850/1‑251 (MQ=255)
                                                           cgccgcTTCCACTTTTTCCCGCGTTTTCGCAGAAACGGGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCAGGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGAt                                                                                                                                                                                  <  2:599822/251‑1 (MQ=255)
                                                            accgtttCCACTTTTTCCCGCGTTTTAGCAGAAACGTGGCTGCACTGGTTCACCACTCGGGAAACGGTCTGCTAAGAAACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑CCCGCGAACGGTTTTGCACCATTCTATGGAGTCAACGCCAATGCATGCCGCTTCGCCTTCCTGCCACCTGCATAGCCTgagctt                                                                                                                                                                                 >  2:466723/6‑246 (MQ=255)
                                                                  tCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAAcc                                                                                                                                                                            <  2:604788/250‑1 (MQ=255)
                                                                     ccTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCCCCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCCCCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCctt                                                                                                                                                                        <  2:71791/250‑1 (MQ=255)
                                                                           tCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCg                                                                                                                                                                   >  1:93503/1‑250 (MQ=255)
                                                                             ccGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCTCCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCGCCAGAATAGCCTGCTATTCAACCCCTTCTTCGATc                                                                                                                                                                <  2:142548/251‑1 (MQ=255)
                                                                                   tttCGCAGAAACGTGGCTGGCATGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGtttt                                                                                                                                                          >  1:463514/1‑251 (MQ=255)
                                                                                      cGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTCTAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCt                                                                                                                                                       >  2:315105/1‑251 (MQ=255)
                                                                                        cAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTAc                                                                                                                                                     <  1:471457/251‑1 (MQ=255)
                                                                                           aaaCGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAAAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCg                                                                                                                                                  <  1:261886/251‑1 (MQ=255)
                                                                                                 ggctggcCGTGTTCGCCACGCGGGAAACGGTCTGATATGAGCCACCGGTAAACCCGGCTACATCGTATAACGGTACTGGTATCACATTCACCACCCTGAATTGCCACTCTTCCTTGCGCTATCAGGCCAT‑CCCGCGAAAGTTTTTGCACCAGTCTATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGTCACCAGAATAGCCTTCGATTCAACCCCTTCGTCGATTTGTTTTGCTACCCGTTGTa                                                                                                                                             <  2:449527/250‑1 (MQ=255)
                                                                                                             tCACCACGCGGGGAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAGCGTTACTGGTTTCCCATTCACCACCCTGAATTGACTCTCTTCCGAGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCTACCCCTTCTTCGATCTGTTTTGCTCCCCGTTGTATCGCCGGCAGATg                                                                                                                                >  2:92238/1‑251 (MQ=255)
                                                                                                             tCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGAt                                                                                                                                 <  2:210876/250‑1 (MQ=255)
                                                                                                               accacGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCt                                                                                                                              <  1:326380/251‑1 (MQ=255)
                                                                                                                ccacGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCtt                                                                                                                             >  1:168186/1‑251 (MQ=255)
                                                                                                                       ggAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCt                                                                                                                      >  1:37335/1‑251 (MQ=255)
                                                                                                                           aCGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCt                                                                                                                  <  1:99933/251‑1 (MQ=255)
                                                                                                                              gTCTGATAAGAGACACCGGCCTACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGtt                                                                                                               <  2:499333/251‑1 (MQ=255)
                                                                                                                                  gATAAGAGTCACCGTCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTGTTCCGCTGCCTGTTCAAt                                                                                                           <  1:428693/251‑1 (MQ=255)
                                                                                                                                          acacCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGCAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTg                                                                                                   >  1:532983/1‑251 (MQ=255)
                                                                                                                                           cacCTGCATACCCTTCGACAGCGTATAACGTTACTGGTTTCACATTCACCACCCTGGATTGACTCTTGTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTTTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCTTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGTAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTgc                                                                                                  <  2:144542/251‑1 (MQ=255)
                                                                                                                                               ggCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTc                                                                                              >  1:35631/1‑251 (MQ=255)
                                                                                                                                                    aCTCTGAGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTATCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATACATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCTCTCGCCat                                                                                         >  2:420718/1‑251 (MQ=255)
                                                                                                                                                     ctctGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCat                                                                                         >  1:158270/1‑250 (MQ=255)
                                                                                                                                                                     aaCGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGa                                                                         <  1:108991/250‑1 (MQ=255)
                                                                                                                                                                             tGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTCTCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATc                                                                >  1:374362/1‑251 (MQ=255)
                                                                                                                                                                              ggTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAAt                                                                 <  1:151566/249‑1 (MQ=255)
                                                                                                                                                                              ggTTTCACAGTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGGTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAAt                                                                 <  1:312902/249‑1 (MQ=255)
                                                                                                                                                                                 ttCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCAccc                                                            >  1:460552/1‑251 (MQ=255)
                                                                                                                                                                                  acacaTTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATAAACCAGATTCAGCCAGCCAATCACCCg                                                           >  2:470577/2‑251 (MQ=255)
                                                                                                                                                                                    acaTTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGATTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCAACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCgtt                                                         >  1:183093/1‑251 (MQ=255)
                                                                                                                                                                                               ccTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTgcgc                                              >  1:377248/1‑251 (MQ=255)
                                                                                                                                                                                                  gAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAg                                            <  1:149910/250‑1 (MQ=255)
                                                                                                                                                                                                      tGACTCTCTTCAGGGCGTTATCAGGCCAT‑ACCGCGAAAAGTTTTGCACCATTCGATGGTGGCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAAGGGTCATTGCGCTCGCCATATCCCCCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCgg                                         <  2:404302/249‑1 (MQ=255)
                                                                                                                                                                                                       gACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGTTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTGCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTGTTCCGCTGCCTGTTCAAGGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTa                                       <  1:298285/250‑1 (MQ=255)
                                                                                                                                                                                                        aCTCTCGTCCGGTCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAACAGCCTGCGATTCAACCCCTGCTTCGATCTGTTTTTCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACg                                     <  2:519485/251‑1 (MQ=255)
                                                                                                                                                                                                             ctTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGt                                                                                                                                                             <  2:116292/126‑1 (MQ=255)
                                                                                                                                                                                                             ctTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGt                                                                                                                                                             >  1:116292/1‑126 (MQ=255)
                                                                                                                                                                                                               tCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCAGATACACCAGATGCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATa                               <  2:534610/250‑1 (MQ=255)
                                                                                                                                                                                                                ccGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGAt                                >  2:628234/1‑248 (MQ=255)
                                                                                                                                                                                                                         aTCATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCATTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGAtc                      <  1:447276/249‑1 (MQ=255)
                                                                                                                                                                                                                           cATGCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGAtcttct                  <  1:150882/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCCAT‑ACCGCGAAAGGTTTTGCACCATTCGATTGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTctcc                <  1:129199/250‑1 (MQ=255)
                                                                                                                                                                                                                                cAT‑ACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCtcctcc             <  1:159012/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gAAAGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAAGAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCGTCTCCTCCTGATACCa     <  2:362112/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aagtttttGCACCATTCGATGGTGGCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCTTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGc   <  2:153874/247‑1 (MQ=255)
                                                                                                                                                                                                                                            aGGTTTTGCACCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGcc  <  1:12902/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
GACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAT‑ACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCC  >  NC_000913/367329‑367814

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