Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,433,665 C→T intergenic (‑105/‑25) smf ← / → def protein Smf/peptide deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,433,6650CT96.4% 168.0 / NA 56intergenic (‑105/‑25)smf/defprotein Smf/peptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/26);  minor base G (0/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTA  >  NC_000913/3433420‑3433908
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
ccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTg                                                                                                                                                                                                                                                >  1:26548/1‑251 (MQ=255)
ccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTg                                                                                                                                                                                                                                                >  1:64300/1‑251 (MQ=255)
   cAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAac                                                                                                                                                                                                                                             >  1:100676/1‑251 (MQ=255)
    aaTGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAaca                                                                                                                                                                                                                                            <  1:553323/251‑1 (MQ=255)
       gTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACAt                                                                                                                                                                                                                                         >  1:323934/1‑251 (MQ=255)
        tAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACAtc                                                                                                                                                                                                                                        >  2:10076/1‑251 (MQ=255)
              cTGTTTGCTGAAATACAACCGCATCAATATGCGACTGTTTTGCCACCCCGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTCTCGACAATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTACTTCCTGCACACATCTCTgg                                                                                                                                                                                                                                   >  2:226799/1‑250 (MQ=255)
                gTTTGCTGAAATACAACCGCACCAATATGCGATTGTTTTGCCACCCAGTGAGCGATACGGACCATATAATCGCCGTACAAGCTGCTGATACTCATTAAACGAAGCCAAATTTCTGTACCGACCTTCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATGCTTGCACACATCTCTGga                                                                                                                                                                                                                                  <  1:200120/249‑1 (MQ=255)
                       gcAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGAGATACGGACCATATCATCGCCGTACACGCTGCTGATACTCATTAAACGCAGCCAACTTTCTGTCTCGACAATCCTTATATCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGCTTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGt                                                                                                                                                                                                                         >  2:616250/1‑251 (MQ=255)
                               aCCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTAGCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTg                                                                                                                                                                                                                 <  1:397925/251‑1 (MQ=255)
                                    aTCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCGGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAgt                                                                                                                                                                                                            <  2:26548/251‑1 (MQ=255)
                                                  gTTTTGCCTCCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAGCGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCGTTGCGATTGGTCAGTGATGCTGTCAATCAGAGTGGGAGTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCGGGAACACATCTCTGGAGATTTATGGCAGTTTTGCAAGTGTTACATATTcc                                                                                                                                                                                                <  1:574957/249‑1 (MQ=255)
                                                       gCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACg                                                                                                                                                                                           >  1:234575/1‑249 (MQ=255)
                                                           cccAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATAATCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGAAACTGTTACATATTCCGGACGAGAGGc                                                                                                                                                                                     <  2:137281/251‑1 (MQ=255)
                                                               gtgtgCGATAAGGACAATATAATCGCCGTACAAGCAGCTGATAATCATTAAACGCAGCCAAATTTCTGTATCGACCCTCCTTATCTCCCTGCCATAAGCAGCCTTAGCATTCTTTGCGATTGGTCAGTGATGCTGTCAATCAGTGGGGGATTTTTCTAGAATAGAAGATATAATCTTTCTAATTCTTGAACACATCTATGGAGATTTATGTCAGTTTTTCAAGTGTTACATATTCCGTACGAGCGTCTTCg                                                                                                                                                                                 >  2:512967/5‑251 (MQ=255)
                                                                 gAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAATCGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACAGATTCCGGACGAGCGGCTTCGc                                                                                                                                                                                <  1:219317/250‑1 (MQ=255)
                                                                     gATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGt                                                                                                                                                                           >  1:502396/1‑251 (MQ=255)
                                                                     gATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGt                                                                                                                                                                           >  1:465283/1‑251 (MQ=255)
                                                                              cATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAAcc                                                                                                                                                                  >  2:509113/1‑251 (MQ=255)
                                                                                              aaGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAAt                                                                                                                                                  >  2:25368/1‑251 (MQ=255)
                                                                                                       aTACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGCTATGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCTGAAGAGCGGCTTCGCAATGTTGCTAAACAGGTAGAAGACGAGAATGCAGAAAtt                                                                                                                                         >  2:447545/1‑251 (MQ=255)
                                                                                                         aCTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCa                                                                                                                                       <  1:487501/251‑1 (MQ=255)
                                                                                                            cATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCCCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCATAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAAGTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTTCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTTAATGCAGAAATTCAGCg                                                                                                                                    <  2:415216/251‑1 (MQ=255)
                                                                                                                aaaCGCAGCCAAATTTCTGGATCGACCATCCTTAGCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGGTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTAtc                                                                                                                                <  1:226799/251‑1 (MQ=255)
                                                                                                                      aGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCgat                                                                                                                          >  2:318563/1‑251 (MQ=255)
                                                                                                                        ccAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGCGGGGGATTTGTCTAGAATAGAAGGAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCgatg                                                                                                                         >  1:578484/1‑250 (MQ=255)
                                                                                                                                       aaacatcattctctcACTGCCATCAGAAGCCTTAGCAATCTTAGCGATTGGTCCGAGATGCTGTCCATCAGCGGGGGATTTGTATAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGCGCGGCTTCGCAAAGTTGATAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGGATCGTCGATGATATGTTCGAGACGAt                                                                                                         >  2:423673/12‑251 (MQ=255)
                                                                                                                                             ccTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGATATTCAGCGTATCGGCGAGGAGATGTTCGAGACGATGTAc                                                                                                     <  2:329443/249‑1 (MQ=255)
                                                                                                                                               ttATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGc                                                                                                   >  2:330299/1‑249 (MQ=255)
                                                                                                                                                   ctcCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCaga                                                                                                >  2:78292/1‑248 (MQ=255)
                                                                                                                                                     cccTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCgatgat                                                                                                                       >  1:348700/1‑223 (MQ=255)
                                                                                                                                                     cccTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCgatgat                                                                                                                       <  2:348700/223‑1 (MQ=255)
                                                                                                                                                         gCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTAt                                                                                       <  1:565395/251‑1 (MQ=255)
                                                                                                                                                                 cagcCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTgg                                                                               >  1:509344/1‑251 (MQ=255)
                                                                                                                                                                   gcattatcAATCTTTGCTATTGGTCAGTCATGCAGTCAATCAGAGGGTGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTTAACCTATGTATGGAGATTTATGTCAGTTTTGCAAGTGGTACATATTCCTTACGAGCGGCTTCTCAAATTTGCTATACCGGTAGAAAAAGTGAATGCAGAAATTCAGCGTATCGTCGATAATAGGTTCGAGACGATGTACGCAGACGAAGGTATTGGCCTggcg                                                                             <  2:346783/244‑1 (MQ=255)
                                                                                                                                                                       tAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGGACGCAGAAGAAGGTAGTGGCCTGGCGGCa                                                                          <  2:410516/250‑1 (MQ=255)
                                                                                                                                                                         gCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAAcc                                                                       >  1:534033/1‑251 (MQ=255)
                                                                                                                                                                                       tGGTCAGTGATTCTGTAAATCAGAGGGGGATTTGTATATAAAAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCGGGAGATTTATGTCATTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTGGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGGATCGTCGATGATATGTTCGAGACGATGTACGCCGAAGAAGGGATTGGCCTGGCGGCAACCCAGGTTGATATcc                                                          <  2:476087/250‑1 (MQ=255)
                                                                                                                                                                                                aTGCTGTCAATCAGAGGGGGATTGGTCTAGAATAGAAGAAATAAACTTTCTAATTCCTGAATACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGATTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTGCGATACGATGTACTCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATTCCCATCAACGTAt                                                <  2:104748/251‑1 (MQ=255)
                                                                                                                                                                                                     gTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTg                                           >  1:214469/1‑251 (MQ=255)
                                                                                                                                                                                                       cAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAAGTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACTTATCATTGtt                                         <  2:100676/251‑1 (MQ=255)
                                                                                                                                                                                                        aaTCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGtt                                         >  2:530366/1‑250 (MQ=255)
                                                                                                                                                                                                           cAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTACTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCATATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTg                                     >  2:622203/1‑251 (MQ=255)
                                                                                                                                                                                                                   gATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCg                             >  1:475334/1‑251 (MQ=255)
                                                                                                                                                                                                                      ttGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCggtg                          >  1:433945/1‑249 (MQ=255)
                                                                                                                                                                                                                      ttGTCTAGAATAGAAGAAATAATCGTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCAGTGTTATTGATGTTTCgg                            <  2:433945/249‑1 (MQ=255)
                                                                                                                                                                                                                       tGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGaa                          <  2:323934/250‑1 (MQ=255)
                                                                                                                                                                                                                            agaatagaaGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAAcc                      <  2:465283/249‑1 (MQ=255)
                                                                                                                                                                                                                            agaatagaaGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACACAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAAc                       <  2:578484/248‑1 (MQ=255)
                                                                                                                                                                                                                            agaatagaaGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATGGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAAcc                      <  2:351341/249‑1 (MQ=255)
                                                                                                                                                                                                                                 agaagaAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACg                <  1:381600/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aagaaATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAAc             <  2:583273/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAgtgt    >  2:106734/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ttCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTa  <  1:298247/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ttCTAATTCCTGAACACATCTCTGGAGATTTATGGCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCGTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTa  <  2:552840/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 cTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATAGTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCATCCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTa  <  2:172856/248‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
CCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTA  >  NC_000913/3433420‑3433908

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