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,9030TG95.5% 144.7 / ‑4.5 44N307T (AAC→ACC) fdoGformate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (20/22);  total (21/23)
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.25e-01

TTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCA  >  NC_000913/4084657‑4085140
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
ttCAGGAACGCGTCTTTTGGCTGACCACAGATGTGTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTTGTTCAGCAAGTCCCACACGCAGCGCGTATGGTTCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTgg                                                                                                                                                                                                                                             <  2:45403/250‑1 (MQ=255)
     gAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAgg                                                                                                                                                                                                                                       <  1:126903/251‑1 (MQ=255)
               tttGGCGTACCTCAGATGTTTTCAACCACATCTGGCGTGTAGCGGTAAACGTGCTGTTTCATCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCTCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCTCGCACGATCAGTCTGGCGGTGGTATAGGCTTCGgtg                                                                                                                                                                                                                               <  1:154065/249‑1 (MQ=255)
                      tACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGCTGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCgg                                                                                                                                                                                                                      >  1:411684/1‑251 (MQ=255)
                        ccACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGtt                                                                                                                                                                                                                    >  1:110165/1‑251 (MQ=255)
                             gATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTGCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCTTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTAATCGCGGTTGAAtt                                                                                                                                                                                                               <  1:288801/251‑1 (MQ=255)
                                  tttCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑c                                                                                                                                                                                                          <  1:404650/250‑1 (MQ=255)
                                        ccACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTc                                                                                                                                                                                                   <  1:32835/251‑1 (MQ=255)
                                           caTCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGGTGGATTTTT‑TATTGTTcagc                                                                                                                                                                                                >  2:264642/1‑251 (MQ=255)
                                                ggCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAgg                                                                                                                                                                                             <  2:347627/249‑1 (MQ=255)
                                                 gCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTGCAGCAGGTAc                                                                                                                                                                                          <  2:63838/251‑1 (MQ=255)
                                                   gggtAGCGGGAAACGGGGTGTTTCAGCAAGTGCCACCCTCAGCCCGGATGGTGCAGCGTGGTATCGCGTTTGGCTAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGGAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCTCGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGTTACAg                                                                                                                                                                                        <  2:229031/249‑1 (MQ=255)
                                                               aCGTGCTGTTTCAGCAAGTGCCACACGCCGCGCGGATGGTGCAGCGTGGTATCGCGGTTGTCGAAGCCGTTTTCGTCCAGTTAATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAATAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCTCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTgc                                                                                                                                                                            <  2:123243/251‑2 (MQ=255)
                                                                                   ccACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAc                                                                                                                                                        <  1:264642/251‑1 (MQ=255)
                                                                                           aGCGCGGATGTTGCAGCGTGGTATCGGGTTTGGCGAAGCCGTTTTCGTCCAGTGCATAAGTCCAGGAGGATTTATCGTACTTTCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCTGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGa                                                                                                                                                <  1:407760/251‑1 (MQ=255)
                                                                                                       gCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCCTCGGGGTATTCGCGGTTGAATTTTTTAATTGT‑CCGCAGGTTCAGCAATACGCCTGGAACGAAAGAAATGTccattccggaccgaaaaagggcatg                                                                                                                                     >  2:190211/1‑226 (MQ=255)
                                                                                                           cttGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATagtagt                                                                                                                                >  1:399735/3‑251 (MQ=255)
                                                                                                           cGTGGAATTGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGTATGATTTATCTTACTTTCTTTTTACCGCGGCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCTGTGTATTCTCGGTGGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACTCCTGACAGGAAAGCAATGGCAGTACCGGAACGAATAGGGGCATtgtagt                                                                                                                                <  2:58644/251‑1 (MQ=255)
                                                                                                                        ttGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTTCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCGGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccg                                                                                                                    <  1:43256/250‑1 (MQ=255)
                                                                                                                             gAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTTACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAACGCCGCCGTAc                                                                                                              >  1:176303/1‑251 (MQ=255)
                                                                                                                              aaGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTAGCGTACTTTCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcg                                                                                                             <  1:38979/251‑1 (MQ=255)
                                                                                                                                                 cATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatcgatc                                                                                          >  1:76669/1‑251 (MQ=255)
                                                                                                                                                 cATAAGGCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatcgat                                                                                           <  1:270198/250‑1 (MQ=255)
                                                                                                                                                          aggaggATTTATCGTACTTGCGGTTTTCCGCGACGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGTGTTTTT‑AATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGAACCGGAACGAATAAGGGCATAGTAGTCAGCTACCGCCGCCGTACGAGTAAAGCGCGGATCGATCACAATCAGc                                                                                 >  1:8186/1‑251 (MQ=255)
                                                                                                                                                             aggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTc                                                                              <  1:1831/251‑1 (MQ=255)
                                                                                                                                                              ggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcg                                                                             >  1:310272/1‑251 (MQ=255)
                                                                                                                                                                    ctCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTACAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                              >  1:206282/2‑148 (MQ=255)
                                                                                                                                                                    aTCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTACAGGCTGCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                              <  2:206282/148‑1 (MQ=255)
                                                                                                                                                                      cGTACTAGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGACCGAATAGGCGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGtt                                                                      <  1:198769/250‑1 (MQ=255)
                                                                                                                                                                               gTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTg                                                            >  1:189590/1‑251 (MQ=255)
                                                                                                                                                                                     ccGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGGTCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTcc                                                      >  2:202309/1‑251 (MQ=255)
                                                                                                                                                                                     ccGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCCGGTACAGCAATACGCCTGGCAGGAAAGCACTGTCCGTACCGGAACGAATAGGGGCATAGGAGTCAGCTACCGCCGCCGTCCGCGTAAAGCGAGGATCGGTCACAATCAGCTTCGCCCCGTTGTGAATTTTGGCTTcc                                                      >  2:384251/1‑251 (MQ=255)
                                                                                                                                                                                      cgcgTCGGACCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGTGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCa                                                     <  2:289260/251‑1 (MQ=255)
                                                                                                                                                                                           cgtATCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATCAGCCTCGGTGTATTCGCGGTTGAATTTTT‑CATTGGTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAAGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGGGGATCGATCACAATAAGCTTCGCGCCGTTTTGACTATTGGCTTCCATCGcc                                                >  2:49932/1‑251 (MQ=255)
                                                                                                                                                                                            gtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGccc                                               >  1:104555/1‑251 (MQ=255)
                                                                                                                                                                                                    tGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                           <  1:297612/215‑1 (MQ=255)
                                                                                                                                                                                                    tGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGGACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                           >  2:297612/1‑215 (MQ=255)
                                                                                                                                                                                                         aGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAc                                  >  1:11950/1‑251 (MQ=255)
                                                                                                                                                                                                               tCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGCGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGGTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTg                            >  1:157169/1‑251 (MQ=255)
                                                                                                                                                                                                                  tCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAATCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCAGCGCCCAGCGGAACCCGACCGGGTGa                           <  1:342688/249‑1 (MQ=255)
                                                                                                                                                                                                                  ccAAAGCCGTAATCCGCCCGCACGATCAGGCGGGCTTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGa                           <  2:399735/248‑1 (MQ=255)
                                                                                                                                                                                                                                    cgcacgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGcc       >  1:146647/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gATCAGGCTGTCGGCGGTATAGGCTTCGGTGTATTCGCGGTTGCATTTTT‑CATTGTTCAGCAGGTACAGCCATACGCCTGACAGGAAAGCCATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGCACGCGTAAAGCGAGGATCGGTCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCCTCGCCCAGCGGAACCCGACCGGGGGCGCTTCAGCGGCGTTACCGcccctcc  >  2:58324/1‑247 (MQ=255)
                                                                                                                                                                                                                                         gATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGTCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATc   >  1:397126/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
TTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTT‑CATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCA  >  NC_000913/4084657‑4085140

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