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,6650CT97.6% 127.7 / ‑4.3 43intergenic (‑105/‑25)smf/defprotein Smf/peptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/21);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

GCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCT  >  NC_000913/3433419‑3433901
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
gccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTcc                                                                                                                                                                                                                                           <  1:60716/250‑1 (MQ=255)
 ccgaaaTGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTg                                                                                                                                                                                                                                         <  2:127034/247‑1 (MQ=255)
  cGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGa                                                                                                                                                                                                                                        <  1:188805/251‑1 (MQ=255)
        gTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACAt                                                                                                                                                                                                                                  <  1:195284/251‑1 (MQ=255)
                        gcAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTAt                                                                                                                                                                                                                    <  1:386002/249‑1 (MQ=255)
                          aaTACAACCGCATAAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTTGAGATTTATGTCa                                                                                                                                                                                                                <  1:107517/251‑1 (MQ=255)
                            tACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGt                                                                                                                                                                                                              >  1:237743/1‑251 (MQ=255)
                            tACAACCGCATCAATATGCGACTGTTTTGCCACCCAGCGCGCGATACGGACAATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCATTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGCGATTTATGTCAGt                                                                                                                                                                                                              >  2:286403/1‑251 (MQ=255)
                                  cGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATATCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGc                                                                                                                                                                                                         >  2:194631/1‑250 (MQ=255)
                                    cATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGaa                                                                                                                                                                                                                                       >  1:261158/1‑218 (MQ=255)
                                                        gCCACCCGGTGAGGAATACGGACAATAATAGCGCCGGACAAGCAGCAGATATACATTAAACACAGCAAATTTTCTGTATCTACAAACCTTATCTCCCTGCCATAAGCAGCCTTAGCCATCTTTGCGATTGGTCAGTGATGCTTGCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAAGTCATGAACACAGCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATAGTCCGGAc                                                                                                                                                                                     <  2:113828/248‑1 (MQ=255)
                                                                     cGATACGGACAATTTCATCGCCGTACAACCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGGCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAg                                                                                                                                                                     <  2:98256/251‑1 (MQ=255)
                                                                           ggACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCt                                                                                                                                                                >  1:174442/1‑250 (MQ=255)
                                                                                  tccatcCCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCgg                                                                                                                                                         <  2:174442/248‑1 (MQ=255)
                                                                                     agcGCCTTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTTTCGACCAGCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTaga                                                                                                                                                     <  2:54596/249‑1 (MQ=255)
                                                                                      cccccGTAAAAGCTGCTGATCCTCATTAAACGCAGCCAAATTTCGGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTagaa                                                                                                                                                    <  2:237743/248‑1 (MQ=255)
                                                                                                        aTACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTg                                                                                                                                                                  >  1:462077/1‑219 (MQ=255)
                                                                                                        aTACTCATTAAACGCAGCCAAATTTCTGTATCCACCAACCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTg                                                                                                                                                                  <  2:462077/219‑1 (MQ=255)
                                                                                                             cATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCg                                                                                                                             <  1:2550/251‑1 (MQ=255)
                                                                                                                                        gcccatcattttcccccTTCCATAAGCAGCCTTAGCAAGCTTTGCGATTGGTCAGTGAAGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTGCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCTGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCTGCGATGATATGTTCGAGACGAt                                                                                                  <  2:331120/237‑1 (MQ=255)
                                                                                                                                         ccctttcttttgccccTGCCCAAAGCAGCCTTAGCAAGCTGTGCTATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTTGAATAGAAGAAAGAATCTTTCTAATGCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGATGAAGCGAATGCAGAAATTCATCGTATCGGCGATGATATGGTCGAGACGATg                                                                                                 <  2:323209/238‑1 (MQ=255)
                                                                                                                                                      tccgtccATAAGCAGCCTTAGCAATCTTTTCGATTGGTCAGTGATGCTGTCAATCAGCGGGGGATTTGGCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTGTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCTGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGTTATGTTCGAGACAATGTACGCagaaga                                                                                      <  1:316465/244‑1 (MQ=255)
                                                                                                                                                         tGCCATAAGCAGCCTTAGCAATCTTTGCGCTTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAgga                                                                                  >  1:379311/1‑249 (MQ=255)
                                                                                                                                                             aTAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTgg                                                                             >  1:22104/1‑251 (MQ=255)
                                                                                                                                                                agcagcCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGGCAAGCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTGCGAGACTATGTACGCAGAATAAGTTATTGGCCt                                                                          <  2:322543/251‑1 (MQ=255)
                                                                                                                                                                     ccTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTa                                                                                                                                                                                                 <  1:198475/127‑1 (MQ=255)
                                                                                                                                                                     ccTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTa                                                                                                                                                                                                 >  2:198475/1‑127 (MQ=255)
                                                                                                                                                                                   gCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGat                                                       <  1:250167/251‑1 (MQ=255)
                                                                                                                                                                                                 aTGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAAGTTTTCTACTTCCTGAACACATCTCGGGAGATTTATGGCAGTTTTGCAAGTTTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACGGGGAGAAGAAGTGAATGCCGAAATTCAGCGTAGCGTCGATGATAGGTTCGCGACGATGTACGCAGAAGAAGGTATTGGCCTGTGGGCAACCCAGGTTGATATCCATCAACGTAt                                         <  2:440344/251‑1 (MQ=255)
                                                                                                                                                                                                    cTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCAt                                      >  1:382699/1‑251 (MQ=255)
                                                                                                                                                                                                     tGTCAATCAGAGGGGGCTTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCCAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGCAATTCAGCGTATCGTCGATGATATGTTCGAGACTATGTACGCAGAAGAAGGTATTGGCCAGGCGGCAACCCAGGTTGATATCCATCAACTTATCAtt                                     >  1:76389/1‑251 (MQ=255)
                                                                                                                                                                                                            cAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCACCGTATCACTGTTAtt                               <  1:235356/250‑1 (MQ=255)
                                                                                                                                                                                                            cAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTAGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTg                              >  1:459963/1‑251 (MQ=255)
                                                                                                                                                                                                             agagGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGAATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTgt                             >  2:16768/1‑250 (MQ=255)
                                                                                                                                                                                                             agagGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGAAAGTGTTACATATTCCGGACGAGCGGCTTAGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTg                              >  1:229086/1‑250 (MQ=255)
                                                                                                                                                                                                                 ggggaggtggcTAGAATATAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGCGATTTATGTCAGTTATGCAAGTGTTACCTATTCCGGACGAGCGGCTTCGAAAAGTTGTTCAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGGTGATATGTTCGTGACGATGTCCGCGGACACAGGTATTGGACTGGCGGCAACCCAGGTTGATATTCATCAACGTATCATTTTTAATGATGtt                         >  2:39043/11‑251 (MQ=255)
                                                                                                                                                                                                                    gATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGGAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGATTCATTGTTATTGATGTTTCg                      >  2:492797/1‑251 (MQ=255)
                                                                                                                                                                                                                        tGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATg                                                                                                             >  1:156081/1‑160 (MQ=255)
                                                                                                                                                                                                                        tGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATg                                                                                                             <  2:156081/160‑1 (MQ=255)
                                                                                                                                                                                                                                 tagaaGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACg         >  1:142960/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tagaaGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACg         >  1:147121/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aagaaATAATCTTTCTAATTCCTTAACACATCTCTGGAGATTTATTTCAGTTTTGCAAGTTATACATATTCCGGTCGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTTAGCTTATCGTCGATGATATTTTCGAGACGATGTACGCAGAAGAAGGTATTGTCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAAc      <  2:447338/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aataatCTTTCTAATTCCTGAACACATCTCTGGCGATTTATGTCCGTTTTGCAAGTGTTACTTATTCCGTACGAGAGGCTTCGCAAAGTTGCTAAACCGGTAGAAGTAGTGAATGCAAAAAATCAGCGTATCGTCGATGATATGTTCGATACGATGTACGCAGAAGAAGGTATTGGACTGGCGGCAACCCAGGTGGATATCCATCAAAGAATCATTGCTATTGATGTTTCGGAAAACCGGTACGAACggtt  >  2:454452/1‑249 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
GCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCT  >  NC_000913/3433419‑3433901

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