Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913245,6620GA49.2% 6.0 / 95.9 65I48I (ATC→ATT) dpaApeptidoglycan meso‑diaminopimelic acid protein amidase A
Reads supporting (aligned to +/- strand):  ref base G (16/17);  new base A (11/21);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATAACCATGCGCACGGTCATAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTT  >  NC_000913/245424‑245912
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          
cATAACCATGCGCACGGTCATAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTa                                                                                                                                                                                                                                                 <  1:405603/250‑1 (MQ=255)
         gcgcACGGTCATAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGa                                                                                                                                                                                                                                        <  1:124388/250‑1 (MQ=255)
         gcgcACGGTCATAGGCAGTGGGGAAACCAATATTAATCGCTTTGTAGGAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTGTGGGCCCTAAGCCGCCGGAATATTTACAGATGTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGAt                                                                                                                                                                                                                                       <  1:217164/251‑1 (MQ=255)
                   aTAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAAct                                                                                                                                                                                                                             <  1:184960/251‑1 (MQ=255)
                          ttGGGGAAACCAATATTAATCGCTTTGTAGTAAAGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCtt                                                                                                                                                                                                                      >  1:409280/1‑251 (MQ=255)
                           tGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCCGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCtt                                                                                                                                                                                                                      >  1:26610/1‑250 (MQ=255)
                                  acctatatTAAGCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACTCTGGACGCTATATAACCCTTCCGGGCTTGTGAAATCGCCCTGACGCTGTTGTGGCCTTAAGCCTCCGGAAGATTTATAGATTTTATAGTTGTCGAGGAGTTGATATTTCTCGCCCTTTTTGACGTATAGATTGAGCGTACTTTCTTGCTTGATGATCTGGATGTAGACTGGGGATCCCATCAACTGCTGCTTATACTCtt                                                                                                                                                                                                              <  1:39191/247‑1 (MQ=255)
                                        atTAATCGCTTTGGAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTc                                                                                                                                                                                                         >  2:148987/1‑250 (MQ=255)
                                            aTCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAgg                                                                                                                                                                                                    <  1:90545/251‑1 (MQ=255)
                                               gCTTTGTAGTAACGGCTGTATGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTCCAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGCTGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGt                                                                                                                                                                                                 >  2:82603/1‑251 (MQ=255)
                                                cTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTg                                                                                                                                                                                                >  2:352463/1‑251 (MQ=255)
                                                  ttggagtaACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCACTTACGCTGTTTTGGCCCTAAGCCGCCGGACTATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGACGTGGt                                                                                                                                                                                              >  2:407865/5‑251 (MQ=255)
                                                  ttGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTg                                                                                                                                                                                                <  1:294484/249‑1 (MQ=255)
                                                       gtaACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTTTTTGCTCACAGGCGTGTTAGAGc                                                                                                                                                                                         <  2:319915/251‑1 (MQ=255)
                                                         aaCGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCt                                                                                                                                                                                        <  1:498298/250‑1 (MQ=255)
                                                           cGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTAct                                                                                                                                                                                     >  1:167511/1‑251 (MQ=255)
                                                              cTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTATAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTCTACTCTTTGCTCACAGGCGGGGAAGAGCTActgct                                                                                                                                                                                  >  2:3567/1‑251 (MQ=255)
                                                                   tggatttAAATGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATGGCCCTGACGCTGTTTTTGACCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTGGATATTGCTCGACCCTTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAACTCTGGATGTAGACAGGGGATCCAATCAACTGCTGCTTATACTCTTTGCTCACAGGCGAGGTAGGGCTACTGCTGccc                                                                                                                                                                              >  2:474627/5‑250 (MQ=255)
                                                                     gTTTTACCTGATTACGCTGGACGCTATAAAAACATTCCGTGCTTTTGAAATCGCCCTGCCGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGACCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCACCTGCTGCTTATAATCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGc                                                                                                                                                                           >  2:52411/1‑251 (MQ=255)
                                                                                 tACGCTGGACGCTATAAATCCCTTCCGGGCTTGTTAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTAGAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCaaa                                                                                                                                                                <  1:506527/250‑1 (MQ=255)
                                                                                      tGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAAc                                                                                                                                                          >  2:37531/1‑251 (MQ=255)
                                                                                           gCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACg                                                                                                                                                     <  1:313324/251‑1 (MQ=255)
                                                                                                 aaaCCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATca                                                                                                                                               <  2:137542/251‑1 (MQ=255)
                                                                                                 aaaCCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATca                                                                                                                                               <  1:508532/251‑1 (MQ=255)
                                                                                                    cccTTCCGGGCGTTTGAAATCGCCCTGACGCTGTTTTGGCGCTAAGCCGCCGGAATAATTACAGTTTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGATATCGAGCGTACGTTCTGCCTTGAAGAGCTGGATGTAGACAGGGGATCCCAGCATCTGCTGCTGATACGCTTTGCTCACAGGCGTCGTAGAGTGACTGCAGCCCAGCAAACCTGCAAATGAAACTCACTTGATcaaca                                                                                                                                            <  2:226529/251‑1 (MQ=255)
                                                                                                        tCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAacaaca                                                                                                                                         <  1:284680/250‑1 (MQ=255)
                                                                                                           gggCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACCGCTTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCg                                                                                                                                                                     >  2:211683/1‑219 (MQ=255)
                                                                                                           gggCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCg                                                                                                                                                                     <  1:211683/219‑1 (MQ=255)
                                                                                                                ttggAAATCGCCCTGACGCTGTTTTTGCCCTAAGTCGCCGTAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTTCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTTGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGa                                                                                                                                 <  1:272254/247‑1 (MQ=255)
                                                                                                                     aaTCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAAt                                                                                                                           >  1:374034/1‑251 (MQ=255)
                                                                                                                      aTCGCCCTGACGCTGGTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATg                                                                                                                          <  2:420961/251‑1 (MQ=255)
                                                                                                                         gCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAACCGCACGGGATCAACAACATCGCAAGAATTAATGCg                                                                                                                        >  1:318497/1‑250 (MQ=255)
                                                                                                                          cccTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGAt                                                                                                                      <  1:214188/251‑1 (MQ=255)
                                                                                                                           ccTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGAGCAACAACATCGCAAGAATTAATGCGAtt                                                                                                                     <  2:318497/251‑1 (MQ=255)
                                                                                                                              gACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTa                                                                                                                  <  1:98242/251‑1 (MQ=255)
                                                                                                                                 gCTGTTTTGGCCCTAAGCCTCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACATGGGATCCCATCAACTGCTGCTTAGACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGc                                                                                                               <  1:332723/251‑1 (MQ=255)
                                                                                                                                  cTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGc                                                                                                               >  1:122960/1‑250 (MQ=255)
                                                                                                                                    gttgtGGCCCTAAGCCGCCTGAATATTTACAGATTTTATAGCTGTCGAGCATTTGATAGTGCTCGCCCATTTTGACGTATAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTTGTATAGCTACTGCTGCCCAGCAAACCGGCAAATGAATCTCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATa                                                                                                            <  1:265410/247‑1 (MQ=255)
                                                                                                                                               aaGCCGCCGTAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCGTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTc                                                                                                  <  1:318226/250‑1 (MQ=255)
                                                                                                                                               aaGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTAtt                                                                                                   <  1:315007/249‑1 (MQ=255)
                                                                                                                                                       ggAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATaaa                                                                                         <  2:409280/251‑1 (MQ=255)
                                                                                                                                                        gAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATaaa                                                                                         >  1:27909/1‑250 (MQ=255)
                                                                                                                                                         aaTATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAAc                                                                                       >  1:413269/1‑251 (MQ=255)
                                                                                                                                                          atatTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGAGCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTGTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACg                                                                                      <  2:501853/251‑1 (MQ=255)
                                                                                                                                                            atTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGt                                                                                    >  2:46775/1‑251 (MQ=255)
                                                                                                                                                             tttACAGATTTTATAGCTGTCGAGCAGTTGCTATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGtt                                                                                   >  1:398452/1‑251 (MQ=255)
                                                                                                                                                               tACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTa                                                                                  <  1:471949/250‑1 (MQ=255)
                                                                                                                                                                  aGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTacac                                                                               >  2:212050/1‑250 (MQ=255)
                                                                                                                                                                   gATTTTATAGCTGTCGAGCAGTTGATGTTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGCTCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTacacac                                                                             >  2:49188/1‑251 (MQ=255)
                                                                                                                                                                        tataGCTGTCGAGCAGTTGATATTGCTCGCCCGTTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGGAGAGCTACTGCTGCCCAGCAAACCGGCAAATTAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATGCCTTGTTAAAACGGTTACACACGCCa                                                                         <  2:88911/250‑1 (MQ=255)
                                                                                                                                                                               ggcGAGCAGTTGTTATTGCTCGCCCATGTTGACGTAGAGATCGAGCGTACGTTCGTCCGTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCCCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACgg                                                                   <  2:327102/247‑1 (MQ=255)
                                                                                                                                                                               gTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCa                                                                 <  1:82603/251‑1 (MQ=255)
                                                                                                                                                                               gTCGAGCAGGTGATATGGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGGTCTTCCTTGAAAATCAGGATGTAGACAGGGGTGCCCATCAACTGCTGCTTATACTCTTTGCTCGCAGGCGTTGTAGAGCTACTGCTGCCCATCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAAGTAATGCGATTTTACGCATACTGCTTATTCCTTGATATAACGGTTACACACGCCAGGACga                                                                   <  2:515741/249‑2 (MQ=255)
                                                                                                                                                                                                       catgttgactgagagATCGAGCGGACGTACTTCGTGGAAAATCTGGATGTAGACAGGGGATCCCCTCAACTGCTGCTTAGACGCTTTGCTCACAGGCGTGGTAGTGCTACTGCTGACCAGCAAACCGGCAAATGAAACGGACGGGATCAACAACAGCGCAATAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACTCGCCAGGACGGCAAAATGGATACCAAATCGGAATAGt                                         <  2:223139/240‑1 (MQ=255)
                                                                                                                                                                                                               cGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAAAATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGAAGGCAAAATGGATCCCAAATCGGAATAGTCTGGGttt                                 >  1:375239/1‑251 (MQ=255)
                                                                                                                                                                                                                  agagATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGCAGACAGGCGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACAGGCAAATGAAACGCACGGGCTCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGAAATAGGCTGGGTTTgg                               >  2:445773/1‑250 (MQ=255)
                                                                                                                                                                                                                  agagATCGAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATCCTCTTTGCTCACAGGCGTTGTAGAGCTACTGCTGCCCAGCCAACCGGCAAATGAAACCCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACTCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAACAGTCTGGATTTGGa                              <  1:220422/251‑1 (MQ=255)
                                                                                                                                                                                                                        cGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTc                        >  2:233716/1‑251 (MQ=255)
                                                                                                                                                                                                                         gAGCGTACGTTCTTCCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCAAGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTAATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGAATTTGGAAGGCTc                        >  1:179459/1‑250 (MQ=255)
                                                                                                                                                                                                                           gCGTACGTTCTTGCTTGAAGAGCTGGATGTAGACAGGGGATCCCCTCAACTGCTGCTTATACTCGTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACTTCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCGGGACGGGAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTcaca                     <  1:494089/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCa               <  1:212050/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCa               <  1:32495/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tcCTTGATAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCaa           >  2:164874/1‑250 (MQ=255)
                                                                                                                                                                                                                                      tcCTTGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCaa           >  2:421500/1‑250 (MQ=255)
                                                                                                                                                                                                                                          tGAAAATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAaga       <  1:421500/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTAATGCTGCCCAGCAAACAGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAACAGAGtttt  >  1:339586/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          
CATAACCATGCGCACGGTCATAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCGGGCTTTTGAAATCGCCCTGACGCTGTTTTGGCCCTAAGCCGCCGGAATATTTACAGATTTTATAGCTGTCGAGCAGTTGATATTGCTCGCCCATTTTGACGTAGAGATCGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTT  >  NC_000913/245424‑245912

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