Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,084,903 T→G N307T (AAC→ACC)  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,9030TG100.0% 251.9 / NA 70N307T (AAC→ACC) fdoGformate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (31/39);  total (31/39)

GGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACG  >  NC_000913/4084661‑4085148
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
ggAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAg                                                                                                                                                                                                                                               <  1:433192/251‑1 (MQ=255)
     gcgTCTTTTGGCGTACCACAGATGTTTTCACCCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGTGGTGGTATAGGCTTc                                                                                                                                                                                                                                          <  2:182789/251‑1 (MQ=255)
       gTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCgg                                                                                                                                                                                                                                        >  1:314972/1‑251 (MQ=255)
           tttGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGgtgt                                                                                                                                                                                                                                     <  1:33881/250‑1 (MQ=255)
           tttGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGgtgt                                                                                                                                                                                                                                     >  1:251555/1‑250 (MQ=255)
              ggCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTAtt                                                                                                                                                                                                                                  <  1:239103/250‑1 (MQ=255)
              ggCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTc                                                                                                                                                                                                                                 <  1:201643/251‑1 (MQ=255)
              ggCGTACAACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATAGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGATTCGGTGTATTc                                                                                                                                                                                                                                 >  1:312059/1‑251 (MQ=255)
                     cacaGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTg                                                                                                                                                                                                                          <  1:381073/251‑1 (MQ=255)
                         gATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCTTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGGATTCGCGGTTGAAtt                                                                                                                                                                                                                      <  1:146242/251‑1 (MQ=255)
                              tttCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTc                                                                                                                                                                                                                  <  1:207420/250‑1 (MQ=255)
                               ttCAACCACATCTGGCGTGTTGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTGTTTCCGCGTCGTAGCCGGTTAACAGGCCATCTTCAAAGCCTTAATCCTCACGCACTATCAGGCGGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCa                                                                                                                                                                                                                 <  1:221931/250‑1 (MQ=255)
                                 cAACCACATCTGGCGTGTAGCGGGAAACTTGCTGTTTCAGCAAGTGCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTg                                                                                                                                                                                                              <  2:186157/251‑1 (MQ=255)
                                       caTCTGGCGTGTAGCGGGAAACGTGCTTTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTcagc                                                                                                                                                                                                        <  1:328538/251‑1 (MQ=255)
                                       caTCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTcagc                                                                                                                                                                                                        <  1:253386/251‑1 (MQ=255)
                                            ggCGTGTGGCGGGAAACGTGCTGTTTCGGCAAGGTCCACACGCAGCGCGGATGTTGCAGCGTGTTTTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTTCTTTCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTGCAAAGCCGTAATCCTCCCGCTCGATCATGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGt                                                                                                                                                                                                    <  2:312059/250‑1 (MQ=255)
                                            ggCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACCCGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTa                                                                                                                                                                                                   <  1:357298/251‑1 (MQ=255)
                                            ggCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGGTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTACGGTGTATTCGCGGTTGAATTTTTCATTGGTcagcag                                                                                                                                                                                                      <  1:89842/248‑1 (MQ=255)
                                                   aGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAAt                                                                                                                                                                                            <  2:304597/251‑1 (MQ=255)
                                                     cGGGAAACGTGCTGTTTCAGCAAGATCCACACGCAGCGCGGATGTTGCCGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCAAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCCCGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCTTTGTTCAGCAGGGACAGCAATAc                                                                                                                                                                                          >  2:256418/1‑251 (MQ=255)
                                                         aaaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                                      >  1:357472/1‑251 (MQ=255)
                                                            cGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAc                                                                                                                                                                                   >  2:248866/1‑251 (MQ=255)
                                                              tGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGATGCAGCGTGGTATCGCGTTTGACGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACCGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATCGGCTTCGGTGTATTCGCGGCTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAc                                                                                                                                                                                   >  1:211623/1‑249 (MQ=255)
                                                                cTGTTTAAGCAAGTTCCACACGCAGCGCGGACGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAAGAGGATTTATAGTACTTGCGTTTTTCCGAGTCGTAGCAGGTGAACTGGCCATCTTCAAAGCCGAACTCCTCACGCACGATCAGGCTGGCGGTGGTAAAGGCTTCGGGGTATTCGCAGTTGAATTTTTCCTTGTTCAGTAGGTACAGCAATACCCCTGTCAgg                                                                                                                                                                                >  1:186514/1‑250 (MQ=255)
                                                                     tCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGACTTTTTCATTGTTACGCAGGTACAACAATACGCCTGACAGGAAAGc                                                                                                                                                                           >  2:76378/1‑250 (MQ=255)
                                                                          aaGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTc                                                                                                                                                                     >  1:424368/1‑251 (MQ=255)
                                                                           aGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAAt                                                                                                                                                                                                                       >  1:12767/1‑200 (MQ=255)
                                                                           aGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAAt                                                                                                                                                                                                                       <  2:12767/200‑1 (MQ=255)
                                                                              tCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGt                                                                                                                                                                  <  1:143061/250‑1 (MQ=255)
                                                                               ccACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCTATACGCCTGACAGGAAAGCAATGTCAGTa                                                                                                                                                                 <  2:306400/250‑1 (MQ=255)
                                                                               ccACACGCAGCGAGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCGGGCGGTGGTATAGGCTGCGGTGTATTCGCGGTTGAAGTTTTCATTGTTCAGCAGGTACAGCAAGTCGCCTGACAGGAAAGCAATGTCAGTAc                                                                                                                                                                <  2:48307/251‑1 (MQ=255)
                                                                                  cacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCgg                                                                                                                                                             >  1:315651/1‑251 (MQ=255)
                                                                                        gcgcGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGaa                                                                                                                                                       <  1:360630/251‑1 (MQ=255)
                                                                                         cgcgGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAt                                                                                                                                                      <  1:70894/251‑1 (MQ=255)
                                                                                          gcgGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAt                                                                                                                                                      >  1:210321/1‑250 (MQ=255)
                                                                                           gtgATGTTGCAGCGTTGTATCGCGTGTGGCGAAGCCGTTTTCGTCCAGTTCATATGTCCAGGAGGATTTATCGTACTTGCGTTGTTCCGCGTCGTAGCCGGTGAACAGGCCATCTACAGAGCCGTAATCCTCACGCACGATCTGGCTGGCTGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGGAGGTGCAGCAGTACGCGTGACAGGAAAGCAATGGCAGTACCGGAACGAATAg                                                                                                                                                    <  1:243307/249‑1 (MQ=255)
                                                                                                gTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGc                                                                                                                                                >  1:129362/1‑250 (MQ=255)
                                                                                                 ttGCAGAGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATAAGCCTGACAGGAAAGCAATGACAGTACCGGAACGAATAGGGGCAt                                                                                                                                              >  1:271314/1‑251 (MQ=255)
                                                                                                        gTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTc                                                                                                                                       >  1:52126/1‑251 (MQ=255)
                                                                                                                        cGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGt                                                                                                                        >  2:162014/1‑250 (MQ=255)
                                                                                                                          aaGCCGTTTTCGTCCAGTTCAGAATTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcg                                                                                                                     <  2:194192/251‑1 (MQ=255)
                                                                                                                             ccGTGTTCGTCCAGTTCAGAAGTCCAGGAGGATTTATCGTACTTGCGTTTTGCCGCGTCGTAGCAGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGGACAGCAATACGCCGGACAGGAAAGCAAGGTCAGTACCGGAACGAAGAGGGGCATAGTGGTCAGCTACCGCCGCCGTACGCGt                                                                                                                  <  2:160720/251‑1 (MQ=255)
                                                                                                                               gTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTa                                                                                                                 >  1:159951/1‑250 (MQ=255)
                                                                                                                                 tttagtccatgtcagaaGTCCAGGAGGATTTCACGTACTTGCGTTTTCCCGCGGCGTAGCCGGTGAACAGGCCAGTTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTGCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATAAGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATATGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTaaa                                                                                                               <  2:96500/235‑1 (MQ=255)
                                                                                                                                     gTCCAGTTCATAAGTCCAGGAGGATTTATCTTTCGTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTAAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTTACAGGAAAGCAATGTCAGTACCGGAACGATTAGGGGCATAGTGGTCAGCTACCGCCGCCGTACGCGTAAAGCGAg                                                                                                          <  2:315651/251‑1 (MQ=255)
                                                                                                                                       ccAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGAAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCCCGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCCTTGTTCAGCAGGGACAGCAATACGCCTGACCGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGGAGTCAGATACCGCCGCCGTACGCGTAAGGCGAggg                                                                                                        >  2:388207/1‑250 (MQ=255)
                                                                                                                                                aaGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATcac                                                                                                <  1:88980/250‑1 (MQ=255)
                                                                                                                                                 aGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGACCTGGCCTTCTTAAAAGCAGTAATCCTCACGCACGATCACGCTGGCGGTGGCATTGGCTTCGGTGTATTCGAGGATGAATTTTTCAt                                                                                                                                                                                                                >  2:191600/1‑137 (MQ=255)
                                                                                                                                                 aGTCCAGGAGGAGTTATCGTACTTGCGTTTTTCCGCGTCGTAGCGGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGAGCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCAt                                                                                                                                                                                                                <  1:191600/137‑1 (MQ=255)
                                                                                                                                                   tCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCCGGTACAGCAATGCGCCTGACAGGAAAGCAATGTCAGTACCGGCACGAATTGGGGCCTAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATc                                                                                            >  2:155314/1‑251 (MQ=255)
                                                                                                                                                     cAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTCGTCAGCTAACGCCGCCGTACGCGTAAAGCGGAGATGGATCACaaac                                                                                            >  2:116407/1‑247 (MQ=255)
                                                                                                                                                       ggaggaGTTATCGTACTTGCGTTTTTCCGCGTCGTGGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCGACCGCCGCCTTACGCGTAAAGCGAGGATCGATCACAATCAGCt                                                                                        <  2:211623/251‑1 (MQ=255)
                                                                                                                                                            aTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATAAGCTTcgcg                                                                                   >  1:253587/1‑251 (MQ=255)
                                                                                                                                                                   gTACTGGCTTTTTTCCGCGGCGTAGTCGTGGAACGGGCCATCTTCAAAGCCTTAATCCCCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTTTATTCTCGGTTGAAGTTGTCATGGTTCAGCAGGTACTGCAATACGCCTTACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATATTAGTCATCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCTCAATCAGCTTCGCGCCGTtgt                                                                            <  2:364746/251‑1 (MQ=255)
                                                                                                                                                                             ttttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGAGAGGATCGCTCACAATCAGCTTCGCGCCGTTGGGATTTTTGGc                                                                  >  2:298949/1‑251 (MQ=255)
                                                                                                                                                                             ttttCCGCGTCGTAGCAAGAACACAGGACATCTTCAAAGCCGTAATCCTCAAGCACGATCAGGCTGGCGGTGGAATAGGATTCGGTGTATTCGCGGTTGCATTTTTCATTGTTCAGCAGGTCCAGCAATACGCCTGGCAGGGACGCAATGTCAGTACCGGAACGCATAGGGGAATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGACCACAATAAGCTTCGCGCCGTTGTGACtttttgc                                                                  >  2:435808/1‑248 (MQ=255)
                                                                                                                                                                                tCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGAGAGGATCGATAACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTc                                                               >  1:264274/1‑251 (MQ=255)
                                                                                                                                                                                cccGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCCCAATCAGTCTTGCGGTGGTATAGGCTTCGGTGTATCCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTc                                                               <  1:376563/250‑1 (MQ=255)
                                                                                                                                                                                cccGCGGCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTACAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCATCTACCGCCGCCGTACGCGTAAAGCGAGGAGCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTc                                                               <  2:129362/250‑1 (MQ=255)
                                                                                                                                                                                  cgcgTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTcc                                                              >  1:405970/1‑250 (MQ=255)
                                                                                                                                                                                         tAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgcc                                                                                                                             <  1:103057/180‑1 (MQ=255)
                                                                                                                                                                                         tAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgcc                                                                                                                             >  2:103057/1‑180 (MQ=255)
                                                                                                                                                                                                   aCAGGCCATCTTCAAAGCCGTAATCCTCACGCCCGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCg                                            <  1:240185/251‑1 (MQ=255)
                                                                                                                                                                                                     aggtcgtattCAAAGCCGTAATCCTCACGGACGATTAGGCTGGCGGTGGTATAGGCTGCGGTGTATGCGCGGTTGAATTTTTCATTGTGCAGCTGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTTCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGTCGGACGCGTAAAGCGAGGATCGATCACAATCATCTTCGAGCCGGTGTGAAGTTTGGCTTCCATCGCCCAGCGGAACCCGAc                                          <  2:274694/243‑1 (MQ=255)
                                                                                                                                                                                                                           ccTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGGATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGtt                     <  1:347907/250‑1 (MQ=255)
                                                                                                                                                                                                                            cTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAACTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGACGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAACCGACCGGGTGAGCTTCAGCGGCGTTAc                   >  2:63205/1‑251 (MQ=255)
                                                                                                                                                                                                                                cgcacgATCAGGCTAGCGGTGGAATAGGCTTCGTTGAATTCTCGGTTGAATTTTTCATTGTTCCGCAGGGACAACAATACACCTGACAGGGAAGCAATGTCATTACCGGAACGAATAAGGGCATAGTAGTCAGCTACCGCCGCCGACCGCGTAAAGCGCGGATCGATCACAATCAACTTCCCGCCGTTGTGTATTTTGGCATCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTACGCTGCGTTACCGcc               >  2:110561/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gcacgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGccc              >  1:88052/1‑251 (MQ=255)
                                                                                                                                                                                                                                    agagcgtgCTGGCGGTGGGATAGGCGTCGGTGTATTCGCGGTTGAATTTTTAATTGTTCAGCAGGTACAGCAATACGCCGGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATATTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATc           <  2:148807/244‑1 (MQ=255)
                                                                                                                                                                                                                                             gggCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCGGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCacgacg  <  2:88052/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
GGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACG  >  NC_000913/4084661‑4085148

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