Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,691 C→T 16.2% intergenic (+45/+20) yegJ → / ← yegK uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6910CT16.2% 100.8 / 3.0 59intergenic (+45/+20)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base C (31/17);  new base T (4/6);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAT  >  NC_000913/2151492‑2151924
                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
aatgtgggaaGGAGTGAAAAAAAAAGAAAGAATATTTTTCAGATTGGATGTAGCCGCGGGGGAATAAAATTTAAGGCGGTTATCCCATTTATCTTTTTTTGGTCACCTATCAAAAAGAAGAAGCCGCCGAGTTAAGAACAAATTTAGTCCCCTAACCAGCCTTTTTTTTTATTGTTTAAACGGGCGGTTGAATTTTTTTCTTCATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCa                                                                                                                                                                                                        <  2:1213417/230‑1 (MQ=255)
                  aaaaaaaaaaaaaaGTTATTCCAGAGTTGGTGTATCCGCGGGGGAAAAAAATTTAGGGCGGTTAAACTATTGATCTTTTATTGGCCCCTTATACAAAAAAAGAAGCCGACGAGAAAAGAGCAAAATTAGTTCGATAATCAGCCTTTATATTTTTTGTTTAAACTGGCGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                                                                                                                      <  2:1484154/231‑1 (MQ=255)
                   aaaaaagaaagaAGATTTTTCAGCTGGGATGTATACGCGCGGGAAAAAAATTTCCGGGGTTTATCCCTTTGATCTTTTATTGGTCCCCTATCAAAAAGAAGAGGCCAACGAGCTAAGAGAAAAATTAGTTCGCAAATCAGCCTTTATATTTATTGTTAAACCGGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatc                                                                                                                                                                                     <  1:1203877/235‑1 (MQ=255)
                                            tgtaTGTATACGGGGGGGAATAAAATTTAGGGCGTTTATCCCATTTATCCTTTATTGCTCCCCTAACAAAAAAAAAAAGCCGACGAGCTAAGAGAAAAATTGGTTCGCAAATAGCCCTTTATATTTATTGTTTAAACGGGCGGTTGAATTGTTTTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCggg                                                                                                                                                          <  1:874035/236‑1 (MQ=255)
                                                                                tataCCATTGACCTTTTATTGGTCACCTATCAAAAAAAAGAAGCCGACGAGCAAAGAGAAAAATTATTCCGCTAATCAGCCTTTATATTTATTGTTAAACCTGGGGGTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATaaa                                                                                                                      <  1:698095/239‑1 (MQ=255)
                                                                                        tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGTATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCCTACGCTCGTTTACCGCCGCGGTGGACAAAAAATGCTTTTACAATTCTGGGAATAAAACGAggtg                                                                                                            >  2:558064/1‑239 (MQ=255)
                                                                                        tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAAATGCTTTAACAATTTAGGTAATAAATCGAgttg                                                                                                            >  1:985986/1‑239 (MQ=255)
                                                                                         gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCTTGAGAAAACCACTCTTCCCCCACCGCCCGTTCCCACCGTTTATTATCATTCCTACGCCCGCTTTTCCCCGGGCGGGGCCAAAACTGCCTTAACAATTCAGCTAATAAATCGAgttgt                                                                                                           >  2:779101/1‑239 (MQ=255)
                                                                                               ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCCAACGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGGGCCAACGGTTTTTCATCCTACGTAACCCCGTTTACCCCCGGGCTGGACAAAAAATCCTTTAACAATTCCGGTAAACAATCGAGTTGTTCCCGc                                                                                                     >  2:978505/1‑239 (MQ=255)
                                                                                               ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAAACCCCCTGAGAAACTCACTCTTCCCACAGCCCCACTGCCAACGTTTTATCATCATCCCTCCGCCCGTTTTCCGCCGGGGTGGGCAAAAAAAGCTTTTACAATTTCGGTTATAAATCGCGTTGTTCCttc                                                                                                     >  1:1114988/1‑237 (MQ=255)
                                                                                                 aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCACGAGAAACTCCCTCTTCCCACAACGCCAGTGCCAACGCTTTATCATCATCCGTAAGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATTGAGTTGTTCCTTCGt                                                                                                   >  1:490302/1‑239 (MQ=255)
                                                                                                     gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTTATTTTTGTCTTCCACCGCCATTAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACCCTCGCTTACCGCCCGGCTGGACAAAAAATGCTTTTACAAATCAGGGAATAAATCGCGTTGTTCCTGCggt                                                                                                  >  2:280196/1‑234 (MQ=255)
                                                                                                     gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATTATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGGAATAAATCGAGGTGTTCCTGCGtt                                                                                                  >  1:753391/1‑236 (MQ=255)
                                                                                                       cACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCCACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCCACATCCGTACGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTTACAATTCCGGGAATAAAACGCGTTGTTCCTGCGTttctt                                                                                              >  2:1236316/1‑237 (MQ=255)
                                                                                                       cACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATGTTGTCTTCAACCGCCATGAAAAACTCACTCTGCCCACAACCCCAGTGCCAACGTTTTATCATCATCCCTACGCTCGGTTACCGCCCGGCGCGCCAAAAACCGCTTTACCAATTCAGGGAATAAATAAGGTTGTCCCCTCGTtgctg                                                                                              >  2:370670/1‑238 (MQ=255)
                                                                                                        aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTCTTTAACCGCCATGAGAAACTCACTCTGCCCACCGCGCCAGTGCCCACGTTTTATCCTCATCCGCACGCTCGTTTTCCGCCGGGCGTGCCAAAAACTTCTTTCACACTTCAGGGAATAAAATCAGGTGTTCCTCCGTttctgc                                                                                             >  2:1227508/1‑238 (MQ=255)
                                                                                                          cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCCAGTTGTTCCTGCGTTGCTgccg                                                                                           >  1:630231/1‑238 (MQ=255)
                                                                                                          cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAACGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTTGGCAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGGTTTTTCTGCGTTGCTgccg                                                                                           >  1:1627048/1‑238 (MQ=255)
                                                                                                          ataacaaaaaaaaaaaGCCAACGAGATAAGAAAAAATTTATTTCTATAACCAGCTTTTATATTTATTGTTAAACCTGGCGTTTGATTTTTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccgc                                                                                          <  1:521797/233‑1 (MQ=255)
                                                                                                            accaaaaaaaaaAGGCCGACGAGCAAAGATAAAAATGAGTCCTTAACTCAGCCTTTATTTTTATTTTTAAAACTGGGGTTTAATTTGTTGTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCa                                                                                        <  1:62107/235‑1 (MQ=255)
                                                                                                            aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGGGCCAACGGTTTATCATCATCCCGACGCTCCGTTACCGCCCGGCCGGGCAAAAAATGCTTTAACAAATCAGGGAATAAAAAGGGGTGTTCCTGCGTTGCTgccccc                                                                                         >  2:1505699/1‑238 (MQ=255)
                                                                                                             tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCCCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTTCCCCCGGGCTGGGCCAAAACTGCTTTAaaaaat                                                                                                                                 >  1:1256588/1‑197 (MQ=255)
                                                                                                             tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATTTTTTTTTCAACCCCCCTTAAAAACACACTCTGCCCACAACGACAATGCCAACGGTTTATCATCATCCGTAAGCTCGTTTACCGCCCGGACGGCCAAAAACTGCTTTAACAATATCGGAAATAAACCGAGTTGTTCCTTCGTTTCTGCCGCCat                                                                                       >  2:393885/1‑238 (MQ=255)
                                                                                                                       gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGATTTATTGTTTTCTTCAACAGCCATGAGAAACTCACTTTGCCCACAGCGCCCGAGCCAACGGTTTATTAACCTTCGGACGCTCCGTTCACCCCGGGCCGGGCCAAACCTGGTTTTAAAAATCATGGAATAAATACGGGTTTTCCTTGCTTGCAGCCGCCCTTCCAttttaa                                                                             >  2:240416/1‑235 (MQ=255)
                                                                                                                        aagCCGGCGAGCTAAGAAAAAATTTATTTCCATAACCAGCTTTTTTTTTTTTTTTTTAAACGGGGGGTTAATTTGTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaa                                                                            <  1:512202/239‑1 (MQ=255)
                                                                                                                         agCCGACGAGATAAGAGAAAAATTAGTTCGTTAATCAGCTTTTATATTTATTGTTTAACCGGGGGTTTAATTTTTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                                           <  1:280196/239‑1 (MQ=255)
                                                                                                                              acgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTTAATTGTTTTCTTTAACCGCCATGAGAAACTCAATCTGCCCACAGCGCCAGTGCCCCCGATTTATCATCATCCGTACGCTCGCTTACCGCCCGGGCGGGACAAAAATTCTTTAACAAATCCGGGAAAAAAACGAGTTGTTCCTCTGTTTGTGCCGCCAGTCCCTTGaaaaaaaggg                                                                      >  2:1406802/1‑234 (MQ=255)
                                                                                                                                          cAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTTAACCGCCATGAGAAAATCACTCTGCCCACAACGCCAGTGCCAAAGTTTTTTCATCATCCGTACGCCCGTTTAACGCCCGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTTCCGCCCGTTCATTGAAAAAAGGGGGAAAGAAAgggc                                                          >  2:120551/1‑237 (MQ=255)
                                                                                                                                            aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAAACGCCATGAGAAAATCACTCTGCCCACAACGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGTTTACCCCCGGGCTGGGCAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCCGGCCCTTGAAAAAAGGGGGAAAGAAAAGcaca                                                        >  1:955784/1‑239 (MQ=255)
                                                                                                                                             aaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTGTTGTCTTCCAACGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAAAGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAAATGCTTTTACAATTCTGGTAATAAATCGAGTTGTTCCTCCGTTGCTGCCGCCCGTCCATTTAAAAAAGGGGGAAAAAAAAGCACAt                                                       >  2:1586408/1‑239 (MQ=255)
                                                                                                                                                tAGTCCTCTAATCAGCCTTTATTTTTTTTTTATAACGGGGCGTTTGATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  1:120551/239‑1 (MQ=255)
                                                                                                                                                tAGTCCGCCAATCGCCCTTTATTTTTATTGTTAAACCTGGTGGTGAATTTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTCCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  1:370670/239‑1 (MQ=255)
                                                                                                                                                tAGTACGATAACCGGCCGTAATATTTGTTGTTAAAACTGGGGGTTAATTGTTTTCCTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  1:823773/239‑1 (MQ=255)
                                                                                                                                                  gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCCCCATGAGAAACTAAATCTGCCCCAAGCGCCCGTGCCAACGTTTTATAATCATCCGTTAGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAATTCAGGTAAAAAATCGAGTTGTTCCTGCGGTGCCGCCGCCCGTCCATTGAACAACCGGGTAAAGAAAAGCACATgagga                                                  >  1:640619/1‑239 (MQ=255)
                                                                                                                                                  gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCCCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAATTCCGGTAATAAATCGCGGTGTTCCTGCGTTGCTGCCGCCAGTCCCTTTAAAAAAGGGGTaaagaa                                                               >  1:923485/1‑226 (MQ=255)
                                                                                                                                                   ttCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCCACGGTTTTTCATCATCCGTACGCCCGTTTACCCCCGGGCCGGACCAAAACTGCTTTTAAAATTTAGGGAATAAATCGAGTTGTTCCTGCGTTGCTGCCCCCAGTCCATTTAAAAACGGGGTAAAGAAAGGCACATgaggag                                                 >  1:753207/1‑239 (MQ=255)
                                                                                                                                                     ctataacgaGCTGTTATTTTTATTGTTAAACCTGGGGGTTAAATTGTTGTCTACACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGaa                                               <  1:1550690/231‑1 (MQ=255)
                                                                                                                                                     cGCTAATCAGCCTTTATTTTTGTTTTTAAAACGGGCGGTTGATTTGTTTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  1:693092/234‑1 (MQ=255)
                                                                                                                                                             aGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCCGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCa                                        >  1:933228/1‑238 (MQ=255)
                                                                                                                                                             aGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACACCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTTACAATTCAGGGAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCCTTGAAAAAAGGGGGAAAGAAAGGCACATGAGGAGAATTATCCa                                        >  1:244877/1‑238 (MQ=255)
                                                                                                                                                              gcttttATTTTTTTTTTTTAAACTGGGGTTTAATTTTTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                      <  1:1505699/235‑1 (MQ=255)
                                                                                                                                                               tccttATTTTTTTTTTTTTAACGGGGCGTTTGATTGTTTGTTTCAACCCGCCATGAGAACCTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                      <  2:630231/236‑1 (MQ=255)
                                                                                                                                                               cgtttATTTTTTTTGTTTAACTGGCTGTTGTATTTTTTACTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                                     <  1:285245/236‑1 (MQ=255)
                                                                                                                                                               cgtttATATTTGTTTTTAAATCTGGTGTTTGAATTGTTGTTTCAACCCGCCATGAGAATCTCACTCTGCCCACAGCGCCAGTCCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGATCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                                     <  1:1630111/236‑1 (MQ=255)
                                                                                                                                                               ccGTTATATTTATTGTTTAAACTGGCGGTTGTATTTTTGTTTTCCACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATTATCCCTACCCCCGTTTACCGCCCGGCTGGGCAAAAACTGCTTTAACAATTCAGGGAATAAAACGAGTTGTTTCTCCGTTTCTGCCGCCCGTCCATTGAAAAACGGGGGAAAGAAAGCGACATGCGGGGAATTATc                                          >  2:693092/1‑234 (MQ=255)
                                                                                                                                                               agaaagatttttttttttAAAACGGGGGGTGAATTTGTTGCTTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                      <  1:1070370/231‑1 (MQ=255)
                                                                                                                                                                  ttATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCCACCGCCCTGAGAAACTCACTCTTCCCACAATCCCAGTGCCCACGGTTTTTTATAAACCGTAAGCCCGGTTACCGCCCGGGTGGGCAAAAACCGCTTTTAAAATTCAGGTAATAAATCGAGTTTGTTCCGTGTTGCTGCCCCCCGTCC‑‑TGAAAAACGGGGGAAAAAAAAGCCCCAGGGGAGAATTATCCAgactaatt                                >  2:967434/1‑232 (MQ=255)
                                                                                                                                                                       gtttttttttAAACCGGGGGGTTGATTTGTTGTCCAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTaa                              <  1:240416/232‑1 (MQ=255)
                                                                                                                                                                         tattGTTTAAACTGGCGGTTGAATTTTTGTCTTTCACCGCCATGAGAAACTCACTCTGCCCACCGCGCCAGTGCCCACGGTTTTTCATCATCCGTACGCTCGTTTACCCCCCGGGTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGGTGTTTCTGCGTTGCTGCCGCCCGTCCATTGTAAAACGGGGTAAAGAAAAGCACATGAGGAGAATTATTCAGCATAttttaag                            >  2:1120763/1‑234 (MQ=255)
                                                                                                                                                                            tGTTTAAACTGGCGGTTGAATTGTTGTTTTCAACCCCCCTGAGAAAATCACTCTTCCCACAGCGCCAGTGCCAAAGTTTTATCATCATACGTAAGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGGTGTTCCTGCGTTGTTGCCGCCAGTCCATTGAAAAAAGGGGTAAAGAAAGGCACATGAGGAGAAATATCCAGCATCTTTAAAGCCAg                        >  2:980512/1‑239 (MQ=255)
                                                                                                                                                                             tattAAACTGGGCGTTTAAGTGTTGTTCTCAACCCGCCATGAGAAACTCACTCTGCCCAAAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg                        <  1:859913/237‑1 (MQ=255)
                                                                                                                                                                              tttAAACGTGGGGGTGAATTTTTTTCTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGt                       <  1:885555/239‑1 (MQ=255)
                                                                                                                                                                               ttAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCCGTGGCAACGTTTTTTCATCATCCGTACCCCCGCTTTCCCCCGGGCGCGACAAAAACTGCTTTAACAATTCATGGAATAAATCGAGGTGTTCCTGCGCTTCTTCCGCCAGTCCATTTAAAAAAGGGGGAAAGAAAAGCACAAGAGGAGAATTATCCAACCTATTTAAAGCCAggc                      >  2:885555/1‑236 (MQ=255)
                                                                                                                                                                               ttAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCCAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGcaca                                                        >  1:581375/1‑204 (MQ=255)
                                                                                                                                                                                    gtGGCGGTAGATTTTTTGTCCTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACTTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGa                <  1:978505/238‑1 (MQ=255)
                                                                                                                                                                                        ggggTTAAATTTTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATa              <  1:1236316/236‑1 (MQ=255)
                                                                                                                                                                                                     ttgtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTaaag                                                                 >  1:207845/1‑175 (MQ=255)
                                                                                                                                                                                                       gtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAt  >  2:474019/1‑236 (MQ=255)
                                                                                                                                                                                                       gtCTTCAACCGCAATGAGAAACTCACTCTAACCACAGCGCCAGTGCCATCGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGCCACAAACTGCTTTAACATTTCCGGAAATAAATCGAATTGTTCCTGAGTTGAAGCCGCCAGTCCATTGAAAAACGGGGTAACGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCTCTGGATAGCATCCGTAAAt  >  2:1538222/1‑236 (MQ=255)
                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
AATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAT  >  NC_000913/2151492‑2151924

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