Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,234,508 T→C 5.2% intergenic (‑53/+231) higB ← / ← rlmG mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,234,5080TC 5.2% 310.1 / 5.1 96intergenic (‑53/+231)higB/rlmGmRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (44/47);  new base C (3/2);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

CTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTT  >  NC_000913/3234360‑3234653
                                                                                                                                                    |                                                                                                                                                 
cTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt                                                                                                                                                   >  3:490124/1‑149 (MQ=255)
cTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt                                                                                                                                                   >  3:394453/1‑149 (MQ=255)
  tAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTg                                                                                                                                                 <  1:360944/149‑1 (MQ=255)
      aaTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAtt                                                                                                                                                  >  3:105304/1‑144 (MQ=255)
      aaTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAtt                                                                                                                                                  <  4:105304/144‑1 (MQ=255)
        tCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACg                                                                                                                                           >  3:259113/1‑149 (MQ=255)
          gtgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGa                                                                                                                                         >  2:677883/1‑149 (MQ=255)
           tgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGAt                                                                                                                                        >  3:496729/1‑149 (MQ=255)
               ccccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAt                                                                                                                                    <  2:454638/149‑1 (MQ=255)
               ccccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAt                                                                                                                                    <  2:604695/149‑1 (MQ=255)
               ccccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAt                                                                                                                                    <  4:430682/149‑1 (MQ=255)
                cccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAtt                                                                                                                                   >  3:430741/1‑149 (MQ=255)
                     agCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATATCATCCGTTGCCAGATTGGCAACGCATGAATTTATCt                                                                                                                              <  4:576287/149‑1 (MQ=255)
                     agCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCt                                                                                                                              <  4:653537/149‑1 (MQ=255)
                     agCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCt                                                                                                                              <  3:71547/149‑1 (MQ=255)
                     agCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCt                                                                                                                              <  3:111412/149‑1 (MQ=255)
                       ccaccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCtgt                                                                                                                            <  2:277675/149‑1 (MQ=255)
                       ccaccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCGTTTTGAGTTGTCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCtgt                                                                                                                            <  4:430058/149‑1 (MQ=255)
                           caACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGa                                                                                                                        >  1:577504/1‑149 (MQ=255)
                            aaCTCCGTTTTATGTGGCGGGTATTTTTCCGCAGCATCTTTGAATGCTTTTTGAGTTAGCAGGTGCGTTCTTCATCACGTCCGTTGGCGAATGGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGAc                                                                                                                       >  2:210213/1‑149 (MQ=255)
                               tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTTTCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTCTCTGTGTGACGga                                                                                                                    >  4:580105/1‑149 (MQ=255)
                               tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGga                                                                                                                    >  2:586607/1‑149 (MQ=255)
                               tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGga                                                                                                                    >  4:271870/1‑149 (MQ=255)
                               tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGga                                                                                                                    >  2:304113/1‑149 (MQ=255)
                                ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                                   >  2:567452/1‑149 (MQ=255)
                                ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                                   <  2:260422/149‑1 (MQ=255)
                                   ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGgaagaa                                                                                                                <  2:456853/149‑1 (MQ=255)
                                        tGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATaacc                                                                                                           <  4:490124/149‑3 (MQ=255)
                                        tGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAAc                                                                                                           <  3:496821/149‑1 (MQ=255)
                                        tGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAAc                                                                                                           <  1:373386/149‑1 (MQ=255)
                                           tGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTa                                                                                                        <  2:315041/149‑1 (MQ=255)
                                                        ccGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTaaa                                                                                           <  3:606574/149‑1 (MQ=255)
                                                        ccGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTaaa                                                                                           >  3:8382/1‑149 (MQ=255)
                                                        ccGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTaaa                                                                                           <  1:26229/149‑1 (MQ=255)
                                                         cGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAAc                                                                                          <  3:208867/149‑1 (MQ=255)
                                                          gcagcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACGGAt                                                                                                                                        <  3:29424/102‑1 (MQ=255)
                                                          gcagcaTCTTTCAATGCTTTTTGAGCTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACGGAt                                                                                                                                        >  4:29424/1‑102 (MQ=39)
                                                            agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATCTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAag                                                                                       >  2:681323/1‑149 (MQ=255)
                                                             gcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                      <  2:230083/149‑1 (MQ=255)
                                                             gcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                      >  2:544694/1‑149 (MQ=255)
                                                              caTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATc                                                                                                                               >  4:84761/1‑107 (MQ=255)
                                                              caTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATc                                                                                                                               <  3:84761/107‑1 (MQ=255)
                                                                tCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGaa                                                                                   >  3:348517/1‑149 (MQ=255)
                                                                tCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGaa                                                                                   <  2:224998/149‑1 (MQ=255)
                                                                   ttCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGt                                                                                <  3:520880/149‑1 (MQ=255)
                                                                   ttCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTACCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGt                                                                                <  4:605213/149‑1 (MQ=255)
                                                                    tCAATGCTTTTTGAGTGATCAGGTGCATTCTTCATCACGTCCGTTGGCAAATTGGCAATATGATCACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTa                                                                               <  4:554214/149‑1 (MQ=255)
                                                                      aaTGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGa                                                                             >  3:361923/1‑149 (MQ=255)
                                                                        tGCTTTTTGAGTTATCAGGTGCATTCTTTATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAACGTAGACg                                                                           >  3:586697/1‑149 (MQ=255)
                                                                               tGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGa                                                                    <  3:568027/149‑1 (MQ=255)
                                                                                gAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGAc                                                                   <  4:429390/149‑1 (MQ=255)
                                                                                  gTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAg                                                                 <  2:291797/149‑1 (MQ=255)
                                                                                         tgTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttt                                                          <  1:269383/148‑1 (MQ=255)
                                                                                          gTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAagag                                                                                     <  3:385363/121‑1 (MQ=255)
                                                                                          gTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGAGACGGAAGAATAAACGTAGAAAAACGCTAAACAagag                                                                                     >  4:385363/1‑121 (MQ=255)
                                                                                             cATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                      <  4:259113/149‑1 (MQ=255)
                                                                                             cATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                      <  2:122123/149‑1 (MQ=255)
                                                                                             cATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                      >  1:539027/1‑149 (MQ=255)
                                                                                              agtcttcATCACGTCCGTTGACAAATTGGAAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtt                                                     <  1:567452/147‑1 (MQ=255)
                                                                                              aTTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtt                                                     <  1:156839/149‑1 (MQ=255)
                                                                                               ttcttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                                    >  4:45275/1‑149 (MQ=255)
                                                                                               ttcttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                                    >  4:207613/1‑149 (MQ=255)
                                                                                               ttcttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                                    >  3:212807/1‑149 (MQ=255)
                                                                                                 cttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTa                                                  <  2:577504/149‑1 (MQ=255)
                                                                                                  ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGCCAGGGATTTTTTGTTTTAt                                                 <  2:95659/149‑1 (MQ=255)
                                                                                                  ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                 >  4:430427/1‑149 (MQ=255)
                                                                                                  ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                 >  1:148495/1‑149 (MQ=255)
                                                                                                       caCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAg                                            >  1:522648/1‑149 (MQ=255)
                                                                                                        aCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAgc                                           >  1:416399/1‑149 (MQ=255)
                                                                                                        aCGTCCGTTGAAAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAAGGAAgc                                           >  3:386703/1‑149 (MQ=255)
                                                                                                          gTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAgcgc                                         >  2:517289/1‑149 (MQ=255)
                                                                                                            ccGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGa                                       <  4:496729/149‑1 (MQ=255)
                                                                                                             cGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAt                                      <  3:229274/149‑1 (MQ=255)
                                                                                                              gTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAtt                                     >  3:190503/1‑149 (MQ=255)
                                                                                                               ttGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTCGTTTTATGGAAGGGCGATTc                                    >  1:387618/1‑149 (MQ=255)
                                                                                                                tGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTcc                                   <  4:396894/149‑1 (MQ=255)
                                                                                                                 gACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCa                                  <  3:119216/149‑1 (MQ=255)
                                                                                                                    aaaTTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGaa                               >  1:103005/1‑149 (MQ=255)
                                                                                                                    aaaTTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGaa                               >  3:47466/1‑149 (MQ=255)
                                                                                                                         ggCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTggg                          <  1:586607/149‑1 (MQ=255)
                                                                                                                          gCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGc                         >  3:59472/1‑149 (MQ=255)
                                                                                                                            aaTATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGTAg                       >  3:453854/1‑149 (MQ=255)
                                                                                                                                    aaCATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTgg                           <  3:619647/137‑1 (MQ=255)
                                                                                                                                    aaCATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTgg                           >  4:619647/1‑137 (MQ=255)
                                                                                                                                    aaCATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTgtc               >  1:145879/1‑147 (MQ=255)
                                                                                                                                       aTCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGaa                                                                                   <  1:587959/78‑1 (MQ=38)
                                                                                                                                       aTCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGaa                                                                                   >  2:587959/1‑78 (MQ=38)
                                                                                                                                       aTCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGa            >  3:35336/1‑149 (MQ=255)
                                                                                                                                        tCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGAt           >  2:354866/1‑149 (MQ=255)
                                                                                                                                          cGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACaa                                                                                        >  1:9206/1‑70 (MQ=255)
                                                                                                                                          cGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACaa                                                                                        <  2:9206/70‑1 (MQ=255)
                                                                                                                                          cGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGc         >  1:175656/1‑149 (MQ=255)
                                                                                                                                           gTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCg        <  2:145879/149‑1 (MQ=255)
                                                                                                                                           gTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCg        >  1:425302/1‑149 (MQ=255)
                                                                                                                                           gTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCg        >  3:181069/1‑149 (MQ=255)
                                                                                                                                           gTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCg        <  3:536863/149‑1 (MQ=255)
                                                                                                                                            ttGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGa       <  1:625743/149‑1 (MQ=255)
                                                                                                                                               ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt                                                        >  4:362930/1‑97 (MQ=255)
                                                                                                                                               ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt                                                        <  3:362930/97‑1 (MQ=255)
                                                                                                                                                 aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCtt  >  4:569752/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTT  >  NC_000913/3234360‑3234653

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