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

GGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACG  >  NC_000913/4084661‑4085145
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
ggAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCCCGATCAGGCTGGCGGTgggct                                                                                                                                                                                                                                               >  1:463621/1‑246 (MQ=255)
 gAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAgg                                                                                                                                                                                                                                            <  1:335271/251‑1 (MQ=255)
   aCGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTGCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCt                                                                                                                                                                                                                                          <  1:15775/251‑1 (MQ=255)
    cgcgTCTTTTGGCGTTCCACAGATGTTTTCATCCACATCTGGCGTGTAGCGGTAAACGTGCTGGTTCAGTAATTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCTTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTACCGGACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACTGGCCATCTT‑CAAAGCCGTAATCCCCACGTACGATCAGGCTGGCTGTGGTATAGGCtt                                                                                                                                                                                                                                         <  1:103512/251‑1 (MQ=255)
                 gTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATTGGCTTTGGTGTATTcgcg                                                                                                                                                                                                                            >  2:42102/1‑251 (MQ=255)
                       caGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTAGGTGTAtttccgtgtga                                                                                                                                                                                                                       >  2:177397/1‑241 (MQ=255)
                         gATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAAtt                                                                                                                                                                                                                    <  1:240436/251‑1 (MQ=255)
                               ttCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTctt                                                                                                                                                                                                              >  1:354833/1‑249 (MQ=255)
                                         tCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTc‑agc                                                                                                                                                                                                     >  1:402109/1‑249 (MQ=255)
                                          cTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAgg                                                                                                                                                                                                  >  1:145176/1‑251 (MQ=255)
                                                   aGCGGGAAACGTGCCGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAAt                                                                                                                                                                                         <  1:64956/251‑1 (MQ=255)
                                                    gCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATa                                                                                                                                                                                        <  2:132796/251‑1 (MQ=255)
                                                          aaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGAGTTTTTCCGCGTCGTAGCCGGTGACCAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGGGGAGGTAAGGGCTTCGGGGTATTCGCAGTTGAATTTTTCATTGGTC‑AGAAgggccccaaaatcgcctg                                                                                                                                                                                  >  2:259985/1‑235 (MQ=255)
                                                              tGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGGATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTAAAGCAATACCCCTGAc                                                                                                                                                                                >  1:446348/1‑249 (MQ=255)
                                                                 tGTTTCAGCAAGTTCCACACGCCGCGCGGATGTTGCAGCGTGGTTTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAGGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTTGTAGCCGGTGAACAGGCAATCTT‑CACAGCCGTAATCCTCACGCCCGATAAGGCTGGCGGTGGTTAAGGCTTCGGTGTTTTCGCGGTGGCATTTTTCCTTGTTC‑AGCAGGTAccgccttaagcccgacagg                                                                                                                                                                             >  2:18084/1‑231 (MQ=255)
                                                                      cAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCaa                                                                                                                                                                      <  1:13605/251‑1 (MQ=255)
                                                                        gCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATg                                                                                                                                                                    <  1:410542/251‑1 (MQ=255)
                                                                           aGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGCTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGGAAGCAATGTc                                                                                                                                                                  >  1:301464/1‑250 (MQ=255)
                                                                              tCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTAt                                                                                                                                                                                                                                 <  1:378047/185‑1 (MQ=255)
                                                                               ccACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGATTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATAAGCCTGACAGGAAAGAAATGTCAgat                                                                                                                                                              >  1:84374/1‑248 (MQ=255)
                                                                                  cacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGc                                                                                                                                                                                     <  1:416810/224‑1 (MQ=255)
                                                                                  cacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGGACAGCAATAccc                                                                                                                                                                                     >  2:416810/1‑222 (MQ=255)
                                                                                                          ggTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACCGGAAAGCAATGTCAGTACCGGCACGAATAGGGGCAGAGGAGTc                                                                                                                                    >  2:298834/1‑249 (MQ=255)
                                                                                                           gTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGc                                                                                                                                 >  1:388363/1‑251 (MQ=255)
                                                                                                             aTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTACTAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCACGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGGGACATCGCGGTTGAATTTTTCATTGTTC‑AGCCGGTACAGCCATACGCCTGACCGGAGAGCAATGTCAGTACCGGGACGGATAGGGGCATAGTAGTCAGCt                                                                                                                                >  1:427609/1‑250 (MQ=255)
                                                                                                                  gTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgcc                                                                                                                          <  1:102889/251‑1 (MQ=255)
                                                                                                                      ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTGGCGTTTTGCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccgcc                                                                                                                       <  1:42102/250‑1 (MQ=255)
                                                                                                                           agacGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcgc                                                                                                                 <  1:164774/248‑1 (MQ=255)
                                                                                                                               gTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCGCGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTAAAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAAGCAGCTACCGCCGCCGTACGCGTa                                                                                                              >  1:73056/1‑250 (MQ=255)
                                                                                                                                     gaccAGTTCAAAAGGCAAGGATGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGGGAACAGGCCATCTT‑AAAAGGCGTAATCCTCACGCACGATCAGGCTGGGGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAAGGTCAGTACCGGAACGAAGAGGGGAATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAg                                                                                                       <  2:189465/249‑1 (MQ=255)
                                                                                                                                       ccAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGga                                                                                                     <  1:318006/251‑1 (MQ=255)
                                                                                                                                        cAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                    <  1:469329/251‑1 (MQ=255)
                                                                                                                                          gtcaaTAATTCCAGGAGGATTGATCGTATTGGTGTTTTTCCGCGTCGTAGCCGGGGGACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTACGGTGTAGTCGCGGTTGAATTTTTCATTTTTC‑AGCAGGGACAGCTATACGCCTGACAGGAAAGCAATGTCAGGACCGGAACGAATAGGGGAAAAGTAGTCAGCTCCCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                    <  2:231013/245‑1 (MQ=255)
                                                                                                                                                  gTCCAGGAGGATTGATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAAt                                                                                          <  1:160114/251‑1 (MQ=255)
                                                                                                                                                      aggaggATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGAGCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCAACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGc                                                                                      <  2:367703/251‑1 (MQ=255)
                                                                                                                                                               tATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGGTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCATGACAGGAAAACACTGTCAGTGCCGGAAAGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCg                                                                             >  1:213714/1‑251 (MQ=255)
                                                                                                                                                                  cgtcctggCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGtt                                                                           <  1:228288/243‑1 (MQ=255)
                                                                                                                                                                  cGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGTCTCGGGGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGACTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGGAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTtg                                                                          >  1:5896/1‑251 (MQ=255)
                                                                                                                                                                     aCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCTTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCATTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGa                                                                       <  1:71632/251‑1 (MQ=255)
                                                                                                                                                                     aCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAACTACCGCCGCCGTACGCGTAAAGCGAGGATCGATAACAAACAGCTTCGCGCCGTtgtg                                                                        >  2:142031/1‑250 (MQ=255)
                                                                                                                                                                        tGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGCATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCGGTCCCGGAACGAATAGGGGCATAGTAGTCCGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGCTCACAATCAGCTTCGCGCCGTTGTGAAtt                                                                    >  1:283594/1‑251 (MQ=255)
                                                                                                                                                                         gCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGGACAGCAATACGCCTGACAGGAAAGCAATTTCAATACCGGAACGAATAGGGGCATAGTAGGCAGCTACCGCCGCCGTAAGCGTAAAACGGGGATCGATCCAAATCAGCTTCGCGCCGTTGTGAAttt                                                                   >  2:322530/1‑251 (MQ=255)
                                                                                                                                                                             ttttCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCTCCGTTGTGAATTTTgg                                                                <  1:419946/250‑1 (MQ=255)
                                                                                                                                                                             gtttCCGCGTCGTAGCCTTTGAACAGGCC‑TCTTTCAAAGCCGTAATCCTCACGCTCGATCAGTCTTGCGTTGGTATAGGCTTCGGTGTATTCGCGTTTGAATTTTTCATTGTTT‑AGCAGCTACATCAATACGCCTGACAGAAAAGCAATTTCAGTACCGGAACGAATAGGGGCATATTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGATGATCGATCACAATCAGCTTCGCGCCTTTGTGAATTTTgg                                                                <  1:419954/249‑1 (MQ=255)
                                                                                                                                                                                    cgtcgtAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCAGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTA‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGGCAGTACCGGAACGGATAGGGGCATAGGAGTCAGCTACCGCCGCCGTAAGCGTAAAGCGAGGAGCGATCACAATCAGCTTCGCGACGTTGTGAATTTTGGCTTCCATc                                                        >  1:435354/1‑251 (MQ=255)
                                                                                                                                                                                      tcgtAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGACTAGGGGCATAGTAGTCAGCTACCGCCGCCGGACGCGGAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGc                                                      >  2:274436/1‑251 (MQ=255)
                                                                                                                                                                                         tAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCCAAGTCAGTACCGGAACGGATAAGGGCATAGTAGTCAGCTAACGCCGCCGTACGCGTAAAGCGAGGGTCGGTCACAAGCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATAGcccc                                                   >  2:295999/1‑250 (MQ=255)
                                                                                                                                                                                         tAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGAATAGGCTTCGGTCTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAACTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGcccc                                                   >  1:66183/1‑250 (MQ=255)
                                                                                                                                                                                         gaGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCa                                                   <  1:309768/250‑1 (MQ=255)
                                                                                                                                                                                          agccgtggAACAGGCCATTGT‑CAAAGCCGTAATCCTCACGCACGATCTGGCTGGCGGTGGTATAGGCTTCGGTTTATTCGCGGTTGAATTTTTCAT‑GTTCCAGCAGGAACAGCAATACGCCTGACAGGAAATCAATGTCAGTACCGGATCGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCa                                                   <  2:127483/243‑1 (MQ=255)
                                                                                                                                                                                              ggTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCCTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCCGTACCGGAACGAATAGGGGCATAGTAGGTAACTACCGCCGCCGCAAGCGGAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGGTGTGAATTTTGGCTTTCATCGCCCAGCgg                                               >  2:172652/1‑250 (MQ=255)
                                                                                                                                                                                                ggAACATGCCATCTT‑CAAAGCCGTAATCCTCCCGCTGGATCAGGCTGGCGGTGGTAGTGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATCCGCCTGACAGGAAAGCAATGTAAGTACCGGAACTAATAGGGGCATATTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAc                                            <  2:378849/250‑1 (MQ=255)
                                                                                                                                                                                                 gAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAc                                            <  1:363485/250‑1 (MQ=255)
                                                                                                                                                                                                                       gTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCTCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcgg                      >  1:316395/1‑251 (MQ=255)
                                                                                                                                                                                                                        tAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAAAGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGACCAGCGGAACCCGACCAGGTGAGCTTCAgcggg                     >  1:225012/1‑250 (MQ=255)
                                                                                                                                                                                                                         aacccTCCCGCACGATCTGGCGGGCGGTGGTATAGGCTGCGGTGTATTCGCGGTTGAATTTTTCATTTTTC‑AGCAGGTACAGAAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACAGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcggc                     <  2:354833/247‑1 (MQ=255)
                                                                                                                                                                                                                         aaTCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATTTCAGTACCGGAAAGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGAACGCGTAAAGCGAGGAGAGATACCAAACAGATTCGCGCCGTTGGGAATTTTGACTTCCATCGCCCAGCGGAACCCGCCCGGGGGGGCTTCaacagtg                    >  2:349287/1‑245 (MQ=255)
                                                                                                                                                                                                                          aTCCTCACTCACTTTCGGGCTGGCGTGTGTATAGGTTCCTGTGTATTCTCGGTTGAATTTTTCATTGTCC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTCCGGGAACGAATAGGGGCATAGTAGTCATCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGt                   <  2:388363/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gATCAGTCTGGCGGTTGAATAGGCTTCGGTGAATTCGCGGTTGAATTTGTCATTATTC‑AGCAGGTACAGCAATCCGCCTGACAGGAAAGCAATGTCAGTTCCGGAACGAATAGGGGCATAGTAGTCAGCTCCCTCCTCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATca       <  1:133123/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGAATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCacg  <  1:295999/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
GGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTT‑CAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTC‑AGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACG  >  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: ‑