Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,171 C→G 5.4% intergenic (+90/+29) ghoT → / ← lysU toxin of GhoTS toxin‑antitoxin pair; membrane‑lytic protein; stimulator of persister cell formation/lysine tRNA synthetase, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,353,1710CG 5.4% 254.7 / 5.8 74intergenic (+90/+29)ghoT/lysUtoxin of GhoTS toxin‑antitoxin pair; membrane‑lytic protein; stimulator of persister cell formation/lysine tRNA synthetase, inducible
Reads supporting (aligned to +/- strand):  ref base C (39/31);  new base G (2/2);  total (41/33)
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

CGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGC  >  NC_000913/4353027‑4353311
                                                                                                                                                 |                                                                                                                                            
cGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGa                                                                                                                                          <  1:279839/149‑1 (MQ=255)
 ggAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGAc                                                                                                                                         <  4:514769/149‑1 (MQ=255)
  gAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACt                                                                                                                                        >  2:409132/1‑149 (MQ=255)
   aaTAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTg                                                                                                                                       >  3:576195/1‑149 (MQ=255)
   aaTAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTg                                                                                                                                       >  1:284282/1‑149 (MQ=255)
      aaCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑AGGCGACTGCtt                                                                                                                                    <  2:82368/149‑1 (MQ=255)
       aCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTc                                                                                                                                   <  3:212164/149‑1 (MQ=255)
       aCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTc                                                                                                                                   >  4:197705/1‑149 (MQ=255)
       aCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTc                                                                                                                                   >  4:454571/1‑149 (MQ=255)
                 tGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTaaa                                                                                                                         <  4:576195/149‑1 (MQ=255)
                     tCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAACACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGa                                                                                                                     >  4:631603/1‑149 (MQ=255)
                                  ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCTTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTTCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                        >  4:508737/1‑149 (MQ=255)
                                  ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                        >  4:665656/1‑149 (MQ=255)
                                  ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                        >  3:725969/1‑149 (MQ=255)
                                    aCTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTgg                                                                                                      >  4:221543/1‑149 (MQ=255)
                                      ttttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                    >  3:1732/1‑149 (MQ=255)
                                             tctctTTTAGGAGCTGTAGTTGTGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAat                                                                                                                  <  1:559449/128‑1 (MQ=255)
                                             tctctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAat                                                                                                                  >  2:559449/1‑128 (MQ=255)
                                               tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCt                                                                                                                                     <  3:250035/107‑1 (MQ=255)
                                               tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCt                                                                                                                                     >  4:250035/1‑107 (MQ=255)
                                               tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTaaa                                                                                                                         <  1:30812/119‑1 (MQ=255)
                                               tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTaaa                                                                                                                         >  2:30812/1‑119 (MQ=255)
                                               tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCg                                                                                           >  3:735268/1‑149 (MQ=255)
                                                      ggAGCTGTAGTTGAGTGGACCCGCAGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTCATTTCTGTGGGCGCATCGCCGGGAagag                                                                                    >  3:640478/1‑149 (MQ=255)
                                                            gTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                     <  1:492373/142‑1 (MQ=255)
                                                            gTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                     >  2:492373/1‑142 (MQ=255)
                                                                ttGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcg                                                                          >  4:647299/1‑149 (MQ=255)
                                                                 tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                         <  3:665656/149‑1 (MQ=255)
                                                                    gTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGaa                                                                      <  1:118770/149‑1 (MQ=255)
                                                                         cccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTa                                                                 <  4:375636/149‑1 (MQ=255)
                                                                          ccGCTGGTTATCGTGAACCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAGTGATGAACGAAGCAGTC‑GCATGACTGCTTCGTTCATTAAAGTAAAATTTATTTCTGTGGCCGCATCGCCGGGAAGAGAATAACGTCATGAATACTAt                                                                >  2:644319/1‑149 (MQ=255)
                                                                          ccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTAt                                                                <  3:370418/149‑1 (MQ=255)
                                                                            gCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATgg                                                              <  3:453255/149‑1 (MQ=255)
                                                                               ggTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                        <  1:348443/104‑1 (MQ=255)
                                                                               ggTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                        >  2:348443/1‑104 (MQ=255)
                                                                                  tATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTa                                                        >  2:749035/1‑149 (MQ=255)
                                                                                  tATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTa                                                        <  2:221798/149‑1 (MQ=255)
                                                                                    tCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGt                                                      <  2:284282/149‑1 (MQ=255)
                                                                                       tGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTa                                                        <  4:358131/144‑1 (MQ=255)
                                                                                       tGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTa                                                        >  3:358131/1‑144 (MQ=255)
                                                                                       tGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                   >  2:454652/1‑149 (MQ=255)
                                                                                            cGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCa                                              <  3:673293/149‑1 (MQ=255)
                                                                                            cGAAGCGTCTCTCCACAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCa                                              <  2:465298/149‑1 (MQ=255)
                                                                                                     tctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTc                                     >  4:446349/1‑149 (MQ=255)
                                                                                                     tctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTc                                     <  2:564497/149‑1 (MQ=255)
                                                                                                        cctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGt                                  <  2:385271/149‑1 (MQ=255)
                                                                                                         ctcAAAACTGGATTTGCTGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAACAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                 <  4:640478/149‑1 (MQ=255)
                                                                                                         ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                 >  3:227855/1‑149 (MQ=255)
                                                                                                         ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                 >  3:574630/1‑149 (MQ=255)
                                                                                                         ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                 <  3:221543/149‑1 (MQ=255)
                                                                                                         ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                 >  2:494715/1‑149 (MQ=255)
                                                                                                          tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAt                                       <  3:580864/142‑1 (MQ=255)
                                                                                                          tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAt                                       >  4:580864/1‑142 (MQ=255)
                                                                                                          tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                                >  3:380565/1‑149 (MQ=255)
                                                                                                          tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                                >  2:536119/1‑149 (MQ=255)
                                                                                                          tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                                >  2:68323/1‑149 (MQ=255)
                                                                                                             aaaCTAGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGATCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTGAACAGCATAATCATTCGGTCGATa                             >  2:775563/1‑149 (MQ=255)
                                                                                                               aCTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAc                            <  1:655263/148‑1 (MQ=255)
                                                                                                               aCTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCAGTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAc                            >  2:655263/1‑148 (MQ=255)
                                                                                                                   gATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATa                       <  1:749035/149‑1 (MQ=255)
                                                                                                                     tttGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                     <  4:735268/149‑1 (MQ=255)
                                                                                                                                  gAACGAAGCAGTCAGGC‑GACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAtt                                      >  3:10169/1‑119 (MQ=38)
                                                                                                                                  gAACGAAGCAGTCAGGC‑GACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAtt                                      <  4:10169/119‑1 (MQ=38)
                                                                                                                                  gAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAtt                                      >  4:387089/1‑119 (MQ=255)
                                                                                                                                  gAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAtt                                      <  3:387089/119‑1 (MQ=255)
                                                                                                                                   aaCGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTg       >  3:489240/1‑149 (MQ=255)
                                                                                                                                    aCGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTgg      <  1:68323/149‑1 (MQ=255)
                                                                                                                                      gAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCTCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                         >  4:495650/1‑128 (MQ=255)
                                                                                                                                      gAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCTCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                         <  3:495650/128‑1 (MQ=255)
                                                                                                                                      gAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGtt                                                         >  2:761577/1‑96 (MQ=255)
                                                                                                                                      gAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGtt                                                         <  1:761577/96‑1 (MQ=255)
                                                                                                                                      gAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg    >  2:671764/1‑149 (MQ=255)
                                                                                                                                      gAAGCAGTC‑AGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGATAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg    >  1:344002/1‑149 (MQ=255)
                                                                                                                                        aGCAGTC‑ACGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCTAGACTAGCGGTTggcggc  >  4:357027/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTC‑GCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGC  >  NC_000913/4353027‑4353311

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