Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,182,574 G→T 100% D444Y (GAT→TAT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,5740GT100.0% 487.6 / NA 132D444Y (GAT→TAT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (69/63);  total (69/63)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 135 total aligned reads displayed.

GACTTGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATAT  >  NC_000913/4182430‑4182719
                                                                                                                                                |                                                                                                                                                 
 aCTTGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATc                                                                                                                                              <  1:651483/149‑1 (MQ=255)
 aCTTGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGACGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATc                                                                                                                                              <  3:646959/149‑1 (MQ=255)
   ttGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGa                                                                                                                                            >  2:567504/1‑149 (MQ=255)
   ttGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGa                                                                                                                                            >  2:500872/1‑149 (MQ=255)
   ttGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGa                                                                                                                                            >  1:244395/1‑149 (MQ=255)
   ttGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGa                                                                                                                                            >  3:351280/1‑149 (MQ=255)
   ttGTCTGCGGATGGTCGTATGAAGTTAAGCCGTTGTCTGCTGCACGAAGAAATCGAAGGTTCCGGTATCCTGAGCAGTGAAGACATCATTGATGATATGTAAGAGATCATCGATATCCGTAACGGTAAAGGCGCAGTCGATTATATCGa                                                                                                                                            >  2:827603/1‑149 (MQ=255)
    tGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGac                                                                                                                                           <  4:677868/149‑1 (MQ=255)
    tGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGac                                                                                                                                           <  1:53998/149‑1 (MQ=255)
         gCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCt                                                                                                                                      <  4:125912/149‑1 (MQ=255)
                gTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAAc                                                                                                                               >  4:303431/1‑149 (MQ=255)
                 tCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCTGGTATCCTGAGCTAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAAcc                                                                                                                              >  4:680017/1‑149 (MQ=255)
                   gTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACcgt                                                                                                                            >  1:347581/1‑149 (MQ=255)
                    tATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACcgtc                                                                                                                           <  4:842092/149‑1 (MQ=255)
                      tGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACcgtcgt                                                                                                                         >  4:899351/1‑149 (MQ=255)
                      tGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACcgtcgt                                                                                                                         >  3:834495/1‑149 (MQ=255)
                      tGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACcgtcgt                                                                                                                         >  1:36154/1‑149 (MQ=255)
                        aaGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAt                                                                                                                       <  2:786740/149‑1 (MQ=255)
                          gTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtcc                                                                                                                     >  3:420399/1‑149 (MQ=255)
                           ttCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccg                                                                                                                    >  4:595216/1‑149 (MQ=255)
                            tCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccgt                                                                                                                   <  3:616129/149‑1 (MQ=255)
                              aaCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                 >  2:821287/1‑149 (MQ=255)
                              aaCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                 >  2:5902/1‑149 (MQ=255)
                              aaCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                 >  4:156946/1‑149 (MQ=255)
                              aaCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                 >  3:46744/1‑149 (MQ=255)
                              aaCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                 <  1:513676/149‑1 (MQ=255)
                                  gTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATccgttccgtt                                                                                                             <  2:345240/149‑1 (MQ=255)
                                   ttCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                                                                                                            <  3:432378/149‑1 (MQ=255)
                                     ctctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGc                                                                                                          >  4:638697/1‑149 (MQ=255)
                                      tctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                                                                                                         >  4:817766/1‑149 (MQ=255)
                                      tctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                                                                                                         >  2:413359/1‑149 (MQ=255)
                                           tgcGCGAAGAAATAGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATg                                                                                                    <  1:500872/149‑1 (MQ=255)
                                             cgcgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGc                                                                                                  <  4:31616/149‑1 (MQ=255)
                                                gaagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATCTCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                               >  2:628061/1‑149 (MQ=255)
                                                gaagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                               >  1:163447/1‑149 (MQ=255)
                                                gaagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                               >  3:649915/1‑149 (MQ=255)
                                                gaagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGCCGAAGTCATTTATATCGACAACCTCGGCAACCGTGGTATCGTTTCCGTTGGGGAAATGGCGGa                                                                                               >  2:73940/1‑149 (MQ=255)
                                                   gaaATCGAAGGTTCCGGTATCCCGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaaa                                                                                            <  3:567766/149‑1 (MQ=255)
                                                    aaaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                           <  4:55521/149‑1 (MQ=255)
                                                         gAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                           <  4:386792/144‑1 (MQ=255)
                                                          aaGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaaa                                                                                                      <  2:200351/132‑1 (MQ=255)
                                                          aaGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaaa                                                                                                      >  1:200351/1‑132 (MQ=255)
                                                              atCCGGTATCCTGAGCCAAGACGACATCATTGCTGTTATGAAAAAGCGCTTCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCGCGGCAACCGTCTTAGACGTTCCGTTAGAGCAATGGCGGAAAACAAGTTCcgcg                                                                                 >  2:90106/2‑149 (MQ=255)
                                                                ccGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGtt                                                                               <  2:244395/149‑1 (MQ=255)
                                                                 cGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTg                                                                              <  4:89791/149‑1 (MQ=255)
                                                                 cGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTg                                                                              <  4:834495/149‑1 (MQ=255)
                                                                       ccTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTgg                                                                        <  1:793712/149‑1 (MQ=255)
                                                                         tGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTa                                                                      >  4:440105/1‑149 (MQ=255)
                                                                          gAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTAc                                                                     >  2:181322/1‑149 (MQ=255)
                                                                            gCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgt                                                                   >  3:873106/1‑149 (MQ=255)
                                                                              aaaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgtgt                                                                 <  3:680017/149‑1 (MQ=255)
                                                                                   cgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGc                                                            <  2:637057/149‑1 (MQ=255)
                                                                                    gacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCg                                                           >  2:300295/1‑149 (MQ=255)
                                                                                     acATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGt                                                          <  2:343305/149‑1 (MQ=255)
                                                                                           ttGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTg                                                    >  4:690970/1‑149 (MQ=255)
                                                                                            tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTTGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGa                                                   <  3:837444/149‑1 (MQ=255)
                                                                                             gATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAtccgttccgt                                                                                                              >  2:532835/1‑89 (MQ=255)
                                                                                             gATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGaa                                                  <  3:899351/149‑1 (MQ=255)
                                                                                             gATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGaa                                                  <  2:295467/149‑1 (MQ=255)
                                                                                                gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAaga                                               <  4:179532/149‑1 (MQ=255)
                                                                                                gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAaga                                               >  1:84992/1‑149 (MQ=255)
                                                                                                gTTATGAAAAAGCTCATCGATAACCGTAACGGTAAAGGAGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTAGTACGTGTAGAGCGTGCGGTGAAaga                                               <  3:346904/149‑1 (MQ=255)
                                                                                                 ttATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAagag                                              >  4:418765/1‑149 (MQ=255)
                                                                                                 ttATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAagag                                              >  3:97043/1‑149 (MQ=255)
                                                                                                 ttATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAagag                                              >  1:343134/1‑149 (MQ=255)
                                                                                                  tATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGc                                             <  1:73940/149‑1 (MQ=255)
                                                                                                     gAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtc                                          >  3:780772/1‑149 (MQ=255)
                                                                                                      aaaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtct                                         <  3:156946/149‑1 (MQ=255)
                                                                                                            cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTctct                                   <  2:163447/149‑1 (MQ=255)
                                                                                                             tCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTg                                  <  2:83060/149‑1 (MQ=255)
                                                                                                             tCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTg                                  >  4:461119/1‑149 (MQ=255)
                                                                                                              cATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTgg                                 <  4:420399/149‑1 (MQ=255)
                                                                                                               aTCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTggg                                <  1:129809/149‑1 (MQ=255)
                                                                                                                 cGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCg                              >  4:778995/1‑149 (MQ=255)
                                                                                                                  gATATCCGTAACGGTAAAGGCGAAGTTGATTATATCGACCATCTCGGCAACCGTCGTATCCGTTCCGTTCGCGAAATGGCGGAAAACCAGATCCGCGTTGGCCTGGTACGTGTAGGGCATGCGGTGAAAGAGCGTCTGTCTCTGGGcca                             <  2:487740/149‑3 (MQ=255)
                                                                                                                     atCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCt                          >  4:210133/1‑149 (MQ=255)
                                                                                                                     atCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCt                          >  2:927811/1‑149 (MQ=255)
                                                                                                                     atCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCt                          >  3:628027/1‑149 (MQ=255)
                                                                                                                        cGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTATCTCTGGGCGATCTGGa                       <  2:3191/149‑1 (MQ=255)
                                                                                                                         gTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                      >  1:659600/1‑149 (MQ=255)
                                                                                                                         gTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                      >  4:76596/1‑149 (MQ=255)
                                                                                                                          tAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATa                     >  3:599251/1‑149 (MQ=255)
                                                                                                                           aaCGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtc                                          >  1:567101/1‑127 (MQ=255)
                                                                                                                           aaCGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtc                                          <  2:567101/127‑1 (MQ=255)
                                                                                                                             cGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAccc                  <  3:778995/149‑1 (MQ=255)
                                                                                                                             cGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAccc                  <  3:281559/149‑1 (MQ=255)
                                                                                                                             cGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAccc                  <  2:726935/149‑1 (MQ=255)
                                                                                                                                tAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGCAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCCGTCTCTGGGCGAt                            <  4:598355/136‑1 (MQ=255)
                                                                                                                                tAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGAt                            >  3:598355/1‑136 (MQ=255)
                                                                                                                                tAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa               >  1:147155/1‑149 (MQ=255)
                                                                                                                                tAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa               >  1:192095/1‑149 (MQ=255)
                                                                                                                                tAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa               >  1:682703/1‑149 (MQ=255)
                                                                                                                                   aGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGc            <  2:580724/149‑1 (MQ=255)
                                                                                                                                     gCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCca          >  4:767721/1‑149 (MQ=255)
                                                                                                                                     gCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCca          <  2:36154/149‑1 (MQ=255)
                                                                                                                                        aaGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAg       <  3:387910/149‑1 (MQ=255)
                                                                                                                                          gTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGa     >  1:447327/1‑149 (MQ=255)
                                                                                                                                           tCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGat    <  1:58226/149‑1 (MQ=255)
                                                                                                                                             gattatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcac         >  1:779107/1‑142 (MQ=255)
                                                                                                                                             gattatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcac         <  2:779107/142‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GACTTGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATAT  >  NC_000913/4182430‑4182719

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: