Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,510,238 T→G intergenic (+445/+452) insO → / ← fecE IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,510,2380TG94.7% 116.7 / ‑1.3 40intergenic (+445/+452)insO/fecEIS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (18/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

TCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGA  >  NC_000913/4509995‑4510476
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
tCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACGCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑c                                                                                                                                                                                                                                          <  1:263620/250‑1 (MQ=255)
          ccTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCCCTTTATGCCAGCTTCGGTAAGCTGGGTGCGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGACTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTATTAGAAAAATCGAGGGCCAGAAACGGGATTACCAATTCAGTCCCA‑GGTTGGCTTTTGATTTAc                                                                                                                                                                                                                               >  1:232774/1‑251 (MQ=255)
              gggCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTg                                                                                                                                                                                                                           <  2:320387/250‑1 (MQ=255)
               ggCGAAAAAGGGAGGAGGACGGGCCCGTCTGCCTTAAATTTATGCCAGCTGCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATACCAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCACACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTgg                                                                                                                                                                                                                          >  1:62396/1‑251 (MQ=255)
               ggCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACGTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGt                                                                                                                                                                                                                             <  2:254591/248‑1 (MQ=255)
                   aaacaGGGAAAAGGACAGGCCCGTCTGCCTTCACTTCATGCCAGCTTCGGAAAGCTAAGTGGGCCACCGACTCATGGCGACATCATTGATACCAATCGATCCATAAACTGGCCTGTAAAAGCTGGCTGACGGTGATTAAGAAAGATGCACCTAACTAAGTAGTAGCAAACGCCCACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTc                                                                                                                                                                                                                      >  2:367021/5‑251 (MQ=255)
                   aaaaaGGGAGAAGGACGGGCCCATCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTc                                                                                                                                                                                                                      <  1:45600/251‑1 (MQ=255)
                      aaGGGCGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTg                                                                                                                                                                                                                    >  2:21177/1‑250 (MQ=255)
                        gggAGGAGGACGGGCCCGTCGGGCTTCACTTTATGCCAGCTTCGGTAAGCTGGGGGAGACACCGACTCATGAAGATATCATTGATAACAAACGATCCATCAACTGGAATGTATAACGTGGATGACGGTGATTAAGAAAGATGCACCTACCTAAGTAGTAGCACACGCACACTTTTTAGACAAATAGATGGACAGACACGGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGAGTCTGATc                                                                                                                                                                                                                 >  2:79864/1‑251 (MQ=255)
                               agACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGGAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCACTTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATAATTGAtt                                                                                                                                                                                                          >  2:225873/2‑251 (MQ=255)
                                gACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                                                                                                                                                                                         >  2:226657/1‑251 (MQ=255)
                                   gggCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTGTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCTTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATc                                                                                                                                                                                                      <  1:356277/251‑1 (MQ=255)
                                      cccGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTc                                                                                                                                                                                                   <  1:10155/251‑1 (MQ=255)
                                                       tATGCCAGCTTCGGGAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGAGCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGAGTTATCCTCAAAAGCCCAACCTCAt                                                                                                                                                                                   <  2:112458/250‑1 (MQ=255)
                                                              gCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATg                                                                                                                                                                           >  1:210216/1‑251 (MQ=255)
                                                                         cTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATcc                                                                                                                                                                                                     <  2:67597/214‑1 (MQ=255)
                                                                         cTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATcc                                                                                                                                                                                                     >  1:67597/1‑214 (MQ=255)
                                                                           gggTGAGACACCGACTCATGAAGATATCAGTGATAACAAGCGATCCATCAACTGGCCTGTATAACGTGGCTGGCGGTGATTAAGAAAGCTGCACCTACCTAATTAGTAGCAAACGCACACTTTTTAGAAAAAGCGATGGTCAGAATCTGGATGAGCAATTCCGTTCCAGGGTTG‑CGTTTGATTTACGTTGGCGTCTTATCATTGATTTATCCTCAGAAGCCCAACCTCATTTGTAGTGAGGCAGCTCCGTg                                                                                                                                                              <  2:385138/251‑1 (MQ=255)
                                                                                       gACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCAATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGAGAATGTCCGcact                                                                                                                                                  >  1:305196/1‑248 (MQ=255)
                                                                                       gACTCATGAAGATATCAGTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCTCCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATGTACGTTGGCGTCGGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGctct                                                                                                                                                  <  2:24152/251‑1 (MQ=255)
                                                                                                  atatCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCg                                                                                                                                       >  1:128628/1‑251 (MQ=255)
                                                                                                    atCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGa                                                                                                                                      >  1:148923/1‑250 (MQ=255)
                                                                                                      cATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGCAAt                                                                                                                                   >  2:144842/1‑251 (MQ=255)
                                                                                                                aaTCGATCCATCCACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAACCGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAtct                                                                                                                         >  2:88217/1‑251 (MQ=255)
                                                                                                                       ccATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCt                                                                                                                  <  1:21177/251‑1 (MQ=255)
                                                                                                                       ccATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCt                                                                                                                  >  2:302685/1‑251 (MQ=255)
                                                                                                                        cATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTGCCTAAGTAGTAGCAAACGCGCACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAGTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTAGCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTg                                                                                                                 <  2:167951/251‑1 (MQ=255)
                                                                                                                         aTCAACTGGCCTGTATAACGTGGCTGCCGGTGATGAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGGTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGATCCAGCTCCGTGAATGTTCGCTCTGTCACAGAGCGAAATTTTTTGATCTCCCGCCTGa                                                                                                                <  1:141261/251‑1 (MQ=255)
                                                                                                                            aaCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAt                                                                                                             <  2:48905/251‑1 (MQ=255)
                                                                                                                             aCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTACCCCCTGAAATc                                                                                                            >  1:124303/1‑251 (MQ=255)
                                                                                                                               tGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCt                                                                                                           >  1:426827/1‑250 (MQ=255)
                                                                                                                                         aaCGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTc                                                                                                 <  1:226657/250‑1 (MQ=255)
                                                                                                                                                  gCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACg                                                                                         <  2:124303/249‑1 (MQ=255)
                                                                                                                                                                                     tagCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGCTCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCg                                                    >  2:159755/1‑251 (MQ=255)
                                                                                                                                                                                                             aaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAAAGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAAt                            >  1:174788/1‑251 (MQ=255)
                                                                                                                                                                                                                       cAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCATCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTa                    >  2:327479/1‑249 (MQ=255)
                                                                                                                                                                                                                              tGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTAAGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTCATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTACAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCTGAGtt           >  1:174669/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTGTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGt            <  1:155992/249‑1 (MQ=255)
                                                                                                                                                                                                                                 aTTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGACCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCaa        <  2:256034/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaTTCCGTTCCAGGGTTG‑CTTTTTATTTACGTTGTCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTTCGCGCTGGCACAGAGCGAAATGTTTTGATCTCCCCCCTGAAATCTAAACTTATTCATGTCTCGTTTTTGGGTTTCGAAAATTTTAACTGCGCGTTTTTCGTTGCCGTAAGGGTTATACAGTAATTTCCGTTAAGCTGAGTTCAAAattga  <  2:405402/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
TCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGA  >  NC_000913/4509995‑4510476

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