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

TTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCA  >  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)
                                  tttCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTc                                                                                                                                                                                                          <  1:404650/250‑1 (MQ=255)
                                        ccACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTc                                                                                                                                                                                                   <  1:32835/251‑1 (MQ=255)
                                           caTCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGGTGGATTTTTTATTGTTcagc                                                                                                                                                                                                >  2:264642/1‑251 (MQ=255)
                                                ggCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAgg                                                                                                                                                                                             <  2:347627/249‑1 (MQ=255)
                                                 gCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTGCAGCAGGTAc                                                                                                                                                                                          <  2:63838/251‑1 (MQ=255)
                                                   gggtAGCGGGAAACGGGGTGTTTCAGCAAGTGCCACCCTCAGCCCGGATGGTGCAGCGTGGTATCGCGTTTGGCTAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGGAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCTCGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGTTACAg                                                                                                                                                                                        <  2:229031/249‑1 (MQ=255)
                                                               aCGTGCTGTTTCAGCAAGTGCCACACGCCGCGCGGATGGTGCAGCGTGGTATCGCGGTTGTCGAAGCCGTTTTCGTCCAGTTAATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAATAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCTCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTgc                                                                                                                                                                            <  2:123243/251‑2 (MQ=255)
                                                                                   ccACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAc                                                                                                                                                        <  1:264642/251‑1 (MQ=255)
                                                                                           aGCGCGGATGTTGCAGCGTGGTATCGGGTTTGGCGAAGCCGTTTTCGTCCAGTGCATAAGTCCAGGAGGATTTATCGTACTTTCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCTGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGa                                                                                                                                                <  1:407760/251‑1 (MQ=255)
                                                                                                       gCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCCTCGGGGTATTCGCGGTTGAATTTTTTAATTGTCCGCAGGTTCAGCAATACGCCTGGAACGAAAGAAATGTccattccggaccgaaaaagggcatg                                                                                                                                     >  2:190211/1‑226 (MQ=255)
                                                                                                           cttGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATagtagt                                                                                                                                >  1:399735/3‑251 (MQ=255)
                                                                                                           cGTGGAATTGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGTATGATTTATCTTACTTTCTTTTTACCGCGGCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCTGTGTATTCTCGGTGGAATTTTTCATTGTTCAGCAGGTACAGCAATACTCCTGACAGGAAAGCAATGGCAGTACCGGAACGAATAGGGGCATtgtagt                                                                                                                                <  2:58644/251‑1 (MQ=255)
                                                                                                                        ttGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTTCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCGGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccg                                                                                                                    <  1:43256/250‑1 (MQ=255)
                                                                                                                             gAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTTACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAACGCCGCCGTAc                                                                                                              >  1:176303/1‑251 (MQ=255)
                                                                                                                              aaGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTAGCGTACTTTCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcg                                                                                                             <  1:38979/251‑1 (MQ=255)
                                                                                                                                                 cATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatcgatc                                                                                          >  1:76669/1‑251 (MQ=255)
                                                                                                                                                 cATAAGGCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatcgat                                                                                           <  1:270198/250‑1 (MQ=255)
                                                                                                                                                          aggaggATTTATCGTACTTGCGGTTTTCCGCGACGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGTGTTTTTAATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGAACCGGAACGAATAAGGGCATAGTAGTCAGCTACCGCCGCCGTACGAGTAAAGCGCGGATCGATCACAATCAGc                                                                                 >  1:8186/1‑251 (MQ=255)
                                                                                                                                                             aggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTc                                                                              <  1:1831/251‑1 (MQ=255)
                                                                                                                                                              ggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcg                                                                             >  1:310272/1‑251 (MQ=255)
                                                                                                                                                                    ctCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTACAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                              >  1:206282/2‑148 (MQ=255)
                                                                                                                                                                    aTCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTACAGGCTGCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                              <  2:206282/148‑1 (MQ=255)
                                                                                                                                                                      cGTACTAGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGACCGAATAGGCGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGtt                                                                      <  1:198769/250‑1 (MQ=255)
                                                                                                                                                                               gTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTg                                                            >  1:189590/1‑251 (MQ=255)
                                                                                                                                                                                     ccGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGGTCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTcc                                                      >  2:202309/1‑251 (MQ=255)
                                                                                                                                                                                     ccGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCCGGTACAGCAATACGCCTGGCAGGAAAGCACTGTCCGTACCGGAACGAATAGGGGCATAGGAGTCAGCTACCGCCGCCGTCCGCGTAAAGCGAGGATCGGTCACAATCAGCTTCGCCCCGTTGTGAATTTTGGCTTcc                                                      >  2:384251/1‑251 (MQ=255)
                                                                                                                                                                                      cgcgTCGGACCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGTGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCa                                                     <  2:289260/251‑1 (MQ=255)
                                                                                                                                                                                           cgtATCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATCAGCCTCGGTGTATTCGCGGTTGAATTTTTCATTGGTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAAGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGGGGATCGATCACAATAAGCTTCGCGCCGTTTTGACTATTGGCTTCCATCGcc                                                >  2:49932/1‑251 (MQ=255)
                                                                                                                                                                                            gtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGccc                                               >  1:104555/1‑251 (MQ=255)
                                                                                                                                                                                                    tGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                           <  1:297612/215‑1 (MQ=255)
                                                                                                                                                                                                    tGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGGACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                           >  2:297612/1‑215 (MQ=255)
                                                                                                                                                                                                         aGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAc                                  >  1:11950/1‑251 (MQ=255)
                                                                                                                                                                                                               tCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGCGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGGTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTg                            >  1:157169/1‑251 (MQ=255)
                                                                                                                                                                                                                  tCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAATCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCAGCGCCCAGCGGAACCCGACCGGGTGa                           <  1:342688/249‑1 (MQ=255)
                                                                                                                                                                                                                  ccAAAGCCGTAATCCGCCCGCACGATCAGGCGGGCTTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGa                           <  2:399735/248‑1 (MQ=255)
                                                                                                                                                                                                                                    cgcacgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGcc       >  1:146647/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gATCAGGCTGTCGGCGGTATAGGCTTCGGTGTATTCGCGGTTGCATTTTTCATTGTTCAGCAGGTACAGCCATACGCCTGACAGGAAAGCCATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGCACGCGTAAAGCGAGGATCGGTCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCCTCGCCCAGCGGAACCCGACCGGGGGCGCTTCAGCGGCGTTACCGcccctcc  >  2:58324/1‑247 (MQ=255)
                                                                                                                                                                                                                                         gATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGTCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATc   >  1:397126/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
TTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCA  >  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: