Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,592,188 C→T H178Y (CAT→TAT)  tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,1880CT98.3% 184.1 / ‑4.4 60H178Y (CAT→TAT) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/30);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

ATTTCGCTGAAACAGGCGGAAAAAAA‑CTGGGCGGATTACGAAGCGTT‑GCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGC  >  NC_000913/4591942‑4592434
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
aTTTCGCTGAAACAGGCGGAAAAAAAACTGGGCAGATTACGAAGCGTT‑GCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGGGCAGGCAAAATCAACGAGTTCTTTGATCAGCAGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGGGCAAAACGATCGGCTCTATg                                                                                                                                                                                                                                                        >  2:186427/1‑251 (MQ=255)
  ttcttataaaCAGGCGGAAAAAAT‑CTGGGCAGATTACGAAGCGTT‑GCCGCGTGACCCGCGTCAGAGCGCCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCTGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGGGGCTTACATGGAGCAAAACGATCGGCTCTATGat                                                                                                                                                                                                                                                      <  2:617465/243‑1 (MQ=255)
             aGGCGGAAAAAAA‑CTGGGCGGATTACGAATCGTT‑GCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGc                                                                                                                                                                                                                                          >  1:346894/1‑251 (MQ=255)
              ggCGGAAAAAAC‑CTGGGCAGATTACGAAGCGTT‑GCCGCGTAACCCGCGTCAGAGCACCGCCGCAGCGGCAGATATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCg                                                                                                                                                                                                                                         <  2:175786/251‑1 (MQ=255)
              ggCGGAAAAAAA‑CGGGGCAGATTAAGAA‑CGTTTGCCGTGCGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTAGCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCg                                                                                                                                                                                                                                         <  1:21888/251‑1 (MQ=255)
                       aaa‑CTGGGCAGATTACGAAGCGTT‑GCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATg                                                                                                                                                                                                                                >  2:597824/1‑251 (MQ=255)
                                        gAAGCGTG‑GCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGc                                                                                                                                                                                                                <  1:604977/251‑1 (MQ=255)
                                               t‑GCCGCGGGCCCCGCGTCAGGGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtgtg                                                                                                                                                                                                          <  2:565967/250‑1 (MQ=255)
                                                      gTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGGGGATTCtg                                                                                                                                                                                                   >  2:200548/1‑251 (MQ=255)
                                                       tGACCCGAGTCAGAGCACCGCCGCAGCGGCAGAGATCAACCGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACACTGCCTCCTACAGCCAGGCGCGGTGGATTCtgg                                                                                                                                                                                                  >  2:627561/1‑251 (MQ=255)
                                                         ccccgtgtcagggacccgccgCAGCGGCAGAGATCAAACGTAATGACGATATTTATCACAATGCGCTGGCGGAGATGATCCAACTATTAGGTGCAGGCAAAAACAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCGATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                                                                                                                                <  2:623581/236‑1 (MQ=255)
                                                         ccccGCGTCAGAGTACCGACGCAGCGGCAGAGATCAAACGTAATTACGGTATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTGTGATCAGCCGACTCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATACCAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                                                                                                                                <  2:474489/250‑1 (MQ=255)
                                                                tCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCCCAATGCGCTGGCGGAGCTGATCCAACGATTAGGGGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCAAGGGATAGCAGGCCGGTTTCGAGAAGCAGTATGTGGCGTACAGGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCAATAACATTGCCACCTACAGCCAGGCGGTGGGGATTCTGGGGGGCggg                                                                                                                                                                                          >  2:214044/1‑248 (MQ=255)
                                                                 cAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGCCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtgat                                                                                                                                                                                        >  1:623099/1‑251 (MQ=255)
                                                                    agCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTgatgat                                                                                                                                                                                     <  1:158947/251‑1 (MQ=255)
                                                                      ccccgccgCAGCGGCAGAGATCAAACGTAATTACGAGATTTATCAAAATGCGCGGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAGTCTGGTGGGCGTGATGAtcg                                                                                                                                                                                   <  2:113084/249‑1 (MQ=255)
                                                                           ccgcAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTa                                                                                                                                                                              <  1:597824/251‑1 (MQ=255)
                                                                            cgcAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCCACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGCCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCCAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTCCAGCCAGGc                                                                                                                                                                                                                >  1:121033/1‑216 (MQ=255)
                                                                            cgcAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCTATGTGGATTCTGGTGGGCGTGATGATCGTCGTAc                                                                                                                                                                             <  1:214044/251‑1 (MQ=255)
                                                                                    cAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTc                                                                                                                                                                     >  1:20592/1‑251 (MQ=255)
                                                                                           cAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTc                                                                                                                                                               <  1:375028/250‑1 (MQ=255)
                                                                                           cAAAAGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATAAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTCCTGGCGGTCCTCTTcc                                                                                                                                                              >  1:470428/1‑250 (MQ=255)
                                                                                                  aaTTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCttttgccgctg                                                                                                                                                       >  1:670529/1‑242 (MQ=255)
                                                                                                       cGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAACTTCGCCGTCT‑GGTTCg                                                                                                                                                 >  1:164928/1‑251 (MQ=255)
                                                                                                                 cacaATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑GATTTAAGc                                                                                                                                      >  1:183639/1‑251 (MQ=255)
                                                                                                                   caATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATTATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCt                                                                                                                                    <  1:625662/251‑1 (MQ=255)
                                                                                                                          cTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTTGG‑GATTAAAGCCTccactgg                                                                                                                             >  1:637739/1‑245 (MQ=255)
                                                                                                                              cGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAgcg                                                                                                                         <  1:563079/251‑1 (MQ=255)
                                                                                                                                      aTCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAACACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCCGGCGATGAGGATTCTGGTGGGCGTGATGCTCGTCGTACTGGCGGTGATCTTATCCGTCG‑GGTTCTG‑TTTTAAAGCCTCCCGGGCAGCGcccgtgc                                                                                                                  >  1:608472/1‑245 (MQ=255)
                                                                                                                                        ccAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATc                                                                                                               <  1:200548/251‑1 (MQ=255)
                                                                                                                                        ccAACTATTAGGGGCAGGTCAAATCAACGAGGTGTTGGAGCAGCCGACCCAGGGATATAAGGACGGTTTCGAGAAGCAGGATGTGGCTTACATGGAGCAAAACGATCGGCTCTAGGATATCGCCGTCAGCGCTAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGGGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAAt                                                                                                                <  2:179052/250‑1 (MQ=255)
                                                                                                                                         cAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGCCT‑GGTTCGG‑TTTTAAAGCCTCGCTGGGAGCGCCAATGaaacc                                                                                                              >  1:241371/1‑248 (MQ=255)
                                                                                                                                            ctattAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGGTCGG‑TATTAAAGACTCGCTGGTAGAGCCAATGAATCTCCt                                                                                                           >  1:325098/4‑251 (MQ=255)
                                                                                                                                                tAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGACATCTTTGCCGGCT‑GGTTCGGGTTTTAAA‑CCTCGCGGGTAGGGCCACTGAATATCCTGAtt                                                                                                       >  1:492113/1‑251 (MQ=255)
                                                                                                                                                  ggTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGCAGCGCCAATGAATCGCCTGATTGa                                                                                                     <  1:116372/251‑1 (MQ=255)
                                                                                                                                                  ggTGAAGGCAAAATCAACGAGTTCTGTGATCAGACGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTTGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGa                                                                                                     <  2:10806/251‑1 (MQ=255)
                                                                                                                                                       aGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTTGCCCGCT‑GGGTCGG‑GGTTAAAGCCTCGCTGGGTGGGCCCATGTATCGCCTGGTTggttgg                                                                                                 >  2:250081/1‑245 (MQ=255)
                                                                                                                                                           aaaaTCAACGAGTTCTTTGAGAAGCCGTCCCAGGGATATCTGGTCGGGTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACTAACGGCTCTATGATATCGGCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTTGCGCCAATGAATCGCCTGATTGACAGCATTCg                                                                                            <  2:64087/251‑1 (MQ=255)
                                                                                                                                                                cAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCCGCT‑GGTTCGG‑GATTAAAAACCCCCCGGTGGAGCCCATGAAACGCCTGATTGACAGCCATCGTCat                                                                                        >  1:238150/1‑250 (MQ=255)
                                                                                                                                                                  aCGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGTGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTAGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCata                                                                                       <  2:240741/249‑1 (MQ=255)
                                                                                                                                                                     aGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTAGGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑a                                                                                 <  2:536509/251‑1 (MQ=255)
                                                                                                                                                                      gTTCTTTGATCAGCCTACCCAGGGATATCAGGACGGTGTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGAGCGGCTCTATGATATCGCAGTCAGCGATAACAATGCCGCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCT‑GGTTCGG‑TATTAAAGCATCGCTGGTAGCGCCAATGAATCGCCGGATTGACAGCATTCGTCAGTTTGC‑Ag                                                                                <  2:86895/251‑1 (MQ=255)
                                                                                                                                                                                 aGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGGGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTTGG‑TATTTAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGGTTTACCACCTTCGTCAAATTGC‑a                                                                                 >  1:438929/1‑239 (MQ=255)
                                                                                                                                                                                 aGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTGACATGGGGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑a                                                                                 <  2:438929/239‑1 (MQ=255)
                                                                                                                                                                                   ccGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGt                                                                   <  1:186427/251‑1 (MQ=255)
                                                                                                                                                                                    cGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCT‑GGTTCGG‑GATTAAAGCCCCGCTGGGTGCGCCAATGGATCGCCTGGTTGACAGCATTTGGCATATTGC‑AGGCGGCGATCTGGTg                                                                  >  1:134973/1‑251 (MQ=255)
                                                                                                                                                                                        ccAGGGATATCAGGTATGTTTCGGGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCGCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCGGATTTACAGCATTCGTCATATTGC‑AGGTGGCGATCTGGTGAAAc                                                              <  2:50273/251‑1 (MQ=255)
                                                                                                                                                                                          aGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCCGCT‑GGTTCGG‑TATTAAAGCCTCGCTGGGAGCGCCCAAGAATCGCCTGATTGTCAGCATTTGTCATATTTC‑AGGCGGCGGTCTgggggaaacc                                                            >  1:382573/1‑243 (MQ=255)
                                                                                                                                                                                            ggATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAACCGAt                                                          <  1:334470/251‑1 (MQ=255)
                                                                                                                                                                                                      aCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTGATTTTCGCGGGCTAGGTTCGG‑‑CTTAAAGCCTCACTGGTAGCGCCACAGAAACGCCTGGTTGAAAACACTCGTCAAATTGCCAGGCGG‑GGTCTGGTGGAACCCATTGAggggggt                                                >  1:119556/1‑246 (MQ=255)
                                                                                                                                                                                                                      gTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGGTCGG‑GATTAAAGCCCCGCCGGTAGCGGCAATGAAATGCCTGAGTGACAACATTCGTCCTATTGC‑AGGCGGCGATCTGGGGAAAACGGTTGGGGGGGGTGGCGCTAATgaggagg                                >  1:332489/1‑247 (MQ=255)
                                                                                                                                                                                                                                tACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTTGG‑GATTTAAGCCTCGCTGGTAGCGCCAATGGATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGGGAGCTTGTGAAAACGATTTAGGGGGGTGGCTCTAATGAGATgggggaaatggg                      >  1:557662/1‑243 (MQ=255)
                                                                                                                                                                                                                                  cATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCag                    <  1:290240/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aGCAAAACGATCGGCCCTAGGATATCGCCGTCAGCGATAACAATGCCTCCGACAGCCAGGCGATTTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAAGCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGt               <  2:557662/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaaCGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAAcc                                                             <  1:540568/202‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaaCGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTg            <  1:130292/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaCGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTa                                         <  1:92396/221‑1 (MQ=255)
                                                                                                                                                                                                                                                aTCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCT‑GGGTCGG‑TATTAAAGCCTCCCTGGTAGCGCCCATGAAACGCCCGCTTGACAGCATTCGTCATATTGC‑AGGCGGCGGTCTGGGGGAACCGATTGAGGGGGGTGGCTCTAATTAGATGGGGGAACTGGCAGAGAGGTTTCGccct      >  1:240439/1‑249 (MQ=255)
                                                                                                                                                                                                                                                    gCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCTATTTGCCGCCT‑GGTTTGG‑GATTAAAGCCTCCCCGGTAGCGCCAACTAATCGCCCGGTTGACAGCAATCGTCATTTTGG‑AGGCGGAGGGCTGGTGAAAACGGTTGCGGTGGGTTATTCTAAAGAGAAGGGGCAACTGGCAGAGAGGTTGCGCCatattc  >  1:632376/1‑249 (MQ=255)
                                                                                                                                                                                                                                                    gCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTTGG‑GAATAAAACCTCGCTGGGAGGGCCAATGGATAGCCTGCTTGAAAGGATTCGTCATATTTC‑CGGCGGGGGTCTGGGGGAACCGGTTGTGGTGGGTGGGTCTAAAGAGAAGGGGGAACTGGCCGAGAGGTTggggcatctgc  >  1:275732/1‑241 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
ATTTCGCTGAAACAGGCGGAAAAAAA‑CTGGGCGGATTACGAAGCGTT‑GCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCT‑GGTTCGG‑TATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGC‑AGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGC  >  NC_000913/4591942‑4592434

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: ‑