Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,084,932 G→T F297L (TTC→TTA fdoG ← formate dehydrogenase O subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,084,9320GT98.6% 228.2 / NA 70F297L (TTC→TTAfdoGformate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (36/32);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCC  >  NC_000913/4084686‑4085173
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
gatgatTTCCACCACATCTGGCGTGTAGCGGGAAACGCGCTGTTTCAGCACGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGAGAAGCCGTTTTCGTCCAGTTCATAAGTACAGGAGGATATATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCATTCTTCAAAGCCGTAATCCTCACGCCCGATCAGGATGGCGTTGGCATCGGTTTCGGTGTTTTCGTGGTTTAAtt                                                                                                                                                                                                                                               >  2:191913/1‑251 (MQ=255)
        cAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGCTGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCAttt                                                                                                                                                                                                                                       >  1:550948/1‑250 (MQ=255)
        cAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTg                                                                                                                                                                                                                                       >  1:564526/1‑251 (MQ=255)
                   gtcGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCCCACGCAGCTCGGATGTTGCAGCGTTGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGGAGGATTTTTCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACGGGCCATCTTCAAAGCCGTAATCCTCACTCACTATCAGGCTGGCGTTGGTATAGGCTTCGGGGTATTCTCGGTTGAATTTTTCACTGTGcagcag                                                                                                                                                                                                                               <  2:370619/246‑1 (MQ=255)
                   ggCGTGTAGCGGGAAACGTGCTGGTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTcagcag                                                                                                                                                                                                                               <  1:500910/248‑1 (MQ=255)
                     cGTGTAGCGGGAAACGTGCTGTTTCAGCATGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACa                                                                                                                                                                                                                          <  1:335331/251‑1 (MQ=255)
                        gtAGCGGGAATCTTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTTGTATAGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGc                                                                                                                                                                                                                        <  1:462584/250‑1 (MQ=255)
                        gtAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGc                                                                                                                                                                                                                        >  1:468562/1‑250 (MQ=255)
                              tttaaCGTGCTATTTCAGCCAGTTCCACACGCAGCGCGGCTGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCCTCCAGTTCATAAGTCCAGGAGGCTTTATCGTACTTGCGTTTTTCCGCGTCTTAGCCGGTGAACAGGCCATCTTCAAGTCCGTACTCCTCACTCACGACCAGGCTGGCGTTGGAGTATGCTTCGGTGACTTCGCGGGTTAGTTTTTCATTTTTctgcagtgacagccataaac                                                                                                                                                                                                                 >  2:304160/4‑237 (MQ=255)
                              ggtaaCGTGCTGTTTCAGCAAGTTCCACACGCTGCTGGGTGGTTGCAGCGTGGTAGCACGTTTGGCGACGCCGTTTTCGTCCAGTTCATAAGCCCAGGAGGATTTCTTGTACGTGAGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCCCGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCTGTTTAATTTTTCATTGTTCAGCAGGCACAGCAATACGc                                                                                                                                                                                                                 <  2:339171/248‑1 (MQ=255)
                                aaaCGTGCTTTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                                                               <  2:170965/251‑1 (MQ=255)
                                aaaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                                                               <  1:67798/251‑1 (MQ=255)
                                    gTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACa                                                                                                                                                                                                           <  1:404653/251‑1 (MQ=255)
                                         gTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGa                                                                                                                                                                                                        <  1:473617/249‑1 (MQ=255)
                                                  aGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCa                                                                                                                                                                                             <  1:347031/251‑1 (MQ=255)
                                                     tCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCCTTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGt                                                                                                                                                                                           >  1:457365/1‑250 (MQ=255)
                                                      ccACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAc                                                                                                                                                                                         <  1:478585/251‑1 (MQ=255)
                                                                  cgGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCTTTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGGACGAAt                                                                                                                                                                               >  2:200715/1‑249 (MQ=255)
                                                                           cAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTACAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCCGGATGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGGACAGCAATAAGCCTGACAGGAAAGCAATGTCAGTCCCGGAACGAATAGGGGCAtagt                                                                                                                                                                    >  2:318640/1‑251 (MQ=255)
                                                                                      cgcgTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAcc                                                                                                                                                         >  1:318017/1‑251 (MQ=255)
                                                                                      cgcgTTTGGCGAAGCCGGTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAcc                                                                                                                                                         <  1:507453/251‑1 (MQ=255)
                                                                                       gcgTTTGGTGAAGCCGTGTTCGTCCAGTTCATAAGTCCAAGAGGATTTATCGTACTAGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccg                                                                                                                                                        <  2:159642/251‑1 (MQ=255)
                                                                                             ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGGTCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccgccg                                                                                                                                                  >  1:148592/1‑251 (MQ=255)
                                                                                             ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccgccg                                                                                                                                                  >  1:397823/1‑251 (MQ=255)
                                                                                                 aaGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGGGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcg                                                                                                                                              <  1:386785/251‑1 (MQ=255)
                                                                                                      gTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTa                                                                                                                                          >  1:155079/1‑250 (MQ=255)
                                                                                                         ttCGTCCAGTTCATAAGTCCAGTAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCACTACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATATTAGTCAGCTACCGCCGCCGTACGCGTAAAg                                                                                                                                       <  1:25405/250‑1 (MQ=255)
                                                                                                          tcctccattcaaTAAGTCCAGGTGGTGTTTTCGTATTTGCTTTATCCCCCGTCGTCGCGGGTGAACAGGCCTGCTTCAAAGCCGTAATCCTCACGCACGACCAGGCTGTCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTCCAGCAGGTCCAGCAATCCGCCTGACAGGAAACCTATTTCAGTACCGGAACGCATAGGGGAATATTTGGCAGCTACCGCCGCCGTACGCGTAAAGc                                                                                                                                      <  2:383398/239‑1 (MQ=255)
                                                                                                            gTCCAGTTCATAAGTCCAGGAGAATTTATCGTACTTGCGATGTTCCGCTGCGTGGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTTTAGGCTTCGGTGTATTCGCGTTTAAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAAAGGGGCATAGGAGTCAGCTACCGCCGCCGTACGCGTAAAGCGa                                                                                                                                    <  2:481381/250‑1 (MQ=255)
                                                                                                              caagtaaaTAAGTCCAGGAGGCTTTATCGTACTTGAGATTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCATGACAGGAAAGCAATGTCAGTACCGTATAGGATAGGGGCATCGAAGTCAGCTAACGCCGCCGTACGCGTAAAGCGAgg                                                                                                                                  >  2:234343/8‑250 (MQ=255)
                                                                                                                       aaGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCaa                                                                                                                                                                                                                      >  2:132199/1‑157 (MQ=255)
                                                                                                                       aaGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCaa                                                                                                                                                                                                                      <  1:132199/157‑1 (MQ=255)
                                                                                                                          tCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATc                                                                                                                     >  1:402671/1‑251 (MQ=255)
                                                                                                                          cccAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATc                                                                                                                     <  1:157930/250‑1 (MQ=255)
                                                                                                                           ccAGGAGGATTTATCGTACTGGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGAACGCGTAAAGCGAGGATCGATCACAATc                                                                                                                     <  1:570195/250‑1 (MQ=255)
                                                                                                                             aggaggATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGc                                                                                                                  >  1:491548/1‑251 (MQ=255)
                                                                                                                              ggaggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTATTCACCTACCGCCTCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCt                                                                                                                 <  1:193134/251‑1 (MQ=255)
                                                                                                                                   aTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTc                                                                                                               <  1:234373/248‑1 (MQ=255)
                                                                                                                                   aTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTc                                                                                                               >  2:234373/1‑248 (MQ=255)
                                                                                                                                            aCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTtgtg                                                                                                    >  1:395678/1‑250 (MQ=255)
                                                                                                                                              ttGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAAt                                                                                                 >  1:34147/1‑251 (MQ=255)
                                                                                                                                                     tttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAt                                                                                                                                                                               >  1:164960/1‑166 (MQ=255)
                                                                                                                                                         cgcgTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGACTTTTGGCTTccc                                                                                      >  2:116187/1‑250 (MQ=255)
                                                                                                                                                         cgcgTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCa                                                                                      <  1:249587/251‑1 (MQ=255)
                                                                                                                                                            gtcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCAtag                                                                                                                                                                     >  1:386579/1‑169 (MQ=255)
                                                                                                                                                            gtcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCAtag                                                                                                                                                                     <  2:386579/169‑1 (MQ=255)
                                                                                                                                                            gtcgtAGCCGGTGAAAAGGCCATCTTCAAAGCCGTAATCCTCCCGCACGATCAGGCTGGCGTTGGTTTAGGTTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCg                                                                                   <  1:435974/251‑1 (MQ=255)
                                                                                                                                                                 agccgtggAACAGGCCTTCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTGCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGACAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCa                                                                               <  2:419957/243‑1 (MQ=255)
                                                                                                                                                                 aGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGccc                                                                                >  1:285068/1‑249 (MQ=255)
                                                                                                                                                                  gCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGc                                                                             >  1:406966/1‑251 (MQ=255)
                                                                                                                                                                     ggTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCgg                                                                           >  1:403476/1‑250 (MQ=255)
                                                                                                                                                                                     ccaacgCCGTAATCCTCACGCACGGGCAGGCGGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATCCGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTCCCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGc                                                          <  2:148592/246‑1 (MQ=255)
                                                                                                                                                                                        aggccttAATCCTCACGCACGATCAGGCTGGCGTTGGTTTAGTCTTCTGGGTATTCGCGGTTTTATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGTCAGGAAAGCAATTTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTCCCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGcgtc                                                       <  2:517878/245‑4 (MQ=255)
                                                                                                                                                                                        aaGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTc                                                       <  1:105614/251‑1 (MQ=255)
                                                                                                                                                                                         agcggTAATCCTCACGAACGATCTGGCTGGCGTTGGTATAGGCATCGGGGTGATCGCGGTTTAATTTTTCATTGTTCGGCAGGGACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAGAGAATAGGGGCATAGTAGACAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAAGCAGCTTCGCGCCGTGGGGAAATTTGGCTTCCATCGCCCAGAGGAACCCGACCGGGAGAGCTTc                                                       >  2:436546/5‑250 (MQ=255)
                                                                                                                                                                                               aaTCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcggcg                                                <  1:339674/251‑1 (MQ=255)
                                                                                                                                                                                                  ccTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTa                                             <  1:191280/251‑1 (MQ=255)
                                                                                                                                                                                                      acgcacgATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGAACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGACTTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGc                                         >  1:517062/1‑251 (MQ=255)
                                                                                                                                                                                                       ggcacgAACAGGCTGGCGTTGTTATAGGCCTCGGTGTATTCGCTGTTTAATTTTTCATTTTTCAGCAGGTACAGCAATACTCCTCACAGGAACGCAATGTCAGTCCCGGAACGAATGGGGGAATCGTAGTCAGTTACCGCCGCCGTACGCGTAAATCGAGGATCGATCTCAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCTGAACCCGACCGGGTGAg                                                           <  2:552301/231‑1 (MQ=255)
                                                                                                                                                                                                       cgcacgATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAg                                                           >  1:552301/1‑232 (MQ=255)
                                                                                                                                                                                                           cgATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATc                                    <  1:77146/251‑1 (MQ=255)
                                                                                                                                                                                                               cAGGCTGGAGTAGGAATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGGACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGACTAGGGGCAAAGTAGTCAGCTACCGCCGCCGTACGCGTAACGCGAGGATCGATCCAAATCAGATTCGCGCCGTTGTGAGTTTTGGCTTCCATCGCCCAGCGGAACCCGTCCGGGGGAGCTTCAGAGTCGTTCCCGCCCATaaccc                                >  2:409053/1‑250 (MQ=255)
                                                                                                                                                                                                                aGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATcaccac                               <  1:15655/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGTTGGTATAGGCTACGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGGAGTCAGCTACCGCCGCCGGACGCGTAAAGCGAGGATAGATCACAATCAGCTTCGCGCCGTTGCGAATTTTGGCTTCCATCGCCCAGCGGACCACGACCGGGTGAGCATCAGCTGCGTTACCGCCAATCACCacgtagagg                        >  2:394111/1‑245 (MQ=255)
                                                                                                                                                                                                                                aGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTc               >  1:539897/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTAACGCCCATCACCACGACGAGGTTGGCGTTCtt             >  1:431822/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAACGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCCTCACCACGACGCGGTTGGCGTTCTTg            >  2:493618/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCAATCGCCCAGCGGAACCCGACCGGGGGAGCTTCAGCGGCGTTACCGCCCCTCCCCACGACGACGTTGGCGTTCTTGATg         >  2:95942/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ttCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGCGGTTGGCGTTCTTGATGTCGAccc  >  2:110875/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ttCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATAACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGAccc  >  1:336242/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
GATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCC  >  NC_000913/4084686‑4085173

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