Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,245,221 C→G V517L (GTA→CTA)  cirA ← iron‑catecholate outer membrane transporter CirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,245,2210CG98.2% 194.5 / NA 56V517L (GTA→CTA) cirAiron‑catecholate outer membrane transporter CirA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base G (27/27);  minor base T (0/1);  total (27/29)
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.97e-01

CTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACAATCTTACATGCGCCACGGC  >  NC_000913/2244976‑2245466
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
cTGCCCGGTATAGTGCCAAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATAGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCggg                                                                                                                                                                                                                                                   >  2:301190/1‑249 (MQ=255)
   cccGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTATTAACGTTGTAGTAGCTAAATAGCGGTATgc                                                                                                                                                                                                                                               >  2:38687/1‑251 (MQ=255)
                 cAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAgcg                                                                                                                                                                                                                                 <  1:26325/251‑1 (MQ=255)
                 cAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAgcg                                                                                                                                                                                                                                 <  2:184657/251‑1 (MQ=255)
                     aaCATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCg                                                                                                                                                                                                                             >  1:356158/1‑251 (MQ=255)
                          aGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGtctc                                                                                                                                                                                                                         <  1:317952/250‑1 (MQ=255)
                             aTGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCaaa                                                                                                                                                                                                                      >  2:270671/1‑250 (MQ=255)
                               gACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTAAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGCCCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGACTTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAAcc                                                                                                                                                                                                                   >  1:56905/1‑251 (MQ=255)
                                 ccAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGATCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCCCGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACcaa                                                                                                                                                                                                                 <  1:85390/251‑1 (MQ=255)
                                            gcgccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTgg                                                                                                                                                                                                      <  1:260071/251‑1 (MQ=255)
                                               ccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAg                                                                                                                                                                                                   <  1:15184/251‑1 (MQ=255)
                                                   cgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCg                                                                                                                                                                                               <  2:341994/251‑1 (MQ=255)
                                                     gTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGAGCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTAGGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGt                                                                                                                                                                                             <  2:153875/251‑1 (MQ=255)
                                                         ccagtcccgCGTACCGTTAGCAGTATGGAACGCCAGATCGGATAGCGGATTGTTTTCACCGTTGCTGACATAACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGACTTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCACCAAAGTTTTGGTAGCCCGGTGc                                                                                                                                                                                           >  1:149612/1‑249 (MQ=255)
                                                         ccagtccagCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa                                                                                                                                                                                          <  2:372114/250‑1 (MQ=255)
                                                           agtccagCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATTGTTGTAGGGGTAGTTGATCGACAGTTTCCGTTCATCGTTGAACGGAATTTTCAGTTCGGATTCCACGCCCTGCATACGAGCTTTGTTAAAGTTGTAGTAGCTAAATAGCGGTATGCGCCGGCCGCTAGCGCACGTCTCAAAACCAACCAAGTTTTGGTAGCCCGGTGCAGCg                                                                                                                                                                                       >  1:281624/1‑251 (MQ=255)
                                                                    tACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCg                                                                                                                                                                                               >  2:30481/1‑234 (MQ=255)
                                                                    tACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGGCCGTTATCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCg                                                                                                                                                                                               <  1:30481/234‑1 (MQ=255)
                                                                      ccGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTc                                                                                                                                                                               >  1:18711/1‑248 (MQ=255)
                                                                        gTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACAAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACg                                                                                                                                                                          >  2:248741/1‑251 (MQ=255)
                                                                           aGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGTAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTAc                                                                                                                                                                       <  1:101022/251‑1 (MQ=255)
                                                                                     aaCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGgctgatgct                                                                                                                                                             <  1:167411/251‑1 (MQ=255)
                                                                                       cGGCAGCTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACCGTTTCCATTCATCGTTGAACGGACTTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGgctgatgctg                                                                                                                                                            >  2:289894/1‑250 (MQ=255)
                                                                                          cAGATCGGATATCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATATCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATAc                                                                                                                                                        <  2:402504/251‑1 (MQ=255)
                                                                                            gATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGa                                                                                                                                                      <  1:51480/251‑1 (MQ=255)
                                                                                               cGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCt                                                                                                                                                   >  1:19611/1‑251 (MQ=255)
                                                                                                  aTAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGCATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTc                                                                                                                                                >  2:235409/1‑251 (MQ=255)
                                                                                                        gtttgtttTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCg                                                                                                                                          >  2:38673/1‑251 (MQ=255)
                                                                                                           tgtttTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGtt                                                                                                                                        <  1:444292/250‑1 (MQ=255)
                                                                                                           tgtttTCGCCGTTGCTGACAGCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAAGAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTTCAGCGTTGACGTCAGGCGTACGGCTGATGCTGATACGATCTTTCACATCGtt                                                                                                                                        <  2:27165/250‑1 (MQ=255)
                                                                                                                      ttGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGt                                                                                                                            >  1:158849/1‑251 (MQ=255)
                                                                                                                             cagcACGACCATCGTTGTAGGTGTAGTGGAGCGACAGTTTCCATGCATCGTTGAACGGAATTTTCAGTTCGGCTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGGGGTATGCGTCGTCCGTGAGCGCCCTTCTCAAAACCAACAAAGTTGTGGTAGCCCGGTGCAGCGTTGACGGCAGACGTTCGGCTGATGCTGATACGATCGTGCACATCGTTACGGAAAACGGTATCGCTg                                                                                                                     <  2:401146/248‑1 (MQ=255)
                                                                                                                                                 gtgtAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTccagcca                                                                                                 <  1:38673/251‑1 (MQ=255)
                                                                                                                                                     aGTTGATCGGCAGTTTCCAGTCAGCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTGAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCct                                                                                             <  2:18711/251‑1 (MQ=255)
                                                                                                                                                            cGACAGTTTCCATTCATCGTGGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACAGCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCa                                                                                       <  1:153293/250‑1 (MQ=255)
                                                                                                                                                             gACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCAcc                                                                                     >  1:328102/1‑251 (MQ=255)
                                                                                                                                                                aGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCAt                                                                                  >  1:60807/1‑251 (MQ=255)
                                                                                                                                                                             cGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGcccc                                                                     >  2:193956/1‑251 (MQ=255)
                                                                                                                                                                             cGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTGTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGTCTGATGCTGATACGAGCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGcccc                                                                     <  1:390647/251‑1 (MQ=255)
                                                                                                                                                                                tGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAg                                                                   <  1:95998/250‑1 (MQ=255)
                                                                                                                                                                                tGAACGGAATTTTCAGTTCGGTTTCCACGCCCTAAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCATATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGc                                                                  >  2:418712/1‑251 (MQ=255)
                                                                                                                                                                                                   ggTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGctggttt                                               <  2:191223/251‑1 (MQ=255)
                                                                                                                                                                                                              ccTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCtggtttctggtttc                                       >  2:217348/1‑248 (MQ=255)
                                                                                                                                                                                                                tGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCa                                                             >  2:122633/1‑224 (MQ=255)
                                                                                                                                                                                                                tGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAACTCCCa                                                             <  1:122633/224‑1 (MQ=255)
                                                                                                                                                                                                                 gAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGTTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTTATGCTGATTCGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATcc                                 <  1:315308/251‑1 (MQ=255)
                                                                                                                                                                                                                 gAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATcc                                 >  1:261395/1‑251 (MQ=255)
                                                                                                                                                                                                                        aGCTTTGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTAc                          <  1:270671/251‑1 (MQ=255)
                                                                                                                                                                                                                           tttGTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAAAGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTAccc                        >  2:165489/1‑250 (MQ=255)
                                                                                                                                                                                                                              gTTAACGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACaa                    <  1:217348/251‑1 (MQ=255)
                                                                                                                                                                                                                                tAACGTTCTCGTAGCTAAATAGCGGTATGAGCCGTACGTTAGCGTTCGTCTAAAAACCAACAAAGATTTGGTGGCCCGGTACAGCGTTGACGTCAGACGTACGGCTGATGCTGCTACGATCTTTCACATCGTTACGGAACACGGTCACGCTGGATTCAACCACTTCCAGCCAGCCTTCTTCACCCATGAAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCAGTGCTACCcaccctc                  >  2:222494/1‑247 (MQ=255)
                                                                                                                                                                                                                                  aCGTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACAATCtt                <  2:356158/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGCTACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACAATCTTAc              >  1:290639/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTTGTAGTAGCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACAATCTTAc              <  2:118688/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tagCTAAATAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGATTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACAATCTTACATGCGcc       >  1:108595/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aaaTAGCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACCATCTTACATGCGCCACGGc  >  2:352877/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCTGGTTTCAGATCCGGGCTACCCACAATCTTACATGCGCCACGGC  >  NC_000913/2244976‑2245466

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