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,9030TG91.7% 182.8 / ‑2.2 60N307T (AAC→ACC) fdoGformate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (27/28);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACG  >  NC_000913/4084661‑4085145
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
ggAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCCCGATCAGGCTGGCGGTgggct                                                                                                                                                                                                                                               >  1:463621/1‑246 (MQ=255)
 gAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAgg                                                                                                                                                                                                                                            <  1:335271/251‑1 (MQ=255)
   aCGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTGCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCt                                                                                                                                                                                                                                          <  1:15775/251‑1 (MQ=255)
    cgcgTCTTTTGGCGTTCCACAGATGTTTTCATCCACATCTGGCGTGTAGCGGTAAACGTGCTGGTTCAGTAATTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCTTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTACCGGACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACTGGCCATCTTCAAAGCCGTAATCCCCACGTACGATCAGGCTGGCTGTGGTATAGGCtt                                                                                                                                                                                                                                         <  1:103512/251‑1 (MQ=255)
                 gTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATTGGCTTTGGTGTATTcgcg                                                                                                                                                                                                                            >  2:42102/1‑251 (MQ=255)
                       caGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTAGGTGTAtttccgtgtga                                                                                                                                                                                                                       >  2:177397/1‑241 (MQ=255)
                         gATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAAtt                                                                                                                                                                                                                    <  1:240436/251‑1 (MQ=255)
                               ttCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTctt                                                                                                                                                                                                              >  1:354833/1‑249 (MQ=255)
                                         tCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTcagc                                                                                                                                                                                                     >  1:402109/1‑249 (MQ=255)
                                          cTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAgg                                                                                                                                                                                                  >  1:145176/1‑251 (MQ=255)
                                                   aGCGGGAAACGTGCCGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAAt                                                                                                                                                                                         <  1:64956/251‑1 (MQ=255)
                                                    gCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATa                                                                                                                                                                                        <  2:132796/251‑1 (MQ=255)
                                                          aaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGAGTTTTTCCGCGTCGTAGCCGGTGACCAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGGGGAGGTAAGGGCTTCGGGGTATTCGCAGTTGAATTTTTCATTGGTCAGAAgggccccaaaatcgcctg                                                                                                                                                                                  >  2:259985/1‑235 (MQ=255)
                                                              tGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGGATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTAAAGCAATACCCCTGAc                                                                                                                                                                                >  1:446348/1‑249 (MQ=255)
                                                                 tGTTTCAGCAAGTTCCACACGCCGCGCGGATGTTGCAGCGTGGTTTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAGGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTTGTAGCCGGTGAACAGGCAATCTTCACAGCCGTAATCCTCACGCCCGATAAGGCTGGCGGTGGTTAAGGCTTCGGTGTTTTCGCGGTGGCATTTTTCCTTGTTCAGCAGGTAccgccttaagcccgacagg                                                                                                                                                                             >  2:18084/1‑231 (MQ=255)
                                                                      cAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCaa                                                                                                                                                                      <  1:13605/251‑1 (MQ=255)
                                                                        gCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATg                                                                                                                                                                    <  1:410542/251‑1 (MQ=255)
                                                                           aGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGCTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGGAAGCAATGTc                                                                                                                                                                  >  1:301464/1‑250 (MQ=255)
                                                                              tCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTAt                                                                                                                                                                                                                                 <  1:378047/185‑1 (MQ=255)
                                                                               ccACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGATTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATAAGCCTGACAGGAAAGAAATGTCAgat                                                                                                                                                              >  1:84374/1‑248 (MQ=255)
                                                                                  cacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGc                                                                                                                                                                                     <  1:416810/224‑1 (MQ=255)
                                                                                  cacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGGACAGCAATAccc                                                                                                                                                                                     >  2:416810/1‑222 (MQ=255)
                                                                                                          ggTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACCGGAAAGCAATGTCAGTACCGGCACGAATAGGGGCAGAGGAGTc                                                                                                                                    >  2:298834/1‑249 (MQ=255)
                                                                                                           gTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGc                                                                                                                                 >  1:388363/1‑251 (MQ=255)
                                                                                                             aTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTACTAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCACGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGGGACATCGCGGTTGAATTTTTCATTGTTCAGCCGGTACAGCCATACGCCTGACCGGAGAGCAATGTCAGTACCGGGACGGATAGGGGCATAGTAGTCAGCt                                                                                                                                >  1:427609/1‑250 (MQ=255)
                                                                                                                  gTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgcc                                                                                                                          <  1:102889/251‑1 (MQ=255)
                                                                                                                      ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTGGCGTTTTGCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccgcc                                                                                                                       <  1:42102/250‑1 (MQ=255)
                                                                                                                           agacGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcgc                                                                                                                 <  1:164774/248‑1 (MQ=255)
                                                                                                                               gTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCGCGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTAAAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAAGCAGCTACCGCCGCCGTACGCGTa                                                                                                              >  1:73056/1‑250 (MQ=255)
                                                                                                                                     gaccAGTTCAAAAGGCAAGGATGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGGGAACAGGCCATCTTAAAAGGCGTAATCCTCACGCACGATCAGGCTGGGGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAAGGTCAGTACCGGAACGAAGAGGGGAATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAg                                                                                                       <  2:189465/249‑1 (MQ=255)
                                                                                                                                       ccAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGga                                                                                                     <  1:318006/251‑1 (MQ=255)
                                                                                                                                        cAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                    <  1:469329/251‑1 (MQ=255)
                                                                                                                                          gtcaaTAATTCCAGGAGGATTGATCGTATTGGTGTTTTTCCGCGTCGTAGCCGGGGGACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTACGGTGTAGTCGCGGTTGAATTTTTCATTTTTCAGCAGGGACAGCTATACGCCTGACAGGAAAGCAATGTCAGGACCGGAACGAATAGGGGAAAAGTAGTCAGCTCCCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                    <  2:231013/245‑1 (MQ=255)
                                                                                                                                                  gTCCAGGAGGATTGATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAAt                                                                                          <  1:160114/251‑1 (MQ=255)
                                                                                                                                                      aggaggATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGAGCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCAACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGc                                                                                      <  2:367703/251‑1 (MQ=255)
                                                                                                                                                               tATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGGTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCATGACAGGAAAACACTGTCAGTGCCGGAAAGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCg                                                                             >  1:213714/1‑251 (MQ=255)
                                                                                                                                                                  cgtcctggCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGtt                                                                           <  1:228288/243‑1 (MQ=255)
                                                                                                                                                                  cGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGTCTCGGGGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGACTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGGAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTtg                                                                          >  1:5896/1‑251 (MQ=255)
                                                                                                                                                                     aCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCTTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCATTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGa                                                                       <  1:71632/251‑1 (MQ=255)
                                                                                                                                                                     aCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAACTACCGCCGCCGTACGCGTAAAGCGAGGATCGATAACAAACAGCTTCGCGCCGTtgtg                                                                        >  2:142031/1‑250 (MQ=255)
                                                                                                                                                                        tGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGCATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCGGTCCCGGAACGAATAGGGGCATAGTAGTCCGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGCTCACAATCAGCTTCGCGCCGTTGTGAAtt                                                                    >  1:283594/1‑251 (MQ=255)
                                                                                                                                                                         gCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGGACAGCAATACGCCTGACAGGAAAGCAATTTCAATACCGGAACGAATAGGGGCATAGTAGGCAGCTACCGCCGCCGTAAGCGTAAAACGGGGATCGATCCAAATCAGCTTCGCGCCGTTGTGAAttt                                                                   >  2:322530/1‑251 (MQ=255)
                                                                                                                                                                             ttttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCTCCGTTGTGAATTTTgg                                                                <  1:419946/250‑1 (MQ=255)
                                                                                                                                                                             gtttCCGCGTCGTAGCCTTTGAACAGGCCTCTTTCAAAGCCGTAATCCTCACGCTCGATCAGTCTTGCGTTGGTATAGGCTTCGGTGTATTCGCGTTTGAATTTTTCATTGTTTAGCAGCTACATCAATACGCCTGACAGAAAAGCAATTTCAGTACCGGAACGAATAGGGGCATATTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGATGATCGATCACAATCAGCTTCGCGCCTTTGTGAATTTTgg                                                                <  1:419954/249‑1 (MQ=255)
                                                                                                                                                                                    cgtcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCAGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTAAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGGCAGTACCGGAACGGATAGGGGCATAGGAGTCAGCTACCGCCGCCGTAAGCGTAAAGCGAGGAGCGATCACAATCAGCTTCGCGACGTTGTGAATTTTGGCTTCCATc                                                        >  1:435354/1‑251 (MQ=255)
                                                                                                                                                                                      tcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGACTAGGGGCATAGTAGTCAGCTACCGCCGCCGGACGCGGAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGc                                                      >  2:274436/1‑251 (MQ=255)
                                                                                                                                                                                         tAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCCAAGTCAGTACCGGAACGGATAAGGGCATAGTAGTCAGCTAACGCCGCCGTACGCGTAAAGCGAGGGTCGGTCACAAGCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATAGcccc                                                   >  2:295999/1‑250 (MQ=255)
                                                                                                                                                                                         tAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGAATAGGCTTCGGTCTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAACTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGcccc                                                   >  1:66183/1‑250 (MQ=255)
                                                                                                                                                                                         gaGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCa                                                   <  1:309768/250‑1 (MQ=255)
                                                                                                                                                                                          agccgtggAACAGGCCATTGTCAAAGCCGTAATCCTCACGCACGATCTGGCTGGCGGTGGTATAGGCTTCGGTTTATTCGCGGTTGAATTTTTCATGTTCCAGCAGGAACAGCAATACGCCTGACAGGAAATCAATGTCAGTACCGGATCGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCa                                                   <  2:127483/243‑1 (MQ=255)
                                                                                                                                                                                              ggTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCCTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCCGTACCGGAACGAATAGGGGCATAGTAGGTAACTACCGCCGCCGCAAGCGGAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGGTGTGAATTTTGGCTTTCATCGCCCAGCgg                                               >  2:172652/1‑250 (MQ=255)
                                                                                                                                                                                                ggAACATGCCATCTTCAAAGCCGTAATCCTCCCGCTGGATCAGGCTGGCGGTGGTAGTGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATCCGCCTGACAGGAAAGCAATGTAAGTACCGGAACTAATAGGGGCATATTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAc                                            <  2:378849/250‑1 (MQ=255)
                                                                                                                                                                                                 gAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAc                                            <  1:363485/250‑1 (MQ=255)
                                                                                                                                                                                                                       gTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCTCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcgg                      >  1:316395/1‑251 (MQ=255)
                                                                                                                                                                                                                        tAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAAAGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGACCAGCGGAACCCGACCAGGTGAGCTTCAgcggg                     >  1:225012/1‑250 (MQ=255)
                                                                                                                                                                                                                         aacccTCCCGCACGATCTGGCGGGCGGTGGTATAGGCTGCGGTGTATTCGCGGTTGAATTTTTCATTTTTCAGCAGGTACAGAAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACAGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcggc                     <  2:354833/247‑1 (MQ=255)
                                                                                                                                                                                                                         aaTCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATTTCAGTACCGGAAAGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGAACGCGTAAAGCGAGGAGAGATACCAAACAGATTCGCGCCGTTGGGAATTTTGACTTCCATCGCCCAGCGGAACCCGCCCGGGGGGGCTTCaacagtg                    >  2:349287/1‑245 (MQ=255)
                                                                                                                                                                                                                          aTCCTCACTCACTTTCGGGCTGGCGTGTGTATAGGTTCCTGTGTATTCTCGGTTGAATTTTTCATTGTCCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTCCGGGAACGAATAGGGGCATAGTAGTCATCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGt                   <  2:388363/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gATCAGTCTGGCGGTTGAATAGGCTTCGGTGAATTCGCGGTTGAATTTGTCATTATTCAGCAGGTACAGCAATCCGCCTGACAGGAAAGCAATGTCAGTTCCGGAACGAATAGGGGCATAGTAGTCAGCTCCCTCCTCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATca       <  1:133123/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGAATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCacg  <  1:295999/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
GGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACG  >  NC_000913/4084661‑4085145

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