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,2580GA98.0% 666.6 / ‑0.9 197E672K (GAA→AAA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (99/93);  total (101/96)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 199 total aligned reads displayed.

AAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTA  >  NC_000913/4183160‑4183354
                                                                                                  |                                                                                                
  ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAac                                                                                               >  2:1105927/1‑100 (MQ=255)
  ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAac                                                                                               <  2:410023/100‑1 (MQ=255)
    cGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAacac                                                                                             >  2:5019114/1‑100 (MQ=255)
     gAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACg                                                                                            >  1:2685104/1‑100 (MQ=255)
      aaTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGa                                                                                           <  2:3346715/100‑1 (MQ=255)
      aaTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGa                                                                                           <  2:4224326/100‑1 (MQ=255)
       aTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGAt                                                                                          >  1:1878315/1‑100 (MQ=255)
         ccAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGa                                                                                        <  2:1057974/100‑1 (MQ=255)
          cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGAc                                                                                       >  1:3006712/1‑100 (MQ=255)
          cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGAc                                                                                       <  1:3533503/100‑1 (MQ=255)
           aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACg                                                                                      >  1:1920529/1‑100 (MQ=255)
           aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACg                                                                                      >  1:1922953/1‑100 (MQ=255)
            gCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGc                                                                                     <  1:2985552/100‑1 (MQ=255)
             tttGTTCAGCCCGGACCAGGTTGACTACGTGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCTTCCCTGAGCCCATTCCTTAATCACGATGACGcc                                                                                    <  2:2657233/99‑1 (MQ=255)
             cTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGcc                                                                                    >  2:3659410/1‑100 (MQ=255)
                  tCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCg                                                                               >  2:4289350/1‑100 (MQ=255)
                    aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTg                                                                             <  1:3101968/100‑1 (MQ=255)
                    aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTg                                                                             <  2:2182748/100‑1 (MQ=255)
                     gCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGc                                                                            >  2:5041038/1‑100 (MQ=255)
                       tacgtccAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCAt                                                                          >  2:2591468/6‑100 (MQ=255)
                       ngcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCAt                                                                          >  1:4110132/2‑100 (MQ=255)
                         cgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTg                                                                        >  1:834493/1‑100 (MQ=255)
                            ccAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATg                                                                     <  2:4179674/100‑1 (MQ=255)
                              aGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATggg                                                                   <  1:2321366/100‑1 (MQ=255)
                               ggTTGACTACATGGACGTATCCACCCAGCAGGTGTTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGt                                                                  >  1:4080162/1‑100 (MQ=255)
                               ggTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGt                                                                  >  1:3387056/1‑100 (MQ=255)
                                cttGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTg                                                                 <  1:1206859/99‑1 (MQ=255)
                                 ttGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGc                                                                <  2:537127/100‑1 (MQ=255)
                                 ttGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGc                                                                <  2:1675599/100‑1 (MQ=255)
                                    aCTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                             <  1:3988868/100‑1 (MQ=255)
                                     cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                            >  1:4840552/1‑100 (MQ=255)
                                     cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                            <  1:2413541/100‑1 (MQ=255)
                                     cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                            >  2:765820/1‑100 (MQ=255)
                                      tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACa                                                           <  2:3230649/100‑1 (MQ=255)
                                       aCATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACAt                                                          >  2:86020/1‑100 (MQ=255)
                                        cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                         >  2:2107245/1‑100 (MQ=255)
                                        cATGGACGTATCCACCCAGCAGGTGGTAGCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                         <  2:3049519/100‑1 (MQ=255)
                                         aTGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGc                                                        <  1:1113494/100‑1 (MQ=255)
                                            gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAAc                                                     >  2:3972783/1‑100 (MQ=255)
                                              cGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                   >  1:3952086/1‑100 (MQ=255)
                                                tATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGCAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                 >  1:4641985/1‑100 (MQ=255)
                                                tATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCCTTGATGGGTGCGAACATGCAACGTCa                                                 >  1:949695/1‑100 (MQ=255)
                                                tATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                 >  1:371706/1‑100 (MQ=255)
                                                tATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                 >  2:1955885/1‑100 (MQ=255)
                                                 aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                                <  2:3011074/100‑1 (MQ=255)
                                                   ccACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                              <  2:4043114/100‑1 (MQ=255)
                                                   ccACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                              >  2:50714/1‑100 (MQ=255)
                                                   ccACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                              <  1:1168496/100‑1 (MQ=255)
                                                    cACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                             >  2:118758/1‑100 (MQ=255)
                                                    cACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                             >  1:2461015/1‑100 (MQ=255)
                                                     aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCg                                            <  2:3908882/100‑1 (MQ=255)
                                                      cccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgatcgg                                           >  1:2535043/1‑95 (MQ=255)
                                                      cccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGt                                           <  2:380028/100‑1 (MQ=255)
                                                       ccAGCAGGTGGTATCCGTCGGTGCTTCCCAGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCTAACATGCAACGTCAGGCCGtt                                          >  1:4809931/1‑100 (MQ=255)
                                                        cagcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTc                                         <  2:2048418/100‑1 (MQ=255)
                                                         agccgGCGGTATCCGTCGGTGCGTCCCTGACCCCGTTCCTGACACACGATGACGCCCACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTcc                                        <  2:982911/100‑1 (MQ=255)
                                                         agcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTcc                                        <  1:1370021/100‑1 (MQ=255)
                                                          gcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCg                                       >  2:1258361/1‑100 (MQ=255)
                                                          gcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCg                                       <  1:1635611/100‑1 (MQ=255)
                                                             ggtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCAATGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAca                                    >  1:777692/1‑99 (MQ=255)
                                                               tggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActct                                  >  2:2524110/1‑100 (MQ=255)
                                                               tggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAAAGTCAGGCCGTTCCGActct                                  >  2:4545682/1‑100 (MQ=255)
                                                                ggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTg                                 >  2:4217775/1‑100 (MQ=255)
                                                                  tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcg                               >  1:2558206/1‑100 (MQ=255)
                                                                  tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcg                               >  1:4697077/1‑100 (MQ=255)
                                                                   aTCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgc                              >  1:2143231/1‑100 (MQ=255)
                                                                    tCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcg                             <  1:377497/100‑1 (MQ=255)
                                                                      cGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCt                           <  2:777688/100‑1 (MQ=255)
                                                                          ggTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATa                       <  1:5106439/100‑1 (MQ=255)
                                                                          ggTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATa                       <  1:1459887/100‑1 (MQ=255)
                                                                            tGCGTCCCTGATCCCGTTCTTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                     >  1:4402229/1‑100 (MQ=255)
                                                                            tGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                     >  1:106216/1‑100 (MQ=255)
                                                                            ngCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                     >  1:4102712/2‑100 (MQ=255)
                                                                             gCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGc                    <  2:1896169/100‑1 (MQ=255)
                                                                             gCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGc                    <  2:1898572/100‑1 (MQ=255)
                                                                              cGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGcc                   <  2:3563920/100‑1 (MQ=255)
                                                                               gTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                  >  1:1729231/1‑100 (MQ=255)
                                                                                 cccTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt                >  2:178269/1‑100 (MQ=255)
                                                                                 cccTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt                <  1:5148408/100‑1 (MQ=255)
                                                                                  ccTGAAACCGTTCCTGAAACACGATGACTCCAACCGTGCATTGATGGATCCGAACATGCAACGTCAGGCAGCTCCGACTCTGCGCTCTGATAAGCCGCtg               >  2:675693/1‑100 (MQ=255)
                                                                                     gATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACCTCATGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt            <  2:4763749/100‑1 (MQ=255)
                                                                                      aTCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttg           >  1:1123519/1‑100 (MQ=255)
                                                                                       tCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg          <  1:2612082/100‑1 (MQ=255)
                                                                                       tCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg          <  1:2856179/100‑1 (MQ=255)
                                                                                       tCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg          >  1:2820223/1‑100 (MQ=255)
                                                                                        cccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt         >  2:2857419/1‑100 (MQ=255)
                                                                                        cccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt         <  2:135353/100‑1 (MQ=255)
                                                                                        cccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt         <  2:3269313/100‑1 (MQ=255)
                                                                                         ccGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa        >  1:4498648/1‑100 (MQ=255)
                                                                                          cGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTAc       <  2:2421706/100‑1 (MQ=255)
                                                                                          cGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTAc       >  2:370479/1‑100 (MQ=255)
                                                                                           gTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt      >  1:2080423/1‑100 (MQ=255)
                                                                                           gTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt      <  2:4121013/100‑1 (MQ=255)
                                                                                             tCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg    <  1:1404322/100‑1 (MQ=255)
                                                                                             tCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg    >  2:546220/1‑100 (MQ=255)
                                                                                              ccTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt   <  1:3190749/100‑1 (MQ=255)
                                                                                              ccTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt   >  2:1904802/1‑100 (MQ=255)
                                                                                              ccTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt   >  1:3311757/1‑100 (MQ=255)
                                                                                               cTGAAACACGATGACGCCAACCGTGCATTGATGGGTTCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa  >  2:4156882/1‑100 (MQ=255)
                                                                                               cTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa  >  1:1064545/1‑100 (MQ=255)
                                                                                                  |                                                                                                
AAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTA  >  NC_000913/4183160‑4183354

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