Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,357,202 C→T intergenic (‑92/+221) sapA ← / ← ymjA putative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2020CT94.3% 161.4 / 0.9 55intergenic (‑92/+221)sapA/ymjAputative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (32/18);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTC  >  NC_000913/1356958‑1357444
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
tttGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGtgg                                                                                                                                                                                                                                                   >  1:167084/1‑251 (MQ=255)
         gTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGtt                                                                                                                                                                                                                                          >  2:144003/1‑251 (MQ=255)
         gTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGtt                                                                                                                                                                                                                                          <  1:144003/251‑1 (MQ=255)
           tAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTg                                                                                                                                                                                                                                         >  2:579931/1‑250 (MQ=255)
                        ttGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGtt                                                                                                                                                                                                                                          <  2:385590/236‑1 (MQ=255)
                        ttGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGtt                                                                                                                                                                                                                                          >  1:385590/1‑236 (MQ=255)
                           cccGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGcc                                                                                                                                                                                                                         >  2:111323/1‑250 (MQ=255)
                                              aaaCCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTTACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTAcc                                                                                                                                                                                                    <  2:188669/251‑1 (MQ=255)
                                                 ccgctggatcgGATATCAGCATGCGGGGGAGATTCAGGCGCGACGATTGCCTGACCACTCACAAGTCAAGCAATCACCAAAAGAGACGAGAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGTGGCAATCTTACTAATATTTAATGACATTTGCCATTACAGTTTGTGTTCAGGGGTAGTAATGAGTGATGTGTGGTTGGACGTTGAACAGGTGACGGCGCCATCTGATCAACA‑TTCGTAACCgg                                                                                                                                                                                                 >  2:511779/10‑250 (MQ=255)
                                                                          ggggAGATTCAGGCGCGGCGATTGCCTGACCACTCACAATTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTTGTCGTTGAACAGGTGACGTCTCCATCTGTTCAACA‑TTCGTACCCTTGATTTCTCTATAACTATAGCTca‑ca                                                                                                                                                                       <  2:71031/251‑1 (MQ=255)
                                                                                    aGGCGCGGCGATTGCCTGACAACTCACAATTCCAGCAATCACCAAAAGAGACGATAATCCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGGTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGGGTTCCGGGATCGTATTGCAGGCTGTGTGGTTGGTCGTTGAACAGGTGAAGTCTCAATTTGATCAAAA‑TTCGTACCCGTGTTTTTTCTAAAACTCTAACT‑A‑CATGAACAACttt                                                                                                                                                            >  2:670481/1‑250 (MQ=255)
                                                                                     ggCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAg                                                                                                                                                           <  2:116430/251‑1 (MQ=255)
                                                                                            cgATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCgagg                                                                                                                                                      >  1:694917/1‑249 (MQ=255)
                                                                                             gATTGCCTTACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGTCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGGCGCCATCTGTGCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGc                                                                                                                                                   <  1:43664/251‑1 (MQ=255)
                                                                                                       ccACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTAc                                                                                                                                         >  1:481105/1‑251 (MQ=255)
                                                                                                        cACTCAGAAGTCCAGCAATCACCAGAAGAGACGATAATACCTGGAGCATAATATTAAGGTATTTTAAGTAAAGAGGCTATCTTACTAATATTTAATAAAATTTGCCATTACCGTTTGTGTTCAGGGGTAGTAATGAGGGCTGTGTGGTTGGTCATTGAACAGGTGAAGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATTACACA‑CA‑GAACAACTTAACGAGGAGCagcgcagaacg                                                                                                                                        >  2:449576/1‑242 (MQ=255)
                                                                                                           tCACAAGTCCAGCAATCACGGAACGATAAGTTAATACCTGTCGCATAATATTAAAGGATTTTATGTAAAGAGGATATCTTACGAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTAATGAGGACTGTGTTGTTGGTCGTTGAACAGGGGACGTCGCCATCTGTACAACA‑TTCGTACCCGTGATTTCTCTATACCTATAACTAA‑CA‑GAACATATTAGCGAGGAGAAGAGCGGTACGATc                                                                                                                                     >  2:539914/1‑251 (MQ=255)
                                                                                                             acaAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt                                                                                                                                   <  1:579931/251‑1 (MQ=255)
                                                                                                                             ccAAAAGCGACGAAAATCCCTGGCGCTTAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATAACATTTGCCATTAACGTTTGTGTTCAGGGGTCTTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCAATCTGTTCA‑CACTTCGTACCAGTGATTTCTCTATAACTATAAATCA‑CA‑GAACAACTTAGCGCGGCACAGAGCGGAACGATCGAGAGAAAcaaattaa                                                                                                                     >  2:289525/1‑242 (MQ=255)
                                                                                                                             ccAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGCCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGCGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGAACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaaa                                                                                                                    >  1:521829/1‑250 (MQ=255)
                                                                                                                               aaaaGAGACGATAATACCTGGCGCATAATATTAAAGGCTTTTATATAATAAGGCTATCTTACGAATATTTAATAACATTTGCCATTAAAGTTTGTGTTCAGGGGTCGTAATGAAGGCTGTGTGCTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCCACA‑TTCGACCCCGAGGTTTCTCTATAACACTAACTCAACA‑‑AACAACTTGCCGAGAAGCAGAGCGCTACGAtaatgatacaccacttaaagc                                                                                                                  >  2:514084/1‑230 (MQ=255)
                                                                                                                                                  ggCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCAtt                                                                                              <  1:221586/251‑1 (MQ=255)
                                                                                                                                                     gcATAATATCAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGATTGTGTTCCGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGACATCTGTTCAACA‑TTCGTACCCGTGTTTTCTCTATAACTATCACTCA‑CA‑GAACACCTTAGCGCGGAGCAGAGCGGTCCGATCGTGCGAAAACCATTAAAGCACGCTGCAGTGCGTGTCATAGtt                                                                                           >  2:697093/1‑251 (MQ=255)
                                                                                                                                                        taataTTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTa                                                                                                                                                            >  2:198199/1‑183 (MQ=255)
                                                                                                                                                        taataTTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTa                                                                                                                                                            <  1:198199/183‑1 (MQ=255)
                                                                                                                                                          atatTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCa                                                                                      <  1:446973/251‑1 (MQ=255)
                                                                                                                                                           tatTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGcc                                                                                       >  1:33372/1‑249 (MQ=255)
                                                                                                                                                                    aTTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGt                                                                            <  2:694917/251‑1 (MQ=255)
                                                                                                                                                                      tttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTg                                                                           <  1:280550/250‑1 (MQ=255)
                                                                                                                                                                       ttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGa                                                                          >  1:157627/1‑250 (MQ=255)
                                                                                                                                                                       ttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGAACGCTGCAGTGCGTGTCATTGTTAGCCAGAGGCGGCGTGa                                                                          >  1:385975/1‑250 (MQ=255)
                                                                                                                                                                               agagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTAt                                                                <  1:37818/251‑1 (MQ=255)
                                                                                                                                                                                gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATc                                                               >  1:166344/1‑251 (MQ=255)
                                                                                                                                                                                     tATCTTACTACTATTTAATGACCTTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTAGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCAAGCTGCAGTGCGGGTCATTGTAAGCAAGATGAGGCGTGACAC‑CTTTCTCCGGAc                                                          >  2:567460/1‑251 (MQ=255)
                                                                                                                                                                                       tCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCg                                                                                                                                                                                                                                            >  1:664698/1‑75 (MQ=255)
                                                                                                                                                                                       tCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCg                                                                                                                                                                                                                                            <  2:664698/75‑1 (MQ=255)
                                                                                                                                                                                         tttcTAATATTTAATGACAGTTGCCATTACCGGTTGTGTTCAGGGGTCGGAAGGAGGGCGGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTAGAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTTCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAAc                                                       <  2:277591/247‑1 (MQ=255)
                                                                                                                                                                                          tACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAAc                                                       <  1:13091/249‑1 (MQ=255)
                                                                                                                                                                                            cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGCTTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGCCAACG‑At                                                   >  2:262215/1‑250 (MQ=255)
                                                                                                                                                                                                   ttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGAcc                                            <  1:111323/250‑1 (MQ=255)
                                                                                                                                                                                                               gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAg                                                                                                         >  2:299655/1‑179 (MQ=255)
                                                                                                                                                                                                               gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAg                                                                                                         <  1:299655/179‑1 (MQ=255)
                                                                                                                                                                                                               gCCAGTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGGGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAAt                               <  2:506690/251‑1 (MQ=255)
                                                                                                                                                                                                                ccATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATATCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAAtt                              >  2:187726/1‑251 (MQ=255)
                                                                                                                                                                                                                  aTTACCGTTTGTGTTCCGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑TTCCTGCCCGATCGTCTGCAATTa                             >  2:92538/1‑250 (MQ=255)
                                                                                                                                                                                                                  aTTACCGTTTGTGTTCAGGGGACGTAATGAGGGCAGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGAGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑AGACTGACCGATCGTCTGCAATTa                             >  2:462885/1‑250 (MQ=255)
                                                                                                                                                                                                                    tACCGTTTGTGTTCAGGGGTCGTAATGTGGGCTGCGTGGTTGGTCGTTGCGCAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACACCTTAGCGCGGAGCAGAGCGGTACGATCGTGAGCAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAACCAGATGCGGCGTCA‑CCGATTTTTCCGGACAACGCAT‑CTGACCTAACCTCAGCAATTAAAt                          >  2:184890/1‑251 (MQ=255)
                                                                                                                                                                                                                          ttGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGGTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCGGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTTA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAATTAAAtcatta                     <  1:670481/250‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGAAATTAAAtcattactc                  >  1:178989/1‑251 (MQ=255)
                                                                                                                                                                                                                             tgtTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGGGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATAATGACAGATCGTCCGCCATTAAAtcattactc                  >  2:488398/1‑250 (MQ=255)
                                                                                                                                                                                                                              gtTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAATTAAAtcattactcat                >  1:672422/1‑251 (MQ=255)
                                                                                                                                                                                                                              gtTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAATTAAAtcattactcat                >  1:183392/1‑251 (MQ=255)
                                                                                                                                                                                                                                tCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGACCAGGTGCCGTCGACATCTGTTCAAAA‑TTCGTACCCGAGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGGGTCATCGTTAGCCAGATGCGGCGGGA‑ACGCTTTATCCGGGCAACG‑ATACTGACCTGTCGTCTGCAATTAAAtcattacccatt               >  2:361826/1‑250 (MQ=255)
                                                                                                                                                                                                                                          tAATGAGGGCTGTGTGGTTGGTCGTTGCACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCAATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGAt    >  1:260406/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTc  >  1:459599/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
TTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACA‑TTCGTACCCGTGATTTCTCTATAACTATAACTCA‑CA‑GAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA‑ACGCTTTATCCGGACAACG‑ATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTC  >  NC_000913/1356958‑1357444

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