Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,369,920 C→T A165T (GCG→ACG)  yjeH ← L‑methionine/branched chain amino acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,369,9200CT92.5% 106.9 / NA 42A165T (GCG→ACG) yjeHL‑methionine/branched chain amino acid exporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (21/16);  minor base G (1/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGA  >  NC_000913/4369678‑4370153
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
aaaggtgtcaaattgtaaaatgacTACCGTACATCACGAGTAGACAAATCCTGCCAGCAGCAGACCAATCATAAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACTCCAGAACATCAATGATAACGCAGAAAATAACCCGGTAAGTCCGATATTACCAGCTGCCGGAAAGGGGTTATTCGCAGTTTTTATATCGCCCGCCCCCCAgc                                                                                                                                                                                                                                   <  2:273274/230‑2 (MQ=255)
aTAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCAGCCAGCAGCAGACCACTCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAACTTCCGAGGCGAGATGGGCAAATGACTCCAGACCGACAAAACACCAGAACCGACCTGATAACGCAGGAAATAACCCGGTAAGATGGATATTAACAGGTGCCGGCAAGGGGATATTCGACGGTTTGATATCGACCGTccacaagc                                                                                                                                                                                                                                   >  2:166582/1‑247 (MQ=255)
  aGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATa                                                                                                                                                                                                                                 <  1:334636/251‑1 (MQ=255)
                         aCCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATa                                                                                                                                                                                                            >  2:185374/1‑249 (MQ=255)
                                                              gACCAATCAACAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTg                                                                                                                                                                     <  2:117233/251‑1 (MQ=255)
                                                              gACCAATCAACAAAGCACGAGGAAAAGCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACCCCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTg                                                                                                                                                                     <  2:380470/251‑1 (MQ=255)
                                                                 cAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCa                                                                                                                                                                  <  1:231708/251‑1 (MQ=255)
                                                                   atcatcAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGCCGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCAcc                                                                                                                                                                >  2:211364/1‑251 (MQ=255)
                                                                     catcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACcgc                                                                                                                                                              >  1:283837/1‑251 (MQ=255)
                                                                            gCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTAcc                                                                                                                                                       <  1:1089/251‑1 (MQ=255)
                                                                               cGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATAGCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGat                                                                                                                                                    <  1:40708/251‑1 (MQ=255)
                                                                                 gggAAAATAGCGCTCTGGAATTTTAAATTACGCCGCGCGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATGACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGGCCACCAGATAGCGACAATCAGCGCGACGAAAAGTCCGGCAATACCGGTATGTAGCTTAGCACGGGAACAAGCACCGCGAggaacgattt                                                                                                                                                  >  2:70869/2‑242 (MQ=255)
                                                                                 aGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGatat                                                                                                                                                  >  1:59487/1‑251 (MQ=255)
                                                                                         cACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAcc                                                                                                                                          >  1:6610/1‑251 (MQ=255)
                                                                                            gCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAccagc                                                                                                                                       >  1:379844/1‑251 (MQ=255)
                                                                                                             tCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACCCCAGCGCCAGCGTACCGAGTTc                                                                                                                      >  2:191650/1‑251 (MQ=255)
                                                                                                                  ggCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCCGAACATCGCTGATAACGCAGCACATAACCCGGTAAGTTCGATCTTACCAGGTGCCGGAAAGGGGCTATTCGCAGGTTTGATATCACCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTg                                                                                                                                                                           >  2:172083/1‑193 (MQ=255)
                                                                                                                  ggCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTg                                                                                                                                                                           <  1:172083/193‑1 (MQ=255)
                                                                                                                     gagaTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCcaac                                                                                                               >  1:135494/1‑250 (MQ=255)
                                                                                                                             cAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATCACCCGGTAAGTTCGATATTACCAGGTGCCGGCAAGGGGATATTCGCAGGTTTGCTATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATCACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCCGCGTACCGGGTTCTGCCAACAACAGTTGc                                                                                                      >  2:49064/1‑251 (MQ=255)
                                                                                                                              aaaTGCCTCCAGACCGACAAAACACCAGAACATCACTGCTAACGCAGCAAATAACCCGGTAAGATCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCACGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTTGAACTGGAACCGCGAGTACCGCCATACAACACCAGCGCCAGCGTACCACGTTCAGCCACCAACAGATGcc                                                                                                     >  2:214884/1‑251 (MQ=255)
                                                                                                                                             gACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCaa                                                                                      <  1:66221/251‑1 (MQ=255)
                                                                                                                                               cAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATACCGGTTTGTAGATTGGCACTGGAACTGGCACCGCGCGTACCGATATACCACACCAGCGCCAGCGt                                                                                                                               >  2:242148/1‑208 (MQ=255)
                                                                                                                                               cAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGt                                                                                                                               <  1:242148/208‑1 (MQ=255)
                                                                                                                                                aaaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGTGGATATTCGCAGGTGTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATa                                                                                                                                                                                                <  1:233711/142‑1 (MQ=255)
                                                                                                                                                aaaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGGTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATa                                                                                                                                                                                                >  2:233711/1‑142 (MQ=255)
                                                                                                                                                             aTCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGa                                                                       >  1:22967/1‑250 (MQ=255)
                                                                                                                                                                aCTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAAccc                                                                   <  1:185374/251‑1 (MQ=255)
                                                                                                                                                                  tGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAACGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCcgg                                                                 >  1:127190/1‑251 (MQ=255)
                                                                                                                                                                         gcagcaAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAAt                                                          >  1:395609/1‑251 (MQ=255)
                                                                                                                                                                          cagcaAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAAtt                                                         >  1:136080/1‑251 (MQ=255)
                                                                                                                                                                               aaTAACCCGGTAAGTGCGATACGACCAGGTGCCGGAAAGGGGATATTCGCATGTTTGATATCGCCCGTCCCCCAGATAGCGTCAATCAGCTCGACGATAAGTCCGGCAATAACTGTTTTTAGATTAGCACTGGAACTGGCTCCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTGGCGAGCTATGCAAGACAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAg                                                    <  2:18636/251‑1 (MQ=255)
                                                                                                                                                                                  aaCCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGTTAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGAGTAGCACTGGAACTGTCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCTATTTCTGCCTACAACATTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGc                                                 <  1:188482/251‑1 (MQ=255)
                                                                                                                                                                                  aaCCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTgg                                                 >  1:58335/1‑250 (MQ=255)
                                                                                                                                                                                   aCCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCg                                                <  2:146343/251‑1 (MQ=255)
                                                                                                                                                                                         tAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAg                                                    >  1:38850/1‑241 (MQ=255)
                                                                                                                                                                                         tAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAg                                                    <  2:38850/241‑1 (MQ=255)
                                                                                                                                                                                          aaGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGGCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCAGAACCCGGCGGCAATTTGTAGTGCggcaggc                                          <  1:191650/250‑1 (MQ=255)
                                                                                                                                                                                             ttCGATATTACCAGGTGCCGGACAGGGGATATTCGCAGGTTTGATATCTCCCGTCCACCAGATAGCGCCAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTACCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAAc                                      >  2:296219/1‑251 (MQ=255)
                                                                                                                                                                                                     tACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGGCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                               <  2:22967/250‑1 (MQ=255)
                                                                                                                                                                                                      aCCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGa                              >  1:40769/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aGGTTTGATATCGCCCGTCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACCCCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGa  <  1:305854/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
ATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGA  >  NC_000913/4369678‑4370153

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