Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,173 T→C 5.2% intergenic (+92/+27) 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,1730TC 5.2% 281.9 / 5.0 96intergenic (+92/+27)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 T (37/54);  new base C (3/2);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

GTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATA  >  NC_000913/4353025‑4353319
                                                                                                                                                    |                                                                                                                                                  
gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCACAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCt                                                                                                                                                    >  2:711269/1‑149 (MQ=255)
gTCGGAATAACCTGGCCAATGAGTCTGCCTGGGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCt                                                                                                                                                    <  2:81653/149‑1 (MQ=255)
   ggAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTTGCCTGAc                                                                                                                                                 <  1:249489/149‑1 (MQ=255)
    gAATACCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTGTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACt                                                                                                                                                >  2:260701/1‑149 (MQ=255)
    gAATAACCTGGCCAATGAGTCTGCCTGTGGGATTACTTTTTTCTCTCTTTTAGCAGCTGTAGTTGAGTGGACCCGCTGGTTATCGAGAAACGAACCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGACGCAGTCgcctgcct                                                                                                                                                >  2:19046/1‑149 (MQ=255)
    gAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACt                                                                                                                                                >  2:464441/1‑149 (MQ=255)
    gAATAACCTGCCAAATGAGGCTGCCTGTGGCATTACTTTTTTCTCTGTTTTATGCGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGAGGAACGATGAACGAAGCAGTCGCCTGACt                                                                                                                                                <  2:230428/149‑1 (MQ=255)
       tAACCTGGCCACTGACTCTGCCTGTTGCATTACTTTTTTCTATCTTTTAGGAGCTCTAGTTGAGTGGACCCGCTGTTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTAGCGGAACGATGAACGAAGCAGTCGCCTGACTGCt                                                                                                                                             <  1:1044777/149‑1 (MQ=255)
       tAACCTGGCCAAGGAGGCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCt                                                                                                                                             <  1:843588/149‑1 (MQ=255)
        aaCCTGGCCTATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGATGATGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGCAGCAGTCGCCTGACTGCtt                                                                                                                                            <  1:984161/149‑1 (MQ=255)
        aaCCTGGCCAATGAGTCTGCCTGTGTCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACACGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCtt                                                                                                                                            <  1:927835/149‑1 (MQ=255)
        aaCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCtt                                                                                                                                            <  2:1130579/149‑1 (MQ=255)
          ccTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCg                                                                                                                                          <  1:141598/149‑1 (MQ=255)
          ccTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCg                                                                                                                                          <  1:486123/149‑1 (MQ=255)
          ccTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGCTATCGTGAAACGAAGCGTCTCTCCTCAAAACAGGATTTGCGGAACGATCAACGAAGCAGTCGCCTGACTGCTTCg                                                                                                                                          <  1:309336/149‑1 (MQ=255)
          ccTGGCCAATGAGTCTGCCTGTGGAATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCg                                                                                                                                          <  2:925165/149‑1 (MQ=255)
           cTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCGt                                                                                                                                         <  2:473154/149‑1 (MQ=255)
            tGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGtt                                                                                                                                        >  1:1106038/1‑149 (MQ=255)
              gCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCa                                                                                                                                      <  1:456080/149‑1 (MQ=255)
               ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCAt                                                                                                                                     <  1:843180/149‑1 (MQ=255)
                cAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAATCAGTCGCCTGACTGCTTCGTGCaat                                                                                                                                    >  1:97345/1‑147 (MQ=255)
                cAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCAtt                                                                                                                                    >  1:713387/1‑149 (MQ=255)
                cAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCAtt                                                                                                                                    >  1:465952/1‑149 (MQ=255)
                  aTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGttaattaa                                                                                                                                  >  1:677758/1‑149 (MQ=255)
                    gAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAg                                                                                                                                <  1:136524/149‑1 (MQ=255)
                       tCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGa                                                                                                                             >  1:1076831/1‑149 (MQ=255)
                             tgtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCATCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAattta                                                                                                                       >  1:28296/1‑149 (MQ=255)
                                ggCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatttattt                                                                                                                    <  1:36776/149‑1 (MQ=255)
                                 gCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGTTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatttattta                                                                                                                   >  1:1075926/1‑148 (MQ=255)
                                 gCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTc                                                                                                                   >  1:383329/1‑149 (MQ=255)
                                  cATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCt                                                                                                                  <  2:1014384/149‑1 (MQ=255)
                                    ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                                <  1:416673/149‑1 (MQ=255)
                                    ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                                <  1:454351/149‑1 (MQ=255)
                                      aCTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTgg                                                                                                              >  2:256739/1‑149 (MQ=255)
                                        ttttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                            >  2:802184/1‑149 (MQ=255)
                                          ttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGACTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgcgg                                                                                                          >  1:664296/1‑148 (MQ=255)
                                               tctctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGc                                                                                                     <  2:901005/149‑1 (MQ=255)
                                                 tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACATGAATTGAGGAACTATGAACGAATCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGAGCATCGCCg                                                                                                   >  1:230261/1‑149 (MQ=255)
                                                  ctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCgg                                                                                                  <  2:961007/149‑1 (MQ=255)
                                                  ctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCgg                                                                                                  <  1:360688/149‑1 (MQ=255)
                                                   ttttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCggg                                                                                                 <  1:690013/149‑1 (MQ=255)
                                                     ttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTGCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGaa                                                                                               <  2:1106038/149‑1 (MQ=255)
                                                       aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                             >  2:464149/1‑149 (MQ=255)
                                                       aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTGCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                             >  1:679923/1‑149 (MQ=255)
                                                        ggAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAagag                                                                                            >  1:677598/1‑149 (MQ=255)
                                                        ggAGCTGTAGTTGAGTGGACCCGCGGGTTATCGTGAAACGAAGCGGCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAagag                                                                                            <  2:92067/149‑1 (MQ=255)
                                                            cTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAAAGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATa                                                                                        <  2:668795/149‑1 (MQ=255)
                                                              gTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAAc                                                                                      <  1:608596/149‑1 (MQ=255)
                                                                 gTTGAGTGGCCCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTc                                                                                   <  1:739843/149‑1 (MQ=255)
                                                                 gTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGTATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTc                                                                                   <  2:728253/149‑1 (MQ=255)
                                                                 gTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTc                                                                                   <  1:902698/149‑1 (MQ=255)
                                                                   tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                                 <  1:900503/149‑1 (MQ=255)
                                                                   tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                                 <  2:180761/149‑1 (MQ=255)
                                                                   tGAGTGGACCCGCTGGTGATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACAGCTTCGTTCATTAAAGTGAAATATATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                                 <  2:642953/149‑1 (MQ=255)
                                                                     aGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGa                                                                               >  1:838385/1‑149 (MQ=255)
                                                                        ggACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATa                                                                            >  1:566000/1‑149 (MQ=255)
                                                                        ggACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATa                                                                            <  1:154461/149‑1 (MQ=255)
                                                                             cGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATg                                                                       >  2:1176015/1‑149 (MQ=255)
                                                                              gCTGGTTCTCGTGACACGACGCGTCTCTCCTCAAACCTGGATTTGCGGACCGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATATCTGTGGGCGCATCGCCGGGAAGAGAATACCGTCGCGAATAGTATgg                                                                      <  2:713387/149‑1 (MQ=255)
                                                                              gCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATgg                                                                      >  2:753603/1‑149 (MQ=255)
                                                                               cTGGGTATCGTGAAACGAAGCGTCTCTACTCAAAACTGGATTTGCGGAACGATGAACGACGCAGTCGCCTGACTTCTTCGTTCAATAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTAGCGATGAGTATGGc                                                                     >  2:770613/1‑149 (MQ=255)
                                                                                gggatctCGTGAAACGAAGCGTCTCTCCTCAACACTGGCTTTGCCGACCGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCt                                                                    <  2:736101/143‑1 (MQ=255)
                                                                                  gTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGt                                                                  <  1:578049/149‑1 (MQ=255)
                                                                                    tATCGTGCAACGAATCTTCTCACGTCTACTCTGGATTTGCTTAAGGATTAACGAAGGATTCGCATGACTGCTGCGTTCATTAAAGTGAAATTTATTTCTGTGGGCTCATCGCCGGGAAGAGAATAACGTCGCGAATCGTATGGCTGTTa                                                                <  2:1124482/149‑1 (MQ=255)
                                                                                       cGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTa                                                             <  1:464441/149‑1 (MQ=255)
                                                                                         tGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGCCTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                           >  2:676083/1‑149 (MQ=255)
                                                                                         tGAAACGAAGCGTCACTCCTCAAAACTGTATTTGCGGAACGATGACCGAAGCATTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGATAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                           <  2:496500/149‑1 (MQ=255)
                                                                                          gAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAAc                                                          <  1:708838/149‑1 (MQ=255)
                                                                                            aaCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAg                                                        >  1:300739/1‑149 (MQ=255)
                                                                                             aCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGc                                                       <  2:838385/149‑1 (MQ=255)
                                                                                             aCGAAGCGTCTCTCCTCAAAACTGGAATTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGc                                                       <  2:569331/149‑1 (MQ=255)
                                                                                               gAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCAt                                                     >  1:73432/1‑149 (MQ=255)
                                                                                                 aGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGATTAGTATGGCTGTTAGTAAACAGCATaa                                                   >  2:518709/1‑149 (MQ=255)
                                                                                                 aGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATaa                                                   >  2:781525/1‑149 (MQ=255)
                                                                                                   cggctctcCTCAAAACTGGATTTGCTGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCGTAATc                                                 <  2:560923/146‑1 (MQ=255)
                                                                                                   cGTCTCTCCTCAAAACTGGATTTCCGTAACGATGAACGAAGCAGTCGACTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCAGTGTGCCCATCGCCGTGAAGAGAATAACGTCGCGAATACAATGGCTGTTAGTAAACAGCATAATc                                                 <  1:177501/149‑1 (MQ=255)
                                                                                                       tctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTc                                             >  1:251609/1‑149 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCgg                                           >  2:21887/1‑149 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAACAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCAACGCAGGGAAGAGAATAACGGCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAGTCgg                                           >  2:60204/1‑149 (MQ=255)
                                                                                                          cctcAAAACTTTATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAATAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGt                                          <  2:344947/149‑1 (MQ=255)
                                                                                                          cctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGt                                          <  2:383329/149‑1 (MQ=255)
                                                                                                          cctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGt                                          <  2:1075926/149‑1 (MQ=255)
                                                                                                          cctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGt                                          <  2:612381/149‑1 (MQ=255)
                                                                                                                         tGCGGAACGATGAACGAAGCAGGCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                           <  2:424778/149‑1 (MQ=255)
                                                                                                                          gCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAg                          >  1:393211/1‑149 (MQ=255)
                                                                                                                              aaCGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCa                      >  1:971819/1‑149 (MQ=255)
                                                                                                                              aaCGAGTAACGAAGCAGGCGCCTGACTGCGGCGTTCATTAAAGTGAAAATTATTGCGGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCa                      <  2:1076831/149‑1 (MQ=255)
                                                                                                                               aCGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTTTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAg                     <  1:1016461/149‑1 (MQ=255)
                                                                                                                                 gATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGTAATTTATTTCTGTGGGCGCATAGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCg                   >  1:654919/1‑149 (MQ=255)
                                                                                                                                      aCGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTgg              >  1:1221/1‑149 (MQ=255)
                                                                                                                                      aCGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAACCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTgg              >  2:235108/1‑149 (MQ=255)
                                                                                                                                       cGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAAAAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGCTTCCAATACCCAGACCAGCGGTTggc             >  1:276999/1‑149 (MQ=255)
                                                                                                                                        gAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg            <  2:3181/149‑1 (MQ=255)
                                                                                                                                        gAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGAGACCAATACCCAGACCAGCGGTTggcg            >  1:757928/1‑149 (MQ=255)
                                                                                                                                         aaGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAATCGGTCGATACCAATACCCAGACCAGCGGTTggcgg           <  2:677598/149‑1 (MQ=255)
                                                                                                                                             aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAAATTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGa       <  2:406862/149‑1 (MQ=255)
                                                                                                                                               tCGCCTGACTGCTTCGTTCATTAAAGTGATATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATACTCATTCGGTCGAAACCAATACCCAGACCAGCGGTTGGCGGCAGAcc     >  2:326551/1‑149 (MQ=255)
                                                                                                                                               tCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAATCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGAcc     <  2:918709/149‑1 (MQ=255)
                                                                                                                                                  ccTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGCCCata  >  1:1075623/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATA  >  NC_000913/4353025‑4353319

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