Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,172 C→G 5.2% intergenic (+91/+28) 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,1720CG 5.2% 272.1 / 7.0 97intergenic (+91/+28)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 (37/54);  new base G (3/2);  total (40/57)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

GGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATA  >  NC_000913/4353024‑4353319
                                                                                                                                                    |                                                                                                                                                   
ggTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGcc                                                                                                                                                     <  1:330581/149‑1 (MQ=255)
 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)
                                                                                                                                                    |                                                                                                                                                   
GGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATA  >  NC_000913/4353024‑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: ‑