Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,642 T→C V256V (GTA→GTG fecC ← ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,6420TC96.0% 220.2 / ‑3.2 76V256V (GTA→GTGfecCferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (34/39);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATG  >  NC_000913/4512406‑4512881
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
aaTTTTCATCCTCGCCTCCTCACAAGCCAGACAAGGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGCTCTCCGGGGAAGTCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGCCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCGGCCCAGAAGCGCGCCCGATGTGGCACCAGCAGACCGCTCAACCCACCCGGACCTGCCACGCtaacccacac                                                                                                                                                                                                                                   >  2:78870/1‑242 (MQ=255)
   tttCATCCGCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACgcacgcac                                                                                                                                                                                                                                 <  1:83857/250‑1 (MQ=255)
        cccccccctcccCACACGCCAGAAAAAGCAAGGGCGGCAATACAGCGCCAGCCCTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCaccaacca                                                                                                                                                                                                                           <  2:406016/245‑1 (MQ=255)
         ccTCGCCTCCTCCCAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCGGCGGGCAGATCTCCGGGAAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAg                                                                                                                                                                                                                          <  2:28929/251‑1 (MQ=255)
            cGCCTCCGCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCAGGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAACCGCCACCGGGCCTGCCACGCTGACGCACGCACCAACCAGAag                                                                                                                                                                                                                       <  1:20986/251‑1 (MQ=255)
              cctcctCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAagc                                                                                                                                                                                                                      >  1:127626/1‑250 (MQ=255)
                 cctcACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCCCGGCGCCGGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGCCCTGCCACGCTGACGCACGCACCAACCAGAAgcagca                                                                                                                                                                                                                  <  2:483227/251‑1 (MQ=255)
                 cccaacaAGCAAGAAAAAGCAAGGGCTGCCAATCAGGGCCAGCACTGCGCCTGCGGGCAGACCTCCGGGGAAGGCCAGTGCGCGGGCGAGTACATCTGCCAGCAGCATAAGCGTGGCCCCCAGCAGCATGCTCACTGGCCGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGAAGACCGATAAACGCCACGGGACCAGCCACGCTGACGCACGCACCAACCAGAAgcagca                                                                                                                                                                                                                  <  2:224703/247‑1 (MQ=255)
                                aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGACGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGCTAAACGCCCCCGGACCTGCCCCGCTGACGCACGCACCACCCAGAAGCCACACTAACCTATtgatc                                                                                                                                                                                                    >  2:400198/1‑249 (MQ=255)
                                aaaGCAAGGGCTGCCAACAAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGGGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgat                                                                                                                                                                                                     <  1:400198/249‑1 (MQ=255)
                                 aaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGAc                                                                                                                                                                                                  <  1:288981/251‑1 (MQ=255)
                                      aGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTTATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAAc                                                                                                                                                                                             <  1:167653/251‑1 (MQ=255)
                                        gccggCCAATCACCGCCAGCACGCAGCCTGCGGGCGATGCCCCGGGGAAGGCCAGTGCGCGTGCGAGCACATCTGCCACCAGCAGCGGCTTGGCCCCCAGCAGAATGCCCACTGGCAGTAGTTTGCGCTGTCCGAAGCCTGCCCAGAAGCGCGCCAGAGGTGGCACCAGCAGACCGATAAACGCCACCGGAACTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGt                                                                                                                                                                                           <  2:305247/247‑1 (MQ=255)
                                        gccggCCAACCACCGCCAGCACTGCGCCGGCGGGCAGACTGCCGGGGAAGGCCAGTGCGCGGGCGAGTCCAGCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGt                                                                                                                                                                                           <  2:417939/247‑1 (MQ=255)
                                             ccAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGccc                                                                                                                                                                                      >  2:417566/1‑250 (MQ=255)
                                               aaTCAGCGCCAGCACTGGCCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGGGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCAGGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCTATAAACGCCTCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCg                                                                                                                                                                                    <  2:480266/251‑1 (MQ=255)
                                                aTCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGCCCTGCCACGCTGACGCCCGCCACAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGt                                                                                                                                                                                   >  2:152084/1‑251 (MQ=255)
                                                   cagcgccagcgcTGCGCCTGCGGGCAGATCTCCGGGGAAGACCAGTGCGCGTGCGAGTACATCTGCCAACAGCATCAGATTGGCCCCCAGCAGCATGCTCACTGGCAGTCCGTTGCGCTGATCGAAGCCTGCCAAGAAGCGCGCAAGATGTGGCACCAGCCGACCGATAAACGCCACCGGACCTCCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACACATTGCTGACCAGACGCACCCTCgcc                                                                                                                                                                                  >  2:352871/1‑248 (MQ=255)
                                                         catcACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCACAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGTACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCa                                                                                                                                                                           <  1:170926/248‑1 (MQ=255)
                                                            cACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTc                                                                                                                                                                        >  2:91177/1‑251 (MQ=255)
                                                             aCTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACCTCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCCCGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                                                                                                       >  2:60617/1‑251 (MQ=255)
                                                                gcgcCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTccc                                                                                                                                                                      <  1:349777/249‑1 (MQ=255)
                                                                       gCGGGCAGACCTCCGGGGAAGGCCAGTGCGCGTGCGAGGACCTCTGCCCGAAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCTGTCCGTTGCGCTGATCGAAGCCTGCCCAGAGGCGCGCCAGATGTGGCTCCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCCCTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTa                                                                                                                                                                <  2:275679/249‑1 (MQ=255)
                                                                        ggtgccagtcctcCGGGGCAGGCCCGTGCGCGTGCGAGTACATCTGCCAACAGCCTCAGCGTGCCCCCCAGCAGCTTGCTCACTGGCAGTCCGTTGCGCTGACCGAAGCCTGCCCAGAAGCGCGCCAGAGGTGGCACCAGCAGACCGATAACCGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATg                                                                                                                                                              <  2:386581/241‑1 (MQ=255)
                                                                                     ccGGGGAGGGCCCGTGCGCGTGCGAGTACATCCGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCACAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTg                                                                                                                                                    <  2:268333/249‑1 (MQ=255)
                                                                                     ccGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc                                                                                                                                                  <  1:505107/251‑1 (MQ=255)
                                                                                     ccGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc                                                                                                                                                  >  1:465644/1‑251 (MQ=255)
                                                                                      cGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCCTGCTCCCTGGCAGTACGTTGCACTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACTCCACCGGCCCTGCCCCGCTGACGCCCGCACCACCCAGAAGCAGCACTAACATATTGAGGACACAACGAAACCTCGTCAGGTTCCCTCCCCGCGTATGGGAGGTgaggtcc                                                                                                                                                 >  2:164414/1‑245 (MQ=255)
                                                                                       ggggAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGc                                                                                                                                                >  1:366981/1‑251 (MQ=255)
                                                                                                 agtgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTg                                                                                                                                        >  1:319440/4‑249 (MQ=255)
                                                                                                 agtgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTg                                                                                                                                        >  1:252803/4‑249 (MQ=255)
                                                                                                   tgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGGAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGc                                                                                                                                     >  1:370654/2‑250 (MQ=255)
                                                                                                   tgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGAACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGCTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGGGGTTGAGcg                                                                                                                                    >  2:311214/2‑250 (MQ=255)
                                                                                                   ggcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATTCTCACTGGCAGTACGTTGCGCTGATCGAAGCCGGCCCATAAGCGCGCCAGATGTGGCGCCAGCAGACCGATAAACGCCCCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCATCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCa                                                                                                                                    <  1:493240/250‑1 (MQ=255)
                                                                                                           cGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGCGGTTGAGCAGGTTCAGt                                                                                                                            >  1:110834/1‑251 (MQ=255)
                                                                                                              gTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGa                                                                                                                         <  1:123359/251‑1 (MQ=255)
                                                                                                                aCATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGAtt                                                                                                                       <  1:60617/251‑1 (MQ=255)
                                                                                                                    ccgacgtccgcagcggcgggccccccTGCAGAATGCTCTCGGGCAGTCCGTGGCTCTGATGAAACCCTCCCCAGAGGCGCGCTATATGTGGCACCAGAAGACAGATAAACGCCACCGTACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGcc                                                                                                                   <  2:5984/230‑1 (MQ=255)
                                                                                                                      gCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGcc                                                                                                                   >  1:498380/1‑249 (MQ=255)
                                                                                                                        cagcagCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGCTGACCAAACGTAGCCTCGTCCGGTTCACCCCCAGCGTATGGGCGGTGCTGTCGCTTGGGTTGAGCAGGTTCAGTTGATTCGCcagc                                                                                                                >  2:18250/1‑250 (MQ=255)
                                                                                                                         agcagcATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGATGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGGTCAGTTGAATCGCcagc                                                                                                                >  2:243341/1‑249 (MQ=255)
                                                                                                                              caTCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAac                                                                                                          >  1:300725/1‑250 (MQ=255)
                                                                                                                                tCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAAACAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACg                                                                                                       >  1:65956/1‑251 (MQ=255)
                                                                                                                                 cAGCGTGGCCCCCAGCAGCAGGCGCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGa                                                                                                      <  2:319440/251‑1 (MQ=255)
                                                                                                                                      tGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGCCCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGAGCGGTCGCTGAGGTTGGGCAGGTTCAGATGCTTCGCCAGCAGCAACACGCCAggg                                                                                                 >  2:335122/1‑251 (MQ=255)
                                                                                                                                       ggCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAgggg                                                                                                >  2:235123/1‑250 (MQ=255)
                                                                                                                                        gCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGCTAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCAATCCCAGCGTAAGGGGGGTGCTGTCGTTGCGGTTGAGCAGGTTCAGTTGCTTCGCCAGCAGCAacccccccgggcc                                                                                               >  2:472665/1‑240 (MQ=255)
                                                                                                                                        gCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGAc                                                                                               <  1:113622/251‑1 (MQ=255)
                                                                                                                                         cccccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAACGTATGGGCGGTGCTGTCGCTGAGGCTGTGCAGGGTCAGTTGATTCGCCAGCAGCAACACGACAggggcg                                                                                              >  1:138489/1‑248 (MQ=255)
                                                                                                                                          ccccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGCTAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGCGGTTGAGCAGGGTCCGTTGATTCGCCAGCCGCAACCCGACAGGGAcgg                                                                                             >  2:186218/1‑249 (MQ=255)
                                                                                                                                              agcagcATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGt                                                                                         <  1:417566/251‑1 (MQ=255)
                                                                                                                                                cagcaTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGGCTGCAGTg                                                                                        >  1:455337/1‑250 (MQ=255)
                                                                                                                                                  gcaTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGacc                                                                                     <  1:152084/251‑1 (MQ=255)
                                                                                                                                                         cacgggcagtactgtgcgcggagagAAGCCTGGCCGGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGCCCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGaccaccac                                                                                <  2:310013/225‑1 (MQ=255)
                                                                                                                                                            tGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAACCGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAagca                                                                                                                                                                                                                     >  1:221519/1‑113 (MQ=255)
                                                                                                                                                            tGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCAAGATGTGGCACCAGCAGACCGATAAACGCCCCCGGACCTGCCACGCTGACGCACGCACCAACCAGAagca                                                                                                                                                                                                                     <  2:221519/113‑1 (MQ=255)
                                                                                                                                                               cAGTACGTTGCGCTGATAGAAGCCTGCCAAGAAGCGCGCCAGATGTGGCCCCAGCAGACCGCTAACCGCCACCGGACCTGCCACGCTGACGCACGCACCAACAAGACACAGCACTAACATATTGATGACCAAACGTAGCCTCGTAAGGTTCACTCCCAGCGTATGGGCGGGGCAGTCGCTGAGGATGAGCAGGTTCAGTTGATTCGCACGCAGCAACACGACAGGGACTGCAGttgccacaacaggcaaa                                                                         >  2:244353/1‑234 (MQ=255)
                                                                                                                                                                            tGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATc                                                           >  1:100059/1‑251 (MQ=255)
                                                                                                                                                                             gATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCGCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACAt                                                            <  1:260891/249‑1 (MQ=255)
                                                                                                                                                                               ccGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTg                                                        <  1:189931/250‑1 (MQ=255)
                                                                                                                                                                                   aGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGcc                                                      >  1:141990/1‑249 (MQ=255)
                                                                                                                                                                                    gccgccccagagacgcgccgaaggtgGCAACAGCAGACCTATAACCGCCACCGCACCAGCCACGCTGACGCACGCACCAACCAGAAGCAGCCCTAACATATTGAGGACCAAACGTAGCCTCGTTAGGTTCACTCCCAGCTTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGc                                                   <  2:344933/228‑1 (MQ=255)
                                                                                                                                                                                       tGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCggg                                                >  1:412878/1‑251 (MQ=255)
                                                                                                                                                                                           cAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTg                                            <  1:235123/251‑1 (MQ=255)
                                                                                                                                                                                                   cgcCAGATGTGGCCCCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTc                                    <  1:360095/251‑1 (MQ=255)
                                                                                                                                                                                                      cAGATGTGGCTCCAGCAGACCGATAAACGCCACCGGACCGGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGc                                 <  1:52706/251‑1 (MQ=255)
                                                                                                                                                                                                             ggcacaaccagACCGAAAACAGCCCCCGGACCTGCCACGCGGACGCCCGCACCAACCAGAAGCAGCATGAACATATTGAGGACCAAACGTAGCCGCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCGGCCAGCGGGCGTGGGACACTCCGCCTgccag                          <  2:377677/243‑1 (MQ=255)
                                                                                                                                                                                                                 cacgcgaaCCGATAACCCCCCCCGGCCCTCCCACGCCGACGCCCGCCCCAACCAGAAGCAGCACTAACATATGGACGCCCAACCGGAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGccagccag                      <  2:397692/244‑1 (MQ=255)
                                                                                                                                                                                                                     cagACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGGGGGACACTCCGCCTGCCCGCCAGTaa                   >  2:102921/1‑250 (MQ=255)
                                                                                                                                                                                                                      agACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTaaaa                 <  1:18250/251‑1 (MQ=255)
                                                                                                                                                                                                                              aaaCGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGt         >  1:473332/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCAt   <  1:91177/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTCGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAgg                                                                                                                                           >  2:454072/1‑114 (MQ=255)
                                                                                                                                                                                                                                      aCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAgg                                                                                                                                           <  1:454072/114‑1 (MQ=255)
                                                                                                                                                                                                                                       cccgaaCCTGCCACGCTGACGCACGCCCCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATg  <  1:491137/246‑1 (MQ=255)
                                                                                                                                                                                                                                       ccGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGGTCAGTTGATTCGCCAGCAGCACCACGACCGGGACTGCAGTGACCACCACAGGCAAGAGCTGCCCGACATCCTGCCAGCCGGCGATGGACACGCCGCCTGCCAGACAGCAAAAGAtaacgaacacatg  >  2:104214/1‑238 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
AATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATG  >  NC_000913/4512406‑4512881

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