Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,183,258 G→A E672K (GAA→AAA)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,2580GA97.0% 531.9 / 3.4 168E672K (GAA→AAA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (76/86);  total (79/89)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 170 total aligned reads displayed.

GGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGT  >  NC_000913/4183162‑4183353
                                                                                                |                                                                                               
ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGaaa                                                                                               <  2:3425343/99‑1 (MQ=255)
ggCGAATCCAGCTTGTTCAACCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGaaa                                                                                               <  1:2749132/99‑1 (MQ=255)
 gCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAac                                                                                              <  1:2768485/99‑1 (MQ=255)
   gAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAacac                                                                                            <  2:1291051/99‑1 (MQ=255)
   gAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAacac                                                                                            >  1:2890443/1‑99 (MQ=255)
    aaTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACg                                                                                           <  1:1838663/99‑1 (MQ=255)
     aTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGa                                                                                          <  1:4019843/99‑1 (MQ=255)
      tCCAGCTTGTTCAGCCGCGACCATGTTGACTACATTGACGTATCCACCCAGCAGGTGCTATCCGTCGGTGCGTCCCTGGTCCCGTTCCTGAAACACGAt                                                                                         <  1:4981884/99‑1 (MQ=255)
        cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGa                                                                                       <  1:2972211/99‑1 (MQ=255)
         aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGAc                                                                                      >  2:53617/1‑99 (MQ=255)
         aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGATGCGTCCCTGATCCCGTTCCTGAAACACGATGAc                                                                                      >  1:1159800/1‑99 (MQ=255)
          gCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACg                                                                                     >  1:2098466/1‑99 (MQ=255)
           cTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGc                                                                                    <  2:1859360/99‑1 (MQ=255)
               ttCAGCCGGGAGCAGGGTGACGACATGGACGTTGCCAGCAAGCAGGTGGTATCCGTCGGTTCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAAc                                                                                <  2:3981936/99‑1 (MQ=255)
               ttCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGGGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAAc                                                                                >  1:2407560/1‑99 (MQ=255)
                tCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAAcc                                                                               >  2:3348456/1‑99 (MQ=255)
                 cAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCg                                                                              >  2:2520233/1‑99 (MQ=255)
                  atccGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGt                                                                             <  1:4448948/97‑1 (MQ=255)
                  aGCCGCGACCAGGTTGACTACATGGTCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGGCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGt                                                                             <  1:589154/99‑1 (MQ=255)
                  aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGt                                                                             <  2:209293/99‑1 (MQ=255)
                   gCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTg                                                                            <  2:3441/99‑1 (MQ=255)
                    ccGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGc                                                                           <  2:4909126/99‑1 (MQ=255)
                      gcgACCAGGTGGACTACATGGACGTATCCACCCAGCAGGTGGTATGCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCAt                                                                         >  2:4790708/1‑99 (MQ=255)
                       cgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGCTCCCGTTCCTGAAACACGATGACGCCAACCGTGCttt                                                                        <  2:304845/99‑1 (MQ=255)
                        gACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTg                                                                       <  1:614900/99‑1 (MQ=255)
                        gACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTg                                                                       >  2:1324447/1‑99 (MQ=255)
                         aCCAGGTTGACTACATGTACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGa                                                                      <  1:2708975/99‑1 (MQ=255)
                         aCCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGa                                                                      <  1:1496718/99‑1 (MQ=255)
                         aCCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGa                                                                      <  1:4737973/99‑1 (MQ=255)
                          ccAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGAt                                                                     <  1:159593/99‑1 (MQ=255)
                            aGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATgg                                                                   >  1:150409/1‑99 (MQ=255)
                             ggTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATggg                                                                  >  1:4195121/1‑99 (MQ=255)
                              gTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGt                                                                 <  2:1467817/99‑1 (MQ=255)
                               ttGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTg                                                                <  2:618625/99‑1 (MQ=255)
                                tGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGc                                                               <  1:129659/99‑1 (MQ=255)
                                   cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                            <  2:298728/99‑1 (MQ=255)
                                   cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                            >  1:1939079/1‑99 (MQ=255)
                                    tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATTACGCCAACCGTGCATTGATGGGTGCGaat                                                           >  1:5279572/1‑98 (MQ=255)
                                    tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                           >  1:2647918/1‑99 (MQ=255)
                                    tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                           <  1:2110560/99‑1 (MQ=255)
                                    tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                           >  1:2974291/1‑99 (MQ=255)
                                     aCATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACa                                                          <  2:1977633/99‑1 (MQ=255)
                                      cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACAt                                                         <  1:97142/99‑1 (MQ=255)
                                      cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACAt                                                         >  1:2015113/1‑99 (MQ=255)
                                       aTGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                        >  1:4630132/1‑99 (MQ=255)
                                         ggACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCa                                                      >  1:3693132/1‑99 (MQ=255)
                                         ggACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCa                                                      >  2:5014876/1‑99 (MQ=255)
                                         ggACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCa                                                      >  1:5141433/1‑99 (MQ=255)
                                          gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                     <  1:3887246/99‑1 (MQ=255)
                                          gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                     <  1:226320/99‑1 (MQ=255)
                                           aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAAc                                                    >  2:2197141/1‑99 (MQ=255)
                                             gTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                  <  2:2643553/99‑1 (MQ=255)
                                             gTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                  <  2:2926573/99‑1 (MQ=255)
                                              tATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTc                                                 <  2:2603522/99‑1 (MQ=255)
                                               aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                <  1:5187887/99‑1 (MQ=255)
                                               aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATAACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                >  1:1373505/1‑99 (MQ=255)
                                                tCCACCCAGGAGGTGGTATACGTCGGTTAGACCCTGTTCCCGTTCCTGAAACAAGATGCCGCCAAGCGTGCATTGATGGGTGCGAACATGCGACGTCAg                                               <  2:3917019/99‑1 (MQ=255)
                                                tCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGCCGGCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                               >  2:4798598/1‑99 (MQ=255)
                                                tCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                               <  1:2136687/99‑1 (MQ=255)
                                                   aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                            >  1:2063239/1‑99 (MQ=255)
                                                   aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                            <  1:342803/99‑1 (MQ=255)
                                                        gcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTcc                                       <  1:4633984/99‑1 (MQ=255)
                                                        gcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGACAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTcc                                       >  2:4472080/1‑99 (MQ=255)
                                                         cagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTTCTGAAACACGATGCCGCCAACCGTTCATTGATGGGTGCGAATATGCAACGTCAGGCCGTTCCg                                      <  2:4163568/99‑1 (MQ=255)
                                                          agGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                     <  2:113419/99‑1 (MQ=255)
                                                           ggtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAc                                    <  1:1412263/99‑1 (MQ=255)
                                                             tggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActc                                  >  2:4348536/1‑99 (MQ=255)
                                                             tggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActc                                  >  1:2557023/1‑99 (MQ=255)
                                                             tggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATGGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActc                                  >  1:5044758/1‑99 (MQ=255)
                                                              ggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTACATTAATAGATACGAACATACAACGTCAGGCCGTTCCGActct                                 >  1:3047156/1‑99 (MQ=255)
                                                                tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgc                               <  2:716008/99‑1 (MQ=255)
                                                                 aTCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcg                              <  1:5007280/99‑1 (MQ=255)
                                                                 aTCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcg                              <  1:2236769/99‑1 (MQ=255)
                                                                    tgTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACTGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                           <  1:1785015/98‑1 (MQ=255)
                                                                    cGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                           <  2:756655/99‑1 (MQ=255)
                                                                    cGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                           <  1:3910981/99‑1 (MQ=255)
                                                                    cGTCGGTGCGTCCCTGATCCCATTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                           <  1:1035524/99‑1 (MQ=255)
                                                                    cGTCGGTGCGTCCCTGATCCCATTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                           <  1:2331063/99‑1 (MQ=255)
                                                                     gTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCt                          <  1:1540206/99‑1 (MQ=255)
                                                                     gTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCt                          >  2:282974/1‑99 (MQ=255)
                                                                      tCGGTGCGTCCCTGCTCCCGTTCCTGAAAGACGATGACGCCAACTGTGCATTGAAGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTg                         >  2:5075466/1‑99 (MQ=255)
                                                                         gTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATa                      >  1:2250857/1‑99 (MQ=255)
                                                                         gTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATa                      >  2:1966388/1‑99 (MQ=255)
                                                                          tGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATaa                     <  1:3263349/99‑1 (MQ=255)
                                                                             gTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGcc                  >  1:3398887/1‑99 (MQ=255)
                                                                              tCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                 >  2:4161781/1‑99 (MQ=255)
                                                                              tCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                 >  1:3747998/1‑99 (MQ=255)
                                                                                 cTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCAAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtg              >  1:422549/1‑99 (MQ=255)
                                                                                  tGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg             >  2:1487410/1‑99 (MQ=255)
                                                                                    aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt           <  1:2003763/99‑1 (MQ=255)
                                                                                    aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt           <  2:1663883/99‑1 (MQ=255)
                                                                                    aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt           <  2:1814222/99‑1 (MQ=255)
                                                                                      cccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg         <  2:1232051/99‑1 (MQ=255)
                                                                                       ccGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt        <  1:577084/99‑1 (MQ=255)
                                                                                       ccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt        >  2:3162497/1‑99 (MQ=255)
                                                                                       ccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt        >  2:4861546/1‑99 (MQ=255)
                                                                                           tCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg    <  1:3339044/99‑1 (MQ=255)
                                                                                            ccTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg   <  1:1747238/99‑1 (MQ=255)
                                                                                             cTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt  <  2:3453412/99‑1 (MQ=255)
                                                                                             cTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt  >  1:3282141/1‑99 (MQ=255)
                                                                                                |                                                                                               
GGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGT  >  NC_000913/4183162‑4183353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: