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% 238.2 / NA 71F297L (TTC→TTAfdoGformate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (36/34);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTG  >  NC_000913/4084687‑4085177
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
aTGTTTTCAACCACATCTGGCGTGTCGCGGGAAACGTGCTGTTTCTGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCTGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAAttt                                                                                                                                                                                                                                                  >  2:43037/1‑251 (MQ=255)
       cAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTAATTg                                                                                                                                                                                                                                           <  1:232055/251‑1 (MQ=255)
        aaCCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGt                                                                                                                                                                                                                                          <  1:335104/251‑1 (MQ=255)
                    cGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGGTTAATTTTTCATTGTTCAGCAGGTACa                                                                                                                                                                                                                              <  2:354180/251‑1 (MQ=255)
                     gtgtAGCGGGAACCGTGTTGTTTCAGAAAGTTCCACACGCAGTGCGGATGTTGCAGCGTGGTATCGTGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCTGGATGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGTCTGCTTCAAAGCCTTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTTTTCAGCCCGTACAg                                                                                                                                                                                                                             <  2:240752/251‑1 (MQ=255)
                         aGCGGGAAACGTGCTGTATAAGAAAGTTCCACACGCAGCGCGGGTGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGACCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGCACAGGCCATCTTCAAAGACGTAATCCTCAAGAACGATCAGGATGGCGTTGGTATAGGCTTCGACGTATTCGCGGTTTAATTTTTCATTGTTCAGCGGGTACAGc                                                                                                                                                                                                                            >  2:520141/1‑248 (MQ=255)
                         aGCGGGAAACGAGCTGTTTCAGCAAGTTCCACACGCAGCGCGGCTGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTCTTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAAt                                                                                                                                                                                                                         >  2:348710/1‑251 (MQ=255)
                               aaaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTAAAGCAATACGCCt                                                                                                                                                                                                                   >  1:283089/1‑251 (MQ=255)
                                   gTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAcg                                                                                                                                                                                                               >  2:453453/1‑250 (MQ=255)
                                       tGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGa                                                                                                                                                                                                            >  1:57109/1‑250 (MQ=255)
                                             aGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAAt                                                                                                                                                                                                     <  2:36423/251‑1 (MQ=255)
                                                  gTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAg                                                                                                                                                                                                <  1:110005/251‑1 (MQ=255)
                                                      cacaCGCAGCGCGGATGTTGCAGCTTGGTATCGCGGTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAcc                                                                                                                                                                                            <  1:372021/251‑1 (MQ=255)
                                                            cAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACCGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACg                                                                                                                                                                                      >  1:341491/1‑251 (MQ=255)
                                                             aGCGCGGATGTTGCAGCGTGGTATCGCGTTTCGCGAAGCCGTTTTCGTCCCGTTCATAAGTCCAGGAGGATTTCTCGTACTTGAGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTCTTCGCGGTTTACTTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAACGGAACg                                                                                                                                                                                      >  1:523195/1‑250 (MQ=255)
                                                              gcgcGGATGTTGCAGCGTGGTAGCGCGTTTGGCGAATCCGTTTTCGTCCAGGTCCTAAGGCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCGTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGaa                                                                                                                                                                                    <  2:518679/251‑1 (MQ=255)
                                                                 cgGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAt                                                                                                                                                                                   >  1:132708/1‑249 (MQ=255)
                                                                   gATGTTGCAGCGTGGTATCGCGGTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAggg                                                                                                                                                                               <  1:90658/251‑1 (MQ=255)
                                                                      gTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCa                                                                                                                                                                            <  2:195256/251‑1 (MQ=255)
                                                                           aGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTGGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGGCAGTACCGGAACGAATAGGGGCAtagta                                                                                                                                                                       <  2:35890/251‑1 (MQ=255)
                                                                               tGGTATAGCGTTTGGCGACGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTCGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTcc                                                                                                                                                                   >  2:112375/1‑250 (MQ=255)
                                                                                ggTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCa                                                                                                                                                                   <  1:499138/250‑1 (MQ=255)
                                                                                 gTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGc                                                                                                                                                                 >  1:519380/1‑251 (MQ=255)
                                                                                            ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccgccg                                                                                                                                                      >  2:261804/1‑251 (MQ=255)
                                                                                            ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGAAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccgccg                                                                                                                                                      >  1:82106/1‑251 (MQ=255)
                                                                                             gCGAAGCCGTTTTCGTCCAGTTCATAAGTCCCGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTCTAGGCTTCGGTGTATTCGCGGTTTAATTTTTCCTTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGAAGTCAGCTACCGCCGCCGt                                                                                                                                                     >  2:38938/1‑251 (MQ=255)
                                                                                                       tttCGTCCAGTTCATAAGTCCAGGCGGATTTATCGTATTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCCGGCTGGCGTTGGTATAGGCTTCGGTGTCTTCGCGGTTTACTTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAACGGCACGACTAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTaa                                                                                                                                             >  2:574195/1‑249 (MQ=255)
                                                                                                          cGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGCTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTCCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGa                                                                                                                                        <  1:420572/251‑1 (MQ=255)
                                                                                                            tCCAGTTCATAAGTCCAGGAGGCTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAgg                                                                                                                                      >  2:101156/1‑251 (MQ=255)
                                                                                                              cAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                                                    <  1:120319/251‑1 (MQ=255)
                                                                                                               aGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatc                                                                                                                                   >  1:150146/1‑251 (MQ=255)
                                                                                                                  tCATAAGTCCCGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCCGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAAGGGCATATTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatcgat                                                                                                                                >  2:406573/1‑251 (MQ=255)
                                                                                                                    aTAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATca                                                                                                                              <  1:359265/251‑1 (MQ=255)
                                                                                                                         tCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAACGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGCGGATCGATCACAATc                                                                                                                         >  2:585971/1‑251 (MQ=255)
                                                                                                                              gaggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCtt                                                                                                                    >  1:49524/1‑251 (MQ=255)
                                                                                                                                 gaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTAAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGCTCACAATCAGCTTcgc                                                                                                                 >  1:263473/1‑251 (MQ=255)
                                                                                                                                     tATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGTTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAATGCGAGGAGCGATCACAATCAGCTTCGCGcc                                                                                                              <  1:165771/250‑1 (MQ=255)
                                                                                                                                      aTCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCg                                                                                                             <  1:561713/250‑1 (MQ=255)
                                                                                                                                          tACTTGCGTTTTTCCGCGTAGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGGACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTtgtg                                                                                                        >  2:126607/1‑251 (MQ=255)
                                                                                                                                              tGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAAtt                                                                                                    >  2:504863/1‑251 (MQ=255)
                                                                                                                                                 gTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTg                                                                                                 >  2:518869/1‑251 (MQ=255)
                                                                                                                                                    tttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGCTATAGCCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGCAACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGCATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCt                                                                                              >  1:239875/1‑251 (MQ=255)
                                                                                                                                                      tCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAAGCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCtt                                                                                             <  1:7257/250‑1 (MQ=255)
                                                                                                                                                               tAGCCGGTGCACAGGCCATCTTCAAAGCCGTAATCATCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGcccc                                                                                   >  1:124726/1‑250 (MQ=255)
                                                                                                                                                                       gAACAGGCCATCTTCAAAGCCGTAATCCTCACGCGCGTTCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGTTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACTCCTGACAGGAAAGCAATGTCAGTGCCGGAACGAATAGGGGCATAGTAGTCACCTACCGCCGCCGTACGCGTAAAGCTAGGATCGATCACAATCAGCTTCGCGCCTTTGTGAATTTTGGCTTCCATCGCCCATCGGAAc                                                                            <  2:558285/250‑1 (MQ=255)
                                                                                                                                                                         aCAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGt                                                                                                                                               <  1:525002/181‑1 (MQ=255)
                                                                                                                                                                         aCAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGt                                                                                                                                               >  2:525002/1‑181 (MQ=255)
                                                                                                                                                                                aTCTTAAAAGCCGTAATCCTACCGCCCGATCAACATGGCGTTAGTATAGGCTTCGGAGAATACGCGCTTTAATTTTTCCTTGTTCAGCAGGTACAGCAATAAGCCTGACAGGAAAGCAATGTCAGTCCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTCAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGAACCCCGCCCgggg                                                                  >  2:549715/1‑250 (MQ=255)
                                                                                                                                                                                    tCAAAGCCGTAATCCTCACGCACGATCATGCTGGCGGGGTTATAGGCTTCGGTGAATGCGCGGTTTAATTTTTCATTGTGCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTATTCATCTATCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTGCCATCGCCCAGCTGATCCCGAccttgtgag                                                               <  2:241553/250‑8 (MQ=255)
                                                                                                                                                                                      aaaGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCtt                                                            >  2:230952/1‑251 (MQ=255)
                                                                                                                                                                                      aaaGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCtt                                                            <  1:585971/251‑1 (MQ=255)
                                                                                                                                                                                       aaGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTCCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTc                                                           <  1:118169/251‑1 (MQ=255)
                                                                                                                                                                                          ccGTAATCCCCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTCCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgc                                                        <  2:509718/251‑1 (MQ=255)
                                                                                                                                                                                               aTCCTCACGCACGATCCGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTCATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgc                                                                                                                 >  2:60925/1‑189 (MQ=255)
                                                                                                                                                                                               aTCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgc                                                                                                                 <  1:60925/189‑1 (MQ=255)
                                                                                                                                                                                                 ccTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGtt                                                  <  2:163440/250‑1 (MQ=255)
                                                                                                                                                                                                  ccccctcacgTTCAGGCTGGCGTTGGTATAGGCTGCGGGGTATTCGCGGTTTAATTTTTCATTGTACAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGTACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGTCGTTAc                                                <  2:471956/245‑1 (MQ=255)
                                                                                                                                                                                                        cacgATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGccc                                           >  1:426920/1‑250 (MQ=255)
                                                                                                                                                                                                            aTCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACGGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCATCGGCGTTACCGCCCATcac                                      <  1:453453/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCacgacg                               >  2:106094/1‑249 (MQ=255)
                                                                                                                                                                                                                      ggTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAgg                            <  2:457917/250‑1 (MQ=255)
                                                                                                                                                                                                                       gTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTCCCGGAACGAATGGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGt                           <  1:531176/251‑1 (MQ=255)
                                                                                                                                                                                                                       gTTGGTATAGGCTTCGGTGTATTAGCGGCTTAATTTTTCATTGTTAAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCCGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGCGGt                           >  2:354060/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGGCTTCGGTGTATTCGCGGTTTCATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTc                   >  2:375145/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggctggggTGTAGTCGCGGTTGAACTTTTCAATGTTCAGCAGGTACGGCAATCAGCCTGACGGGAAAGCAATGTCAGTGCCGGAGGTAATAGGGGCACAGTAGTAATCGTCCGCCGGCGGACGCGTAAAGCGAGGAACGATCACAAGCAGCTTCGCGCCGTTGTGAATTTTGGCGTCCATCGCCCAGCGGCAACCGACCGGGTGCGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTGGGCGTTCt                  <  2:136219/245‑1 (MQ=255)
                                                                                                                                                                                                                                       gtgtATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTc           <  1:54706/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGAc        >  1:400944/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tCGCGGTTTAATGTTTCATGTTTCAGCAGGTACAGCAATACTCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCTCCCATCACCACGACGATGTTGGCGTTCTTGATGTCGAccc      <  1:316620/250‑1 (MQ=255)
                                                                                                                                                                                                                                                cgGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATca                                                                                                                              <  1:401207/127‑1 (MQ=255)
                                                                                                                                                                                                                                                cgGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATca                                                                                                                              >  2:401207/1‑127 (MQ=255)
                                                                                                                                                                                                                                                cgGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGGCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGtg  <  1:230952/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTG  >  NC_000913/4084687‑4085177

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