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,2380TG97.8% 151.2 / ‑4.8 46intergenic (+445/+452)insO/fecEIS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (19/26);  minor base . (0/1);  total (19/27)
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.91e-01

TCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCC  >  NC_000913/4509995‑4510481
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
tcacgcTTTCCCTCTGGCGAAAAAGGGAGAAGTACGAGCCCGTCTGCCTTCACTTTAGTCCCGCTGTCG‑TAAGCTGGGTGAAACCCCGACTCATGAAGATATAATTGATAACAAGCGATCCATCAACTGGCCTTTATAACGTGGCTGCCCGTGATTAAGAAAGCTGCACCTACCTAAGTAGGAGCAAACGCACACTTTTTAGAGAAATCGATGGTAAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑Tg                                                                                                                                                                                                                                                <  2:219008/246‑1 (MQ=255)
        tCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGAttt                                                                                                                                                                                                                                      >  1:51147/1‑251 (MQ=255)
          ccTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTAc                                                                                                                                                                                                                                    <  1:358983/251‑1 (MQ=255)
                gCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGc                                                                                                                                                                                                                              <  1:185871/251‑1 (MQ=255)
                         ggAGAAGGACGGTCCCGTCTGCCTTCACTTTATGCCAGCA‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCATCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCCCCCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTTGCGTCTGATCa                                                                                                                                                                                                                     <  2:513561/251‑1 (MQ=255)
                         ggAGAAGGACGGGCCCGTCTTCCTTCACTTTATGCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCa                                                                                                                                                                                                                     <  2:77982/251‑1 (MQ=255)
                                         gTCTGCCTTCACTTTATGCCAGAT‑TCGGTAAACTGCGTGAGACACCGACTCATGCAGATATCATTGATAACAATCGATCCATCAACTGGCATGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCCCACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCCATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCaa                                                                                                                                                                                                      >  2:361971/1‑250 (MQ=255)
                                             gCCTTCACTTTATGCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAg                                                                                                                                                                                                   <  1:40554/249‑1 (MQ=255)
                                                   aatttATACCGGCT‑TCGGTAAGCTGGGTGTGACACCGACTCATCAAGATATCATTAATAACAATCGATCCATCAACTGGCCTGATTAACGTGGCTTCCAGAACATAAGAAAGCTGCACAGACATAAGTAGTAGCAAACGCACACTTTTTATAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCAAGGAT‑AGCATTTACTTTACGTTGGCGACTTATCATTGGTTTATCCTCAAAAGccctatct                                                                                                                                                                                           >  2:164155/3‑246 (MQ=255)
                                                          gCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGt                                                                                                                                                                                    >  1:130149/1‑251 (MQ=255)
                                                              gCT‑TCGGTAAGCTGGGTGAGACACCGACACATGAAGATATCATTGATATCAATCGATCCATCAACTGGCCTGTATTACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑GGCTTTTTATGTACGTTGGCGTCGGATCATTGATTTATTTTCAAAATCTCATCCCCATGGGTAATg                                                                                                                                                                                <  2:75934/251‑1 (MQ=255)
                                                                   cGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAAc                                                                                                                                                                             <  1:186595/250‑1 (MQ=255)
                                                                    ggTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGCTTTATCCTCAAAAGCCCAACCTCATTGGTCATGAAcc                                                                                                                                                                            >  1:452491/1‑250 (MQ=255)
                                                                                gagaCACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAAt                                                                                                                                                                <  1:56481/250‑1 (MQ=255)
                                                                                gagaCACCGACTCATGAAGATAGCACTGATAACAATCGATCCATCAACTGGCCGGTATAACGTGGCTGCCGGTGAGTAAGAAAGCTGCACCGTCCGAAGTAGTAGCAAACGCACACTTTGTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGTT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATg                                                                                                                                                               <  2:472704/251‑1 (MQ=255)
                                                                                       cGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTTGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCATCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGctc                                                                                                                                                        >  2:334177/1‑251 (MQ=255)
                                                                                         aCTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTg                                                                                                                                                      >  2:23694/1‑251 (MQ=255)
                                                                                                       cATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCCCCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAAt                                                                                                                                        <  1:60103/251‑1 (MQ=255)
                                                                                                       cATTAATAGCAAACGATACAGCAACTGGCCTGTATAACGTGGATGCCGGTGATTAAGAAAGCTGCACCTAACTAAGTAATAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAGATGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATAATTGATTTATCCTCAAACGCCCAACCTCATTGGTAATGAACCAGCTCCGTGCATGTCCGCGCTGGCGCGGAGCGaact                                                                                                                                        >  2:86581/1‑249 (MQ=255)
                                                                                                        aTTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGACTGTCCGCTCTGGCACAGAGCGAAAtt                                                                                                                                       >  1:344422/1‑251 (MQ=255)
                                                                                                          tGATAACAATCGATCCATCAACTGGCCTGTATCACGTGGCTGCAGGTGATTACGAAAGATGCACCTACCTAAGTAGTAGCAACCGCACACTTTTTAAAAAAATCGATGGTCAGACACTGGATTAGCACTTCCGTTCCAGGGT‑TGCTTTTGATTTTCGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCACACTCTTTGTTAATGAACCAGCTCCGTGAATGTCCGCTCTGACACAGAGCGAAAtttt                                                                                                                                     >  2:185315/1‑251 (MQ=255)
                                                                                                                    cGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAACGATGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGATGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACATCATTTGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCAAAGAGCGCACTTTTTTGATCTccc                                                                                                                           >  2:43202/1‑251 (MQ=255)
                                                                                                                                  gCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCaa                                                                                                                                                                                              >  1:109171/1‑170 (MQ=255)
                                                                                                                                  gCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCaa                                                                                                                                                                                              <  2:109171/170‑1 (MQ=255)
                                                                                                                                                   gCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTATAAAAATCGATGGTCAGAAACTTGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGtt                                                                                            <  2:12368/251‑1 (MQ=255)
                                                                                                                                                   gCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGtt                                                                                            <  1:23694/251‑1 (MQ=255)
                                                                                                                                                            tAAGAAAGCTGCACCTACCGAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCATCCTCATTGGTAGTGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGtt                                                                                    <  1:529342/250‑1 (MQ=255)
                                                                                                                                                               gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTa                                                                                <  1:475782/251‑1 (MQ=255)
                                                                                                                                                               gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTGTCTa                                                                                <  1:143651/251‑1 (MQ=255)
                                                                                                                                                                  aGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTaaag                                                                             <  1:376314/251‑2 (MQ=255)
                                                                                                                                                                   gCTGCACCTACCTAAGTAGTAGCAAACGCACACTGTTTAGAAAAATCGATGGTCAGAAACTTGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCGCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAt                                                                            <  2:167866/251‑1 (MQ=255)
                                                                                                                                                                      gccactccctaAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCCTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAGCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAtttt                                                                         <  1:145856/244‑1 (MQ=255)
                                                                                                                                                                                 agtagtATCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTGCGCGtt                                                              <  2:215222/251‑1 (MQ=255)
                                                                                                                                                                                     gtagCAAACGCTCACTTTTTAGAAAAATCGATGGTCAGAAACGGGATTAGCAATTCCGTTCCCGGGT‑TGATTTTGATTTACGTTTGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTc                                                          <  2:130149/251‑1 (MQ=255)
                                                                                                                                                                                     gtagCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTTCGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTc                                                          <  1:185315/251‑1 (MQ=255)
                                                                                                                                                                                           aaCGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGcc                                                    >  1:345921/1‑251 (MQ=255)
                                                                                                                                                                                                  acTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCATCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAggg                                             <  1:334177/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACCGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTc                              >  2:14589/1‑251 (MQ=255)
                                                                                                                                                                                                                       tCAGAAACTGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCTTGACTGTCCGCTCTGGCACATAGCGAAACTTTTTAATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACATTTTTGTGTTTATAAAATTTTCAATTCGCGATTTTAGTTGCCGTAAGGGTTATACATAAATTTCAGTTa                         >  2:424526/1‑250 (MQ=255)
                                                                                                                                                                                                                             aatGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCCTGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTCTACAGAAATGTCCGTTAAGCag                    >  2:111066/3‑249 (MQ=255)
                                                                                                                                                                                                                               tGGATTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTACGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCCCGTTTTTGGGTTTCTAAAATTTTAACTTAGCGTTTTTCTTTGCCTTATGGGGTATACAAAAATGTACGTTAAGCAGAGtt                >  2:451746/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aTTAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCa              >  1:541649/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tAGCAATTCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGGCAGGTCACGTTTTTg                                                                                        <  1:127626/174‑1 (MQ=255)
                                                                                                                                                                                                                                           tCCGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGCAATGTCCGTTAAGCAGAGTTCAAAattgattg    >  1:103702/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tCCGTTCCA‑GGTATGCTTTTGAATTGCGTTGGCGGCTTATCATTGATTTAACCTCAAAAGCACAACCTCATTGGTAATGAATCAGCTCCGTGAATGTCCGCTTTGGCACAGAGCGAAATGTTTTGATCTCCCCCCTTAAATCGAAACTTAGTCATGTCACGTTTTTGGGTTTATAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTTTACAGAAATGTCCGTTAAGCAGAGTTCAAAattgatt     <  2:38547/250‑1 (MQ=255)
                                                                                                                                                                                                                                             cGTTCCAGGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGcc  <  1:47878/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
TCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCT‑TCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGT‑TGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCC  >  NC_000913/4509995‑4510481

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