Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913540,5580TG33.3% 114.8 / 29.4 69intergenic (+50/‑7)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base T (9/37);  new base G (22/1);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 6.43e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTG  >  NC_000913/540344‑540795
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCAATCTCCGGGTAAGGGGGGTTTTTGGTCCATTTTTCAGt                                                                                                                                                                                                                         >  1:1418862/1‑238 (MQ=255)
  cAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                                                                                                                                                                                      <  1:1277070/239‑1 (MQ=255)
   aaTGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGGTAAGGGGGGGTTTTTTTTAATTTTTCAGttca                                                                                                                                                                                                                      >  1:1421280/1‑235 (MQ=255)
     tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGGTAAGGTGTGTTTTTGTTCAATTTTTCAGTccggc                                                                                                                                                                                                                    >  1:491132/1‑234 (MQ=255)
     tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTCCCGGTTTAAGGGTGGTTTTGTTTAAtttttccgtgagtg                                                                                                                                                                                                                    >  2:1341422/1‑228 (MQ=255)
                              gTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGGCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCAGTTTCCGGGTAAGGGGGGGTTATGTTCAATTTTTCAGTCCGCCGCGAAAGCCCGCTATCAGGGtt                                                                                                                                                                                             >  1:997821/1‑235 (MQ=255)
                                    ccttATGTTGGCCGTACCATTGGCGCGCGTAGCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCGTAAACATCAGCAGTAAGCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                                                                                                                                                                                          <  2:907317/229‑1 (MQ=255)
                                        aTGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGGGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCCGTTTTTTGGATTTCTGCGGTTTTTCt                                                                                                                                                                               >  1:598148/1‑239 (MQ=255)
                                             ggCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCCGGTGAGGGGTGTTTTTGTTCACTTTTTAAGTCAGCGGCTGGCGCCTGATTTCAGGATTTCGGGGGTTGTTCTTTAtt                                                                                                                                                                          >  2:1109988/1‑239 (MQ=255)
                                             ggCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGTGTGTTTTTGTTCAATTTTTCAGTCAGCCGCGGAAGCCTGTTATACGGGTTTTTGTGGGTtttttttttt                                                                                                                                                                          >  1:22738/1‑239 (MQ=255)
                                                 gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCCCCGGTTAAGGGGTGTTTTTGTTTAATTTTTCCGTCAGCCGCGAAAGCCCGTTATTCGGATTTCAGGGGtttttttttttttttt                                                                                                                                                                      >  1:1350765/1‑239 (MQ=255)
                                                  tACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGGTTAGGGGGGTTTTTGTTCAACTTTTTCGTCAGCCGCGAAAGGCTCTTATCAGGGGTTTTATGGttttttttttttttttt                                                                                                                                                                     >  2:620790/1‑238 (MQ=255)
                                                  tACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGTGTTTTTTTTCTATTTTGCAGGCCGCCGCGCAAGACGTTTTTCAGGCTATCATGGGTTTTTCTTTAttttttt                                                                                                                                                                     >  2:990832/1‑238 (MQ=255)
                                                   aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGCTAAGGGGTGTTTATGTTTAATTTTTCAGTCAGCCGCGGAAGCCTGGTTTTCGGATTTTGGGGGGTTTTTTTTAtttttttg                                                                                                                                                                    >  2:934973/1‑237 (MQ=255)
                                                   aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTCCCGGGTTAGGGGGGGTTATGTTTCATTTTTTAGGCAGCCGCCGAAACCCGTTAGTAGGATTTCAGTGGTGTTTTTTTCtttttttc                                                                                                                                                                    >  2:1420392/1‑237 (MQ=255)
                                                   aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTTCGGGTTTGGTGTGGTTTTTTTCAATTTTTTAGGGAGCCGCGAAAGCCTGGTTTTAGGGGTTTAGTGGGTtttttttttttttttc                                                                                                                                                                    >  2:429003/1‑237 (MQ=255)
                                                   aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTg                                                                                                                                                                     <  1:200364/238‑1 (MQ=255)
                                                   aCCATTGGCACGCGTATCACGAAAACCAGCTTACGTGGGGATGTGATTTACGACATTGACCAGAGCTTCCCTGTTGCGCCGCAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTACAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTAAGACGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGc                                                                                                                                                                    <  2:1023269/239‑1 (MQ=255)
                                                    ccAGTGGCGCGGGTAGCACGAAAACCATCTTACGTGGTGATGGGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAGTTTTTGCa                                                                                                                                                                   <  2:1283080/239‑1 (MQ=255)
                                                     cATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCaa                                                                                                                                                                  <  1:275824/239‑1 (MQ=255)
                                                       ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                                                <  1:479848/239‑1 (MQ=255)
                                                            gcgcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                                           <  2:167505/239‑1 (MQ=255)
                                                                  aTCACGAAAGCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTGAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACGCGGTTGTGGt                                                                                                                                                      <  2:868920/238‑1 (MQ=255)
                                                                                        gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTTTTTTCAATTTTGCAGTCCGCCCCGACAGCCTGTTATTCGGAGTTTAGTGGTTTTTCTTTATTTTTTGTAACACGAGTTTGGTTTCTCCCCCGCtccttcct                                                                                                                                   >  2:12592/1‑235 (MQ=255)
                                                                                         tgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCaa                                                                                                                                                                  <  1:93161/203‑1 (MQ=255)
                                                                                              tgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGATTAAGTGGGTTTTTGTTCAATTTTTCAGTCAGCCGCGCAAGCCTGTTATCAGGGTTTCAGTGGGTTTTTTTTATTTTTTGTAACACGGGTGTGGTTTCTCCTACGCTACTTTTTGCTTTTTAGt                                                                                                                         >  2:475207/1‑239 (MQ=255)
                                                                                                atTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTTCTGCTAAGGGGGGTTTATGTTCAAATTTTCCGGCCGCCCCGCAACCCTGGTATCAGGGTTTTATGGGGTTTTTTTTTCTTTTTGTCTCACAGTTTGGGTTCCTCCTACTCCCCTTTTTTCtttttattt                                                                                                                        >  2:874460/1‑233 (MQ=255)
                                                                                                atTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAGGTGGGGTTTTGTTCCATTTTTTAGTCAGCCGCGAAAGCCGGTTATCAGGGTTTCAGTGGGTTTTTTTTTTTTTTTTCAACACCGGTGTTGTTCTCCCCACCGTACTTTATTCTTTTccttt                                                                                                                        >  2:939026/1‑234 (MQ=255)
                                                                                                 tttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgtgt                                                                                                                                                                                                                                            <  1:981426/122‑1 (MQ=255)
                                                                                                 tttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgggtgt                                                                                                                                                                                                                                            >  2:981426/1‑122 (MQ=255)
                                                                                                 tttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTg                                                                                                                       <  1:1357278/238‑1 (MQ=255)
                                                                                                  ttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGcc                                                                                                                     <  1:934973/239‑1 (MQ=255)
                                                                                                   tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTTAGGGGTGTTTTTTTTTCATTTTTCGGTTAGCCGCGGAAGCCCGTTTTTCGGGTTTTTGTGGTTTTTCTTTTTTTTTTTTAACACGGgtgtt                                                                                                                                                        >  2:93161/1‑202 (MQ=255)
                                                                                                          ttGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTATCCGGTTAAGGGGTGTTTATGTTTAATTTTGTAGTCAGCCGCGCAAGCCTGTTCTCTGGGGTTTCTGGGGTTTTCTTTTTTTTTTTCAACACGGTTGTGGGTTCTCCTACGCTACTTTCTGTTTTTTTGTTGCCGCCAAATa                                                                                                             >  2:431379/1‑239 (MQ=255)
                                                                                                           tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCCGCTTAGGTGTGGTTATGTTTCATTTTGTAGTCAGCCGCGAAAGCCCGGTATCAGGATTTTAGTGGGTTTTCTTTTTTTTTTTTAACAACGGTGTGGGTCCTCCTACACTCCTTTCTTCTTTTTTTTTGCCGCaacata                                                                                                             >  2:1357278/1‑234 (MQ=255)
                                                                                                           tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAg                                                                                                            <  1:724814/239‑1 (MQ=255)
                                                                                                            gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAg                                                                                                               <  1:916945/235‑1 (MQ=255)
                                                                                                            gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGGCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                           <  1:1247387/239‑1 (MQ=255)
                                                                                                                       cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGttt                                                                                                                                                                                                                                          <  1:1092198/102‑1 (MQ=255)
                                                                                                                       cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGttt                                                                                                                                                                                                                                          >  2:1092198/1‑102 (MQ=255)
                                                                                                                        ccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTc                                                                                               <  1:431379/239‑1 (MQ=255)
                                                                                                                              gcgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCCCCGGGTTAGGGGTGGTTTTGTTCCATTTTTCCGTCAGCCGCGGAAGACTGGTTTCAGGGTTTTTGTGGGTTTTTTTTATTTTTTTTCACACGGGTGTGGGTCTTCCTACGCCACTTTTTGCTTTTTTGTTTTCGCCAAGTAGTCTGCTTACGCGCACGCaa                                                                                         >  2:492135/1‑239 (MQ=255)
                                                                                                                                 ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGTGTTTTTGTTCCATTTTGCGGTCAGCCGCGAAAGCCTGTTCTCTGGGTTTTAGTGGTTTTTTTTTTTTTTTTGCAAAACGGTTGTGGGGCCTCCTACTGTACTTTCTGTGTTTTTGTTTCCGCAAAGTAGCCTGCCTACGCTCACTCActat                                                                                      >  2:560802/1‑239 (MQ=255)
                                                                                                                                        gTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTg                                                                                                                                                                     <  1:444870/153‑1 (MQ=255)
                                                                                                                                              ttATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTa                                                                         <  1:560802/239‑1 (MQ=255)
                                                                                                                                              ttATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTa                                                                         <  1:990832/239‑1 (MQ=255)
                                                                                                                                                  ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCCGGTAAGGGGTGTTTTTGTTCAATTTTTCAGTCAGCCGCGGAAGGCCGGTATCAGGATTTTAGTGGTTTTTTTTTTTTTTTTTGAAAACgggt                                                                                                                                                           >  2:444870/1‑151 (MQ=255)
                                                                                                                                                  ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGGTAAGGGGGGTTTATGTTTAATTTTTCCGTCAGCCGCGGAAGCCCGCTATCAGGGTTTTTGTGGGTTTTTTTTATTTTTTGTTACACGGTTGTGGGTCCTCCTACGCCACTTTTTGCTTTTCTGTTGCCGCCAAGTAGCCTGCTTACGCTCCCGCAATATTCTTTTCTTGTTAcctc                                                                     >  2:871881/1‑237 (MQ=255)
                                                                                                                                                   cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                    <  1:871881/239‑1 (MQ=255)
                                                                                                                                                               cagTAATCTGGCCCCTGCAATGCCTGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGGCGGGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac                                                                                                                                                               <  1:308792/136‑1 (MQ=255)
                                                                                                                                                               cagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgtgt                                                                                                                                                                                                                                            <  1:574674/60‑1 (MQ=255)
                                                                                                                                                                  tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                                                                                                                                                                                 <  1:185762/84‑1 (MQ=255)
                                                                                                                                                                         gCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCGGGTAAGGGGGGGTTTTGtttca                                                                                                                                                                                                                                  >  2:574674/1‑57 (MQ=255)
                                                                                                                                                                         gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCt                                                                                                                                                                                                          >  1:724596/1‑84 (MQ=255)
                                                                                                                                                                         gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCt                                                                                                                                                                                                          <  2:724596/84‑1 (MQ=255)
                                                                                                                                                                           cccTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                                                <  1:774556/123‑1 (MQ=255)
                                                                                                                                                                            ccggCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGtttt                                           <  2:1359137/236‑1 (MQ=255)
                                                                                                                                                                            ccTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGTGTTTTTGTTCAATTTTTGAGGCAGCCGCGAAAGACTGtt                                                                                                                                                                                                       >  2:185762/1‑84 (MQ=255)
                                                                                                                                                                            ccTGCAATGCCCGTCCTGGCGGCGGGCATTCTCCGGTTAAGGTGTTTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGtttt                                           <  2:1421280/239‑1 (MQ=255)
                                                                                                                                                                                     cccGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGGGTTTTTGTTCCATTTTTCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTTTGTGGGTTTTTTTTTTTTTTTTCAACACGGTTGTGGt                                                                                                                                                      >  2:774556/1‑123 (MQ=255)
                                                                                                                                                                                        gTCCTTGCGGCGGGCATTCTCCGGTTATGGGGGGTTTTTGTTTAATTTTTCCGGCAGCCGCGAAAGCCTGTTTTCAGGATTTCTGTGGTTTTTTTTTTTTTTTTTCAACACGGGTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCTGTTTCCGCAAAGTATCTCGCTTACGCTCACGCAATATGCTTTTTTTGCTACCGCACACGCCTGCTTCGCTCTGGGGGTTTTCGGTCCtcct                                >  2:1200307/1‑236 (MQ=255)
                                                                                                                                                                                          ccTTGGGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                              <  1:686728/238‑1 (MQ=255)
                                                                                                                                                                                          ccTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                              <  1:367011/238‑1 (MQ=255)
                                                                                                                                                                                          ccTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                              <  1:1341422/238‑1 (MQ=255)
                                                                                                                                                                                          ccTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                              <  1:475207/238‑1 (MQ=255)
                                                                                                                                                                                           tttGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                            <  2:1304225/238‑1 (MQ=255)
                                                                                                                                                                                           gttGCGGCGGTCATTCTCCGGTTAAGTTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                            <  2:997821/238‑1 (MQ=255)
                                                                                                                                                                                           cTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                            <  1:492135/239‑1 (MQ=255)
                                                                                                                                                                                           cTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                            <  1:1424306/239‑1 (MQ=255)
                                                                                                                                                                                                     gCATTCTCCGGTTAAGGTGTTTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATgg                   <  2:37575/238‑1 (MQ=255)
                                                                                                                                                                                                              ggTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAgg         <  1:1200767/239‑1 (MQ=255)
                                                                                                                                                                                                                  aaGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGc     <  1:620790/239‑1 (MQ=255)
                                                                                                                                                                                                                      tgtttTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCtg  <  1:429003/238‑1 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
ACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTG  >  NC_000913/540344‑540795

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