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,2380TG95.0% 132.1 / ‑4.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 (20/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 = 1.00e+00

CTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGAT  >  NC_000913/4509993‑4510477
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
cTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGCGATTAAGAAACCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTg                                                                                                                                                                                                                                            >  2:150690/1‑251 (MQ=255)
    gcgcTTTCCCTCTGGCGAAAACGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTACGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGCTTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCttt                                                                                                                                                                                                                                        >  2:250965/1‑251 (MQ=255)
         ttCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGAt                                                                                                                                                                                                                                    <  1:30909/250‑1 (MQ=255)
                 ggCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTg                                                                                                                                                                                                                            <  1:322089/250‑1 (MQ=255)
                   cGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCg                                                                                                                                                                                                                         >  1:378430/1‑251 (MQ=255)
                           ggAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCa                                                                                                                                                                                                                 <  2:64013/251‑1 (MQ=255)
                            gagaAGGACGGGCCAGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGCAGATATCATTGATAACAATAGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAACCGAACACTTTTTAGAAAACTCGCTGGACAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCAt                                                                                                                                                                                                                >  2:240189/1‑251 (MQ=255)
                                aGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                                                                                                                                                                                            <  2:225654/251‑1 (MQ=255)
                                   acgggcaagTCTGCCTTCACTTTATGGCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGAAAACCATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCCAATACCTAAGTCGTAGCAAACGCACACTTTTTAAAAAAATCGATGGTCAGAACCAGGATTAGCAATTCCGTTCCAGGGTTGATTTTGCTATACGTTGGCGACTGATCATTGGttt                                                                                                                                                                                                          >  2:100692/9‑250 (MQ=255)
                                          cGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATAATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                                                                                                                                                                                   >  2:201827/1‑250 (MQ=255)
                                                           tGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTgg                                                                                                                                                                                 >  1:91034/1‑251 (MQ=255)
                                                                  ttCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGa                                                                                                                                                                           >  1:146157/1‑250 (MQ=255)
                                                                    cGGTAAGCTGGGTGAGACACCGACTCATGAAGTTATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAGCGGGATTAGCAATTCCTGTCCATGGTTGCTTTTGATTTACGTTGGCGTATGATCATTGATTTATCCTCAAAAGCCCATCGTCATTAGTAATGAAcc                                                                                                                                                                        <  2:139315/251‑1 (MQ=255)
                                                                    cGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAAc                                                                                                                                                                         <  1:123072/250‑1 (MQ=255)
                                                                                           cTCATGAAGATATCATTGCTAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTa                                                                                                                                                                                                                                 >  2:128537/1‑171 (MQ=255)
                                                                                           cTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTa                                                                                                                                                                                                                                 <  1:128537/171‑1 (MQ=255)
                                                                                             cATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGATGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGATACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGc                                                                                                                                                >  1:17704/1‑250 (MQ=255)
                                                                                                gAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACag                                                                                                                                            <  2:222193/251‑1 (MQ=255)
                                                                                                gAAGATATCATTGATAACAAGCGATCCAGCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCATACTTTTTAGAAAAATCGATGGTCAGAATCTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGAGTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACag                                                                                                                                            <  1:270418/251‑1 (MQ=255)
                                                                                                  aGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAATGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACag                                                                                                                                            >  1:225068/1‑249 (MQ=255)
                                                                                                            gATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAAtttt                                                                                                                                 <  2:365445/250‑1 (MQ=255)
                                                                                                                   aTCGATCCATCAACTGGCCTGTATAACGTGTCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTTAtctc                                                                                                                         <  1:90128/251‑1 (MQ=255)
                                                                                                                           aTCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGa                                                                                                                 <  1:68592/251‑1 (MQ=255)
                                                                                                                            tCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGACACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGACATTTTTTGATCTCCCCCCTGa                                                                                                                 >  2:127738/1‑250 (MQ=255)
                                                                                                                                  ggCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCGCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCt                                                                                                            <  1:241211/249‑1 (MQ=255)
                                                                                                                                           aaCGTGGCTGCCGGTGATTAAGAAAGCTGCACATACCTAATGAGTAGCAAACGCACCCTTTTTAGAAAAAGCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCGCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAAGGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAggca                                                                                                 <  2:378430/251‑4 (MQ=255)
                                                                                                                                                 gCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGCAGTAGCAAACGCAACCTTATTAGAAAAATCCATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGCGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGACAAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGtaac                                                                                           >  1:392059/1‑248 (MQ=255)
                                                                                                                                                               aGAAAACTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTACGTTCCAGGGTTGATTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCt                                                                             >  2:191623/1‑251 (MQ=255)
                                                                                                                                                                  aaGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTAATGTCACGTTTTTGGGTTTCTaa                                                                           >  1:345060/1‑250 (MQ=255)
                                                                                                                                                                       gCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGTGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAATCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAttt                                                                      <  1:71659/250‑1 (MQ=255)
                                                                                                                                                                         acctacctAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTa                                                                    >  1:66705/1‑250 (MQ=255)
                                                                                                                                                                            tacctaAGTAGTAGCAAACGCACACTTTTTAGAACAATCGCTGGTCAGAAACTGGATTAGCAATTCCGATCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACtt                                                                >  1:416294/1‑251 (MQ=255)
                                                                                                                                                                                taAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTACTGAACCAGCTCCGt                                                                                                                                                                >  1:300827/1‑151 (MQ=255)
                                                                                                                                                                                taAGTAGTAGCAAACGCACACTGTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCGCAAAAGCCCATCCTCAGTGTTAATGAACCATCTCCGt                                                                                                                                                                <  2:300827/151‑1 (MQ=255)
                                                                                                                                                                                                  tacTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAATGCCCAACCTCATTGGTATTGAACCAGCTACGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGCTCTCCCCCCTGAAATCTAAACTTAGTAATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCAGTAAgg                                          >  2:229139/2‑251 (MQ=255)
                                                                                                                                                                                                           gAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAg                                 <  1:250965/251‑1 (MQ=255)
                                                                                                                                                                                                                  cGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTc                          <  1:157208/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTCAGGGTTATACAGAAATGTCCg                        >  2:115338/1‑251 (MQ=255)
                                                                                                                                                                                                                          aGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTAATTAGTAATGATCCAGCTCCGTGAATGTCCGCGCTGGCACAGAGCGAAATTTTTTGATCCACTCCCTGAAATCTAAACTTAGTAATGTCGCGTTGTTGGGTATCTAAAAATTTTGCTTCGCGCTTGTCGTTGCCGTAAGGGTTTTCCAGAACTGTCCGTTaggc                  <  2:72331/251‑4 (MQ=255)
                                                                                                                                                                                                                                          aTTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAattgat  >  1:225067/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
CTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGAT  >  NC_000913/4509993‑4510477

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: