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.5% 515.3 / 0.4 158E672K (GAA→AAA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (69/85);  total (71/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 162 total aligned reads displayed.

GGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGG  >  NC_000913/4183162‑4183352
                                                                                                |                                                                                              
ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGaaa                                                                                              <  2:1841612/99‑1 (MQ=255)
ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGaaa                                                                                              <  1:1600616/99‑1 (MQ=255)
   gAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAacac                                                                                           <  1:892321/99‑1 (MQ=255)
   gAATCCAGCTTGTTCAGCCGCGACCAGCTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAacac                                                                                           <  1:3578157/99‑1 (MQ=255)
      tCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGAt                                                                                        <  2:1010921/99‑1 (MQ=255)
      tCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGAt                                                                                        <  1:1559359/99‑1 (MQ=255)
        cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCTTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGa                                                                                      >  2:3585700/1‑99 (MQ=255)
        cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGa                                                                                      <  1:3941872/99‑1 (MQ=255)
        cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGAGCCCGTTCCTGAAACACGATGa                                                                                      <  2:76700/99‑1 (MQ=255)
         aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGAc                                                                                     >  1:3119720/1‑99 (MQ=255)
            ttGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGcc                                                                                  <  1:1681104/99‑1 (MQ=255)
              gTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCaa                                                                                >  1:1093823/1‑99 (MQ=255)
              gTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCaa                                                                                >  1:1093847/1‑99 (MQ=255)
                tCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAAcc                                                                              <  1:1959666/99‑1 (MQ=255)
                   gCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTTGGTGCGTCCCTGATCCCGTTCATGAAACACGATGACGCCAACCGTg                                                                           <  1:3690739/99‑1 (MQ=255)
                    ccGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGc                                                                          <  1:2318517/99‑1 (MQ=255)
                     cgcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCa                                                                         <  1:2781301/99‑1 (MQ=255)
                      gcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCAt                                                                        <  1:1341749/99‑1 (MQ=255)
                       cgACCAGGTTGACTACATGGACGTAGCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCAtt                                                                       <  1:4030077/99‑1 (MQ=255)
                        gACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTg                                                                      >  1:317833/1‑99 (MQ=255)
                          ccAGGTTGTCTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGAt                                                                    <  2:960298/99‑1 (MQ=255)
                          ccAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTAGGTGCGGCCCTGATCCCGTTCCTGAAACACGATGACGCCCACCGTGCATTGAt                                                                    <  2:3794287/99‑1 (MQ=255)
                           cAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATg                                                                   >  1:2869763/1‑99 (MQ=255)
                            aGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATgg                                                                  <  1:2986975/99‑1 (MQ=255)
                             ggTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATggg                                                                 >  1:1802651/1‑99 (MQ=255)
                              gTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGt                                                                <  2:3609802/99‑1 (MQ=255)
                                   ccccATAGACGTCCCCCCCCAGCGGGAGGTATCCGTCGGTGCGTCTCTGATACCGTACCTGAAACACGATGACGACAAGCGTGCATTGATTGGGGCGac                                                           <  2:2439878/96‑2 (MQ=255)
                                   cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                           >  1:937143/1‑99 (MQ=255)
                                   cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                           <  1:4066587/99‑1 (MQ=255)
                                      cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACAt                                                        <  2:439941/99‑1 (MQ=255)
                                      cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACAt                                                        <  2:2960879/99‑1 (MQ=255)
                                        tGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGc                                                      <  1:1166300/99‑1 (MQ=255)
                                         ggACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCa                                                     >  1:2710933/1‑99 (MQ=255)
                                          gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                    <  1:249367/99‑1 (MQ=255)
                                          gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                    <  1:2158339/99‑1 (MQ=255)
                                          gACCTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                    <  2:3139287/99‑1 (MQ=255)
                                           aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGCCGCCAACCGTGCATTGATGGGTGCGAACATGCCAc                                                   >  2:3851781/1‑99 (MQ=255)
                                           aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAAc                                                   >  2:432298/1‑99 (MQ=255)
                                           aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTAATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAAc                                                   >  2:1124385/1‑99 (MQ=255)
                                             gTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                 >  1:3984542/1‑99 (MQ=255)
                                               aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAACCACGATGACTCCAACCGTCCATTGATGGGTGCGAACATGCAACGTCa                                               >  2:2293996/1‑99 (MQ=255)
                                               aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                               >  1:3913760/1‑99 (MQ=255)
                                                 ccACCCAGCAGGTGGTATCCGTCGGTGCGTCCATGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                             <  2:1287519/99‑1 (MQ=255)
                                                  cACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                            >  1:3697197/1‑99 (MQ=255)
                                                   aCCCAGCAGGTGGTATCCGTGGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCTTTGATGGGTGCGAACATGCAACGTCAGGcc                                           >  2:2003950/1‑99 (MQ=255)
                                                   aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                           >  1:261776/1‑99 (MQ=255)
                                                   aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                           >  1:2393353/1‑99 (MQ=255)
                                                    cccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCg                                          >  1:2061289/1‑99 (MQ=255)
                                                     ccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGt                                         >  1:53635/1‑99 (MQ=255)
                                                       gctgggtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTc                                       >  2:3870565/5‑99 (MQ=255)
                                                       agcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTc                                       <  2:3286125/99‑1 (MQ=255)
                                                        gcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGTGAACATGCAACGTCAGGCCGTTcc                                      <  1:459395/99‑1 (MQ=255)
                                                          agGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                    <  2:3076464/99‑1 (MQ=255)
                                                          agGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                    <  2:1584578/99‑1 (MQ=255)
                                                          agGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                    <  1:975004/99‑1 (MQ=255)
                                                          agGTCGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                    <  1:2260187/99‑1 (MQ=255)
                                                            gtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                  >  2:2158874/1‑99 (MQ=255)
                                                            gtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                  >  2:2411210/1‑99 (MQ=255)
                                                            gtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                  >  2:2801091/1‑99 (MQ=255)
                                                            gtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                  >  1:1074447/1‑99 (MQ=255)
                                                             tggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActc                                 >  2:51762/1‑99 (MQ=255)
                                                             tggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActc                                 <  2:883869/99‑1 (MQ=255)
                                                               gtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTg                               >  1:3879047/1‑99 (MQ=255)
                                                               gtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTg                               >  1:3446065/1‑99 (MQ=255)
                                                                tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgc                              >  2:1667234/1‑99 (MQ=255)
                                                                tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgc                              >  2:1668520/1‑99 (MQ=255)
                                                                 aTCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcg                             >  1:2433493/1‑99 (MQ=255)
                                                                    cGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                          <  1:429461/99‑1 (MQ=255)
                                                                      tCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGCACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTg                        <  2:3966544/99‑1 (MQ=255)
                                                                        ggTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGAt                      >  1:1349355/1‑99 (MQ=255)
                                                                        ggTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGAt                      <  1:537220/99‑1 (MQ=255)
                                                                         gTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATa                     >  2:481132/1‑99 (MQ=255)
                                                                          tGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATaa                    <  1:504420/99‑1 (MQ=255)
                                                                          tGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATaa                    >  1:2210429/1‑99 (MQ=255)
                                                                          tGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATaa                    <  1:287766/99‑1 (MQ=255)
                                                                          tGCGTCCCTGATCCCGTGCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGACGTTCCGACTCTGCGCGCTGATaa                    <  2:915115/99‑1 (MQ=255)
                                                                              tCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                >  1:1077956/1‑99 (MQ=255)
                                                                               cccGGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc               >  2:3428414/1‑99 (MQ=255)
                                                                                ccTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt              >  1:3288327/1‑99 (MQ=255)
                                                                                ccTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt              >  2:1639745/1‑99 (MQ=255)
                                                                                ccTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt              >  2:3142041/1‑99 (MQ=255)
                                                                                ccTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt              >  1:150035/1‑99 (MQ=255)
                                                                                 cTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtg             >  2:1988124/1‑99 (MQ=255)
                                                                                  tGATGCCGTCCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg            <  2:2809739/99‑1 (MQ=255)
                                                                                  tGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg            >  1:1425301/1‑99 (MQ=255)
                                                                                   gATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt           >  1:3332436/1‑99 (MQ=255)
                                                                                   gATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt           >  1:676842/1‑99 (MQ=255)
                                                                                    aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt          <  1:1878442/99‑1 (MQ=255)
                                                                                    aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt          >  2:259423/1‑99 (MQ=255)
                                                                                    aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt          <  2:2843955/99‑1 (MQ=255)
                                                                                    aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt          <  2:3385267/99‑1 (MQ=255)
                                                                                     tCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttg         <  1:1275738/99‑1 (MQ=255)
                                                                                     tCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttg         >  1:3041370/1‑99 (MQ=255)
                                                                                       ccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt       <  2:1504217/99‑1 (MQ=255)
                                                                                        cGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa      <  1:1935860/99‑1 (MQ=255)
                                                                                         gTTCCTGAAACACGATGACGCCAACCGTGCACTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACACTGCGCGCTGAGAAGCCGCTGGTTGGAAc     <  2:2237923/99‑1 (MQ=255)
                                                                                          ttCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt    <  1:3812020/99‑1 (MQ=255)
                                                                                           tCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg   <  1:322822/99‑1 (MQ=255)
                                                                                            ccTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg  <  2:1331321/99‑1 (MQ=255)
                                                                                            ccTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCCGATAAGCCGCTGGTTGGTACTgg  <  1:3853685/99‑1 (MQ=255)
                                                                                                |                                                                                              
GGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGG  >  NC_000913/4183162‑4183352

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