Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,173 T→C 5.3% 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.3% 418.1 / 6.6 127intergenic (+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 (60/59);  new base C (2/5);  total (63/64)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

TCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCG  >  NC_000913/4353026‑4353309
                                                                                                                                                    |                                                                                                                                        
tCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCt                                                                                                                                          >  4:102837/1‑148 (MQ=255)
 cGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTg                                                                                                                                         >  3:457658/1‑148 (MQ=255)
     aTAACCTGGCCGATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGTTGAACGAAGCAGTCGCCTGACTGc                                                                                                                                    >  1:797499/1‑149 (MQ=255)
       aaCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCtt                                                                                                                                  <  2:432900/149‑1 (MQ=255)
        aCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGTAACGATGAACGAAGCAGTCGCCTGACTGCTTc                                                                                                                                 <  2:23084/149‑1 (MQ=255)
          cTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGt                                                                                                                               >  2:1070520/1‑149 (MQ=255)
            ggCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTc                                                                                                                                 <  1:35686/145‑1 (MQ=255)
             gcgaaTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCa                                                                                                                            <  1:1038984/146‑1 (MQ=255)
              ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTATTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTGCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCAt                                                                                                                           <  1:453962/149‑1 (MQ=255)
               cAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCAtt                                                                                                                          <  1:447508/149‑1 (MQ=255)
                aaTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCAtt                                                                                                                          >  3:57602/1‑148 (MQ=255)
                aaTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTa                                                                                                                         >  2:30135/1‑149 (MQ=255)
                aaTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTa                                                                                                                         >  1:863234/1‑149 (MQ=255)
                aaTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTa                                                                                                                         >  1:254193/1‑149 (MQ=255)
                 aTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTa                                                                                                                         >  3:225923/1‑148 (MQ=255)
                 aTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTa                                                                                                                         >  3:400646/1‑148 (MQ=255)
                  tGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTaa                                                                                                                        <  3:723006/148‑1 (MQ=255)
                      tCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTg                                                                                                                    >  4:215240/1‑148 (MQ=255)
                       cTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGCAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGTCTGACTGCTTCGTTCATTAAAGTGa                                                                                                                   >  3:640937/1‑148 (MQ=255)
                        tGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGaaa                                                                                                                 <  1:300337/149‑1 (MQ=255)
                            tgtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAattt                                                                                                              <  4:52211/148‑1 (MQ=255)
                             gtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCCGACTGCTTCGTTCATTAAAGTGAAatttat                                                                                                            <  1:22898/149‑1 (MQ=255)
                                 cATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTc                                                                                                         >  3:917063/1‑148 (MQ=255)
                                  aTTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtg                                                                                                       >  2:829388/1‑149 (MQ=255)
                                   ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGtt                                                                                                                              >  2:66623/1‑125 (MQ=255)
                                   ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTaaa                                                                                                                       <  2:404851/132‑1 (MQ=255)
                                      cTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTggg                                                                                                   <  2:951436/149‑1 (MQ=255)
                                       ttttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                  >  2:653331/1‑149 (MQ=255)
                                        tttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                  >  3:664680/1‑148 (MQ=255)
                                        tttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                  >  4:979476/1‑148 (MQ=255)
                                              tctctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCg                                                                                            <  4:811005/148‑1 (MQ=255)
                                                tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGcc                                                                                          >  3:221697/1‑148 (MQ=255)
                                                tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCg                                                                                         <  2:797499/149‑1 (MQ=255)
                                                tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCACAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGcc                                                                                          >  4:935276/1‑148 (MQ=255)
                                                 ctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCg                                                                                         <  4:457658/148‑1 (MQ=255)
                                                  ttttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCgg                                                                                        >  4:943789/1‑148 (MQ=255)
                                                    ttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGa                                                                                      <  4:225923/148‑1 (MQ=255)
                                                      aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAag                                                                                    <  4:89059/148‑1 (MQ=255)
                                                      aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATAGTGAAACGAAGCGTCTCTCCTCAAAAATGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                   <  2:815645/149‑1 (MQ=255)
                                                         aGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAgaga                                                                                 <  4:30389/148‑1 (MQ=255)
                                                          gCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAAt                                                                               <  1:653331/149‑1 (MQ=255)
                                                           cTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATa                                                                              <  2:102711/149‑1 (MQ=255)
                                                                gTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTc                                                                         >  2:906380/1‑149 (MQ=255)
                                                                gTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTc                                                                         <  1:891073/149‑1 (MQ=255)
                                                                  tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                       <  2:957490/149‑1 (MQ=255)
                                                                  tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcg                                                                        <  4:57602/148‑1 (MQ=255)
                                                                       ggACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAAt                                                                   <  4:640357/148‑1 (MQ=255)
                                                                         aCCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAg                                                                 >  3:363024/1‑148 (MQ=255)
                                                                         aCCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGt                                                                >  1:393292/1‑149 (MQ=255)
                                                                          cccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTa                                                               <  1:1062879/149‑1 (MQ=255)
                                                                           ccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAag                                                                                    >  4:574330/1‑127 (MQ=255)
                                                                           ccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAag                                                                                    <  3:574330/127‑1 (MQ=255)
                                                                           ccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTAt                                                              <  1:1070520/149‑1 (MQ=255)
                                                                            cGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTAt                                                              <  4:58129/148‑1 (MQ=255)
                                                                              cTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATgg                                                            >  4:253777/1‑148 (MQ=255)
                                                                                 gTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCta                                                         <  4:363024/148‑2 (MQ=255)
                                                                                   tATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTa                                                      <  1:819847/149‑1 (MQ=255)
                                                                                       gTGAAACGAAGCGTCTCTCCTCATAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaa                                                  >  2:857533/1‑149 (MQ=255)
                                                                                       gTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaa                                                  >  1:658923/1‑149 (MQ=255)
                                                                                          aaaCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTACAGTg                                                                                                                    >  4:685652/1‑80 (MQ=255)
                                                                                          aaaCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAAc                                                <  4:1035786/148‑1 (MQ=255)
                                                                                             cGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCa                                            >  1:992719/1‑149 (MQ=255)
                                                                                               aaGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGc                                                           >  3:264811/1‑132 (MQ=255)
                                                                                                  cGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAAt                                        <  3:941606/148‑1 (MQ=255)
                                                                                                  cGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATc                                       <  2:220443/149‑1 (MQ=255)
                                                                                                   gTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCa                                      >  2:955820/1‑149 (MQ=255)
                                                                                                     ctctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAt                                     <  4:917063/148‑1 (MQ=255)
                                                                                                       ctcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTc                                   <  4:664680/148‑1 (MQ=255)
                                                                                                          ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGt                                                        <  2:45934/124‑1 (MQ=255)
                                                                                                          ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGt                                                        >  1:45934/1‑124 (MQ=255)
                                                                                                          ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGt                                <  4:529554/148‑1 (MQ=255)
                                                                                                           tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                               >  4:974473/1‑148 (MQ=255)
                                                                                                           tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                              >  1:51128/1‑149 (MQ=255)
                                                                                                           tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                              >  2:219375/1‑149 (MQ=255)
                                                                                                            cAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATaa                                                                             >  4:372158/1‑101 (MQ=255)
                                                                                                            cAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTTTGGCTGTTAGTAAACAGCATAATCATTCGGTCg                              >  4:200582/1‑148 (MQ=255)
                                                                                                            cAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAAc                                                <  3:624067/130‑1 (MQ=255)
                                                                                                            cAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAAc                                                >  4:624067/1‑130 (MQ=255)
                                                                                                            cAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGa                             >  2:39959/1‑149 (MQ=255)
                                                                                                             aaaaCTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGAt                            <  2:658923/149‑1 (MQ=255)
                                                                                                              aaaCTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATa                           <  1:829388/149‑1 (MQ=255)
                                                                                                               aaCTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAc                          <  2:51128/149‑1 (MQ=255)
                                                                                                                    gATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATa                     <  1:39959/149‑1 (MQ=255)
                                                                                                                        tGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAccc                  <  3:979476/148‑1 (MQ=255)
                                                                                                                             tttaaTGAACGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCa            >  1:314098/5‑149 (MQ=39)
                                                                                                                             aaCGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGa                             <  3:542766/132‑1 (MQ=255)
                                                                                                                                gATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGc          <  3:580081/148‑1 (MQ=255)
                                                                                                                                  tGAACGAAGCAGTCGCCTGACTTCTTCGTTCAGTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGt                                                                >  2:632564/1‑92 (MQ=255)
                                                                                                                                  tGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGt                                                                <  1:632564/92‑1 (MQ=255)
                                                                                                                                   gAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCAGTGGGCGCATCGCCGGTAAGAGAATAACGTCACGAATATTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGt       >  3:480159/1‑148 (MQ=255)
                                                                                                                                   gAACGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGtt      >  2:435093/1‑149 (MQ=255)
                                                                                                                                    aaCGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTg     >  1:723731/1‑149 (MQ=255)
                                                                                                                                     aggAAGCAGTCCCCTGAGTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCCGGCAGAGAATAACGTCGCGAATTGTATGGCTGTTAAGAAACAGCAGAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTg     >  4:348611/3‑148 (MQ=255)
                                                                                                                                     aCGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTgg    <  2:314098/149‑1 (MQ=255)
                                                                                                                                     aCGAAGCAGTCAGGCGACTGCTTCGATCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGc                                             >  3:778403/1‑108 (MQ=39)
                                                                                                                                       gAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAagag                                                                                  <  1:204738/69‑1 (MQ=255)
                                                                                                                                       gAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAagag                                                                                  >  2:204738/1‑69 (MQ=255)
                                                                                                                                       gAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg  <  1:857533/149‑1 (MQ=255)
                                                                                                                                       gAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg  >  2:345435/1‑149 (MQ=255)
                                                                                                                                       gAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAc                    <  2:1064347/131‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
TCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCG  >  NC_000913/4353026‑4353309

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