Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,234,507 A→T 6.9% intergenic (‑52/+232) 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,5070AT 6.9% 701.0 / 31.8 202intergenic (‑52/+232)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 A (86/102);  new base T (7/7);  total (93/109)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

Only 100 of 205 total aligned reads displayed.

TGCCCCCAACATTGAAAACATAATGCTTATCCAGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAAT  >  NC_000913/3234263‑3234737
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
           ttGAAAACATAATGCTTATCCAGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCTTGTTCCACAGGGCCACCAAATCCGTTTTATGTTGCGTGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTTTCTGTGTGACGGAAGAATAAAc                                                                                                                                                                                               <  2:217644/275‑1 (MQ=255)
                              ctaGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACaa                                                                                                                                                                            <  4:169740/273‑1 (MQ=255)
                                       tGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑Tg                                                                                                                                                                                                                                     >  4:299551/1‑209 (MQ=255)
                                           gTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCaa                                                                                                                                                                                                                                 <  2:299551/209‑1 (MQ=255)
                                              aTCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGTCAATATGATAACATCCGTTGCCAGAT‑TTGCAACGGATGCATTTATCTGTTTGACGGAAGAATAACCGTAGAAAAACGCTAAACAAGAGAATGTAAACGATg                                                                                                                                                            >  2:249195/1‑275 (MQ=255)
                                                                     ttAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     >  2:29248/1‑275 (MQ=255)
                                                                       aaTGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTAGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACTATGGTGACCGGGATTTTTTGTTTTATgg                                                                                                                                   >  2:9516/1‑275 (MQ=255)
                                                                                                ccTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGa                                                                                                                                                                                                                                         >  8:657690/1‑149 (MQ=255)
                                                                                                 cTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCTGTTGCCAGAt                                                                                                                                                                                                                                        <  7:1318619/149‑1 (MQ=255)
                                                                                                 cTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt                                                                                                                                                                                                                                        <  6:791952/149‑1 (MQ=255)
                                                                                                 cTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt                                                                                                                                                                                                                                        <  6:451540/149‑1 (MQ=255)
                                                                                                       aaTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCaa                                                                                                                                                                                                                                 >  6:266619/1‑149 (MQ=255)
                                                                                                        aTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGGGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGcc                                                                                                   >  1:14419/1‑274 (MQ=255)
                                                                                                           gtgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGa                                                                                                                                                                                                                             <  7:841353/149‑1 (MQ=255)
                                                                                                           gtgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGa                                                                                                                                                                                                                             <  7:21894/149‑1 (MQ=255)
                                                                                                            tgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGAt                                                                                                                                                                                                                            >  6:565873/1‑149 (MQ=255)
                                                                                                              ttCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATg                                                                                              >  2:330290/1‑273 (MQ=255)
                                                                                                               tCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACGGATgtt                                                                                                                                                                                                                         <  7:179525/149‑3 (MQ=255)
                                                                                                                ccccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAt                                                                                                                                                                                                                        <  7:1097136/149‑1 (MQ=255)
                                                                                                                  ccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAttt                                                                                                                                                                                                                      <  8:295176/149‑1 (MQ=255)
                                                                                                                      agCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCt                                                                                                                                                                                                                  <  6:718381/149‑1 (MQ=255)
                                                                                                                        ccaccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCATTATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgt                                                                                                                                                                                                                <  8:707263/149‑1 (MQ=255)
                                                                                                                        ccaccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑Tg                                                                                                                                                                                                                                     <  8:1398471/128‑1 (MQ=255)
                                                                                                                        ccaccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgt                                                                                                                                                                                                                <  5:809912/149‑1 (MQ=255)
                                                                                                                          accaACTCCGTTTTATATTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgtgt                                                                                                                                                                                                              >  8:1388570/1‑149 (MQ=255)
                                                                                                                           ccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgtgtg                                                                                                                                                                                                             >  7:1392307/1‑149 (MQ=255)
                                                                                                                           ccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGAAAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAAAACTGGGCAGACATTGCCGGATGCGAAGTTTGCcgcct                                                                               >  4:315338/1‑273 (MQ=255)
                                                                                                                             aaCTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGAc                                                                                                                                                                                                           <  7:677533/149‑1 (MQ=255)
                                                                                                                             aaCTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGAc                                                                                                                                                                                                           >  7:935653/1‑149 (MQ=255)
                                                                                                                             aaCTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGAc                                                                                                                                                                                                           >  8:876485/1‑149 (MQ=255)
                                                                                                                                tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGc                                                                                                                                                                                                                                   >  7:1199451/1‑122 (MQ=255)
                                                                                                                                tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGc                                                                                                                                                                                                                                   <  8:1199451/122‑1 (MQ=255)
                                                                                                                                tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGga                                                                                                                                                                                                        >  6:105691/1‑149 (MQ=255)
                                                                                                                                tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGga                                                                                                                                                                                                        >  6:814640/1‑149 (MQ=255)
                                                                                                                                 ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                                                                                                                       <  6:990067/149‑1 (MQ=255)
                                                                                                                                 ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                                                                                                                       <  5:266619/149‑1 (MQ=255)
                                                                                                                                 ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                                                                                                                       <  5:1428786/149‑1 (MQ=255)
                                                                                                                                  cGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAGATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaag                                                                                                                                                                                                      >  6:974350/1‑149 (MQ=255)
                                                                                                                                   gTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACGGAt                                                                                                                                                                                                                            <  5:845306/126‑1 (MQ=255)
                                                                                                                                   gTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACGGAt                                                                                                                                                                                                                            >  6:845306/1‑126 (MQ=255)
                                                                                                                                    ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaagaa                                                                                                                                                                                                    <  6:666341/149‑1 (MQ=255)
                                                                                                                                    ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaagaa                                                                                                                                                                                                    <  8:192725/149‑1 (MQ=255)
                                                                                                                                     tttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAAt                                                                                                                                                                                                   <  5:672854/149‑1 (MQ=255)
                                                                                                                                        aTGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATaaa                                                                                                                                                                                                <  5:1449576/149‑1 (MQ=255)
                                                                                                                                        aTGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATaaa                                                                                                                                                                                                <  7:554628/149‑1 (MQ=255)
                                                                                                                                          gTTGCGGGTATTTTTCCGCAGCCTCTTCCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATACTATCCGTTGCCCGAT‑TGGCAACGGATGAATTTATCTGTGTTACGGAAGAATTAACGTCGAAAAACGCTAAACAATAGAATGTAGACGATGGTGACAGGGATCTTGTGTTTTATGGAAGCGCTATTCTCGAACTGGGCAGACATTGCCGGAATCTACGCATTCCGCGTTTTAACAGGCCTggc                                                                >  4:287073/1‑274 (MQ=255)
                                                                                                                                             gCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAg                                                                                                                                                                                           <  7:1066045/149‑1 (MQ=255)
                                                                                                                                                 gTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaa                                                                                                                                                                                       >  7:769003/1‑149 (MQ=255)
                                                                                                                                                 gTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaa                                                                                                                                                                                       >  6:637045/1‑149 (MQ=255)
                                                                                                                                                 gTATTTATCCGCAGCATCTTTCAATGCTCTTTGAGTTATCAGGTGCATTCTTCATCTCGACCGTTGACAAATTGGCAATATGTTAACATCCCTGGCAAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaat                                                                                                                                                                                       >  6:200730/1‑148 (MQ=255)
                                                                                                                                                  tATTTTTCCGCAGCATTTTTCAATTCTTTTTGAGTTATCAGGGTCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTTACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTg                                                        <  4:226175/275‑1 (MQ=255)
                                                                                                                                                  tATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaaa                                                                                                                                                                                      <  5:177171/149‑1 (MQ=255)
                                                                                                                                                   aTTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  5:974350/149‑1 (MQ=255)
                                                                                                                                                   aTTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  5:40819/149‑1 (MQ=255)
                                                                                                                                                   aTTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  8:1454282/149‑1 (MQ=255)
                                                                                                                                                   aTTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  6:483806/149‑1 (MQ=255)
                                                                                                                                                    tttttCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACg                                                                                                                                                                                    >  8:1349881/1‑149 (MQ=255)
                                                                                                                                                              gcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          <  8:820657/149‑1 (MQ=255)
                                                                                                                                                               caTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAagag                                                                                                                                                                         >  8:552624/1‑149 (MQ=255)
                                                                                                                                                               caTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCt                                                                             >  2:139777/1‑241 (MQ=255)
                                                                                                                                                                 tCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGATGAATAAACGTAGAAAAACGCTAAACAAGAGaa                                                                                                                                                                       <  8:1248722/149‑1 (MQ=255)
                                                                                                                                                                   tttCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATc                                                                         <  4:139777/241‑1 (MQ=255)
                                                                                                                                                                         tGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          >  8:1075244/1‑138 (MQ=255)
                                                                                                                                                                            tttttGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAt                                                                                                                                                                                                                        >  5:1286349/1‑89 (MQ=255)
                                                                                                                                                                            tttttGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAt                                                                                                                                                                                                                        <  6:1286349/89‑1 (MQ=255)
                                                                                                                                                                              tttGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGt                                                                                                                                                          <  5:647980/149‑1 (MQ=255)
                                                                                                                                                                                 gAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGAc                                                                                                                                                       >  8:685217/1‑149 (MQ=255)
                                                                                                                                                                                      aTCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          >  7:187273/1‑125 (MQ=255)
                                                                                                                                                                                      aTCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGa                                                                                                                                                  <  6:385293/149‑1 (MQ=255)
                                                                                                                                                                                           gTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCTGTTGCCGGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAttttt                                                                                                                                             >  6:864452/1‑149 (MQ=255)
                                                                                                                                                                                            tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt                                                                                                                                            <  6:1114720/149‑1 (MQ=255)
                                                                                                                                                                                              cATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                                                                                                          >  6:1174875/1‑149 (MQ=255)
                                                                                                                                                                                               aTTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtt                                                                                                                                         >  5:1314408/1‑149 (MQ=255)
                                                                                                                                                                                               aTTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtt                                                                                                                                         >  6:1453357/1‑149 (MQ=255)
                                                                                                                                                                                                 tcttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                                                                                                                       >  6:295528/1‑149 (MQ=255)
                                                                                                                                                                                                 tcttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCTGGCCTGGAGTCACGTCCTGAACGTAGGCCGGATAAGGCGTTTCCGCCGCACCCGGCAGGTaa          >  2:36310/1‑274 (MQ=255)
                                                                                                                                                                                                   ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCTGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt                                                                                                                                            <  7:1337056/142‑1 (MQ=255)
                                                                                                                                                                                                   ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCTGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt                                                                                                                                            >  8:1337056/1‑142 (MQ=255)
                                                                                                                                                                                                   ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     <  8:557733/149‑1 (MQ=255)
                                                                                                                                                                                                       tcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGaa                                                                                                                                 <  7:326129/149‑1 (MQ=255)
                                                                                                                                                                                                         aCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          <  8:1028442/106‑1 (MQ=255)
                                                                                                                                                                                                         aCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          >  7:1028442/1‑106 (MQ=255)
                                                                                                                                                                                                         aCGTCCGTTGACAAATTGGCAATATGATAACAGCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAAt  <  4:29248/274‑1 (MQ=255)
                                                                                                                                                                                                            tCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGcgcg                                                                                                                            >  5:619143/1‑149 (MQ=255)
                                                                                                                                                                                                                  gACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCa                                                                                                                      >  7:1241449/1‑149 (MQ=255)
                                                                                                                                                                                                                     aaaTTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGaa                                                                                                                   <  6:1314408/149‑1 (MQ=255)
                                                                                                                                                                                                                         tGGCAATATGATAACATCCGTTGACAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTgg                                                                                                               <  5:1174875/149‑1 (MQ=255)
                                                                                                                                                                                                                                  gATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACAtt                                                                                                      >  8:1008057/1‑149 (MQ=255)
                                                                                                                                                                                                                                    tAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGCCGGAAGAATAAACGTAGAAAATCGCTAAACAAGAGCATGTAGACGATGGTGACAGGGATTTTTTGTCTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGc                                                                                                    <  8:371823/149‑1 (MQ=255)
                                                                                                                                                                                                                                    tAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACTGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGc                                                                                                    >  6:1219989/1‑149 (MQ=255)
                                                                                                                                                                                                                                      aCAGCCGTTGCCAGAT‑TGGCAACGGATGAAATTATCGGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCg                                                                                                  <  8:1282184/149‑1 (MQ=255)
                                                                                                                                                                                                                                        aTCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGa                                                                                                <  7:876485/149‑1 (MQ=255)
                                                                                                                                                                                                                                         tCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGAt                                                                                               >  5:721515/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTgg                                                                                                               >  5:1083923/1‑133 (MQ=255)
                                                                                                                                                                                                                                         tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGc                                                                                                    >  5:582943/1‑144 (MQ=255)
                                                                                                                                                                                                                                         tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGc                                                                                                    <  6:582943/144‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     <  7:460926/110‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGTTGCCA‑ATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     >  8:460926/1‑110 (MQ=255)
                                                                                                                                                                                                                                            gTTGCCAATC‑TGGCAACGGATGAATTGATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCg                                                                                            >  7:1222769/1‑149 (MQ=255)
                                                                                                                                                                                                                                              tGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGAc                                                                                          <  5:565873/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
TGCCCCCAACATTGAAAACATAATGCTTATCCAGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAAT  >  NC_000913/3234263‑3234737

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