Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,183,802 A→G D853G (GAC→GGC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,8020AG98.3% 223.0 / NA 61D853G (GAC→GGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (31/28);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGC  >  NC_000913/4183655‑4183948
                                                                                                                                                   |                                                                                                                                                  
tCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCa                                                                                                                                                  <  2:360076/150‑1 (MQ=255)
  aTGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATcc                                                                                                                                               >  2:425609/1‑151 (MQ=255)
    gCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCg                                                                                                                                             <  2:486020/151‑1 (MQ=255)
     ccGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCg                                                                                                                                             >  2:526255/1‑150 (MQ=255)
         ggAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGt                                                                                                                                        <  1:950956/151‑1 (MQ=255)
            aTGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTggg                                                                                                                                     >  2:702339/1‑151 (MQ=255)
                 tACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGaa                                                                                                                                 <  1:62471/150‑1 (MQ=255)
                     aCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTgc                                                                                                                            >  2:555174/1‑151 (MQ=255)
                        tCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGct                                                                                                                         >  1:891966/1‑151 (MQ=255)
                         cGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGct                                                                                                                         <  2:914863/150‑1 (MQ=255)
                          gAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGctct                                                                                                                       >  2:146956/1‑151 (MQ=255)
                           aaGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCtctc                                                                                                                      >  1:149996/1‑151 (MQ=255)
                            aGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGctct                                                                                                                       <  2:442497/149‑1 (MQ=255)
                             gACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTcc                                                                                                                     >  2:346837/1‑150 (MQ=255)
                              aCTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCaa                                                                                                                   <  1:823639/151‑1 (MQ=255)
                                tCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAc                                                                                                                 >  1:67820/1‑151 (MQ=255)
                                  cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTg                                                                                                               <  1:544292/151‑1 (MQ=255)
                                        cGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATg                                                                                                           <  1:269262/149‑1 (MQ=255)
                                        cGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGaa                                                                                                         <  2:784694/151‑1 (MQ=255)
                                        cGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGa                                                                                                          <  1:643203/150‑1 (MQ=255)
                                               gAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGCGAAGCTGCGCTCTCCAAACTGGATGAATCCGGt                                                                                                   >  2:934962/1‑150 (MQ=255)
                                               gAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCTCTCTCCAAACTGGATGAATCCGGTa                                                                                                  <  2:774024/151‑1 (MQ=255)
                                                 gCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATc                                                                                                <  1:579312/151‑1 (MQ=255)
                                                 gCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATc                                                                                                >  1:111010/1‑151 (MQ=255)
                                                       tgttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTAc                                                                                          >  2:735827/1‑151 (MQ=255)
                                                        gttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTAc                                                                                          <  2:96515/150‑1 (MQ=255)
                                                           cAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTg                                                                                      >  2:333513/1‑151 (MQ=255)
                                                               aaGACCGATTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGCGAAGCTGCGCCCCCCAAACTGGATGAATCCGGTATCGTTTACCTTGGTGc                                                                                  >  2:854390/1‑151 (MQ=255)
                                                                aGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCg                                                                                 >  2:7074/1‑151 (MQ=255)
                                                                      tttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTg                                                                           >  1:120186/1‑151 (MQ=255)
                                                                       ttCAGCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGAGACACCAAGCTGGGTCCGGGAGCGATCACCGCAGGCCCCCCGAACGCGGGCGCAGCAGCGCTCTCCAAACTGGGTGAACCCGGTATCGTTTACATTGGGGCggta                                                                              >  2:302601/1‑145 (MQ=255)
                                                                       ttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTg                                                                           <  1:22054/150‑1 (MQ=255)
                                                                         caccacCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAcc                                                                        <  1:152611/151‑1 (MQ=255)
                                                                            caccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGt                                                                     <  2:216838/151‑1 (MQ=255)
                                                                                aTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCg                                                                 >  2:469566/1‑151 (MQ=255)
                                                                                  ccACATTCAGGAACTGGCGTGTGTGACCCGTGACACCAAGCTGGGTCCGGAAGAGCTCACCGCTGTCATCCCGAACCTGGGCGAAGCTGCGCTCTCCGAACTGGATGAATCCGGTATCGCTTACATTGGTGCGGCAGTGACCGGTGGCGAc                                                               >  1:606853/1‑151 (MQ=255)
                                                                                   cacaTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACa                                                              <  1:766746/151‑1 (MQ=255)
                                                                                     caTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                                                            <  1:295813/151‑1 (MQ=255)
                                                                                       ttCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCATCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGCATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTc                                                           >  2:900102/1‑150 (MQ=255)
                                                                                       ttCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTGTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                                                            <  2:277282/149‑1 (MQ=255)
                                                                                        tCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCGGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtg                                                         >  2:389571/1‑151 (MQ=255)
                                                                                         cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtg                                                         <  2:160571/150‑1 (MQ=255)
                                                                                           ggAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggtt                                                      >  1:278660/1‑151 (MQ=255)
                                                                                           ggAACTGGCGTGTGTGGCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggtt                                                      <  2:399111/151‑1 (MQ=255)
                                                                                                tGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggta                                                  >  1:706247/1‑150 (MQ=255)
                                                                                                 ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggta                                                  >  1:265939/1‑149 (MQ=255)
                                                                                                 ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggta                                                  >  1:438312/1‑149 (MQ=255)
                                                                                                  gCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGTGACATTCTGGTTggtaagg                                               <  2:902599/151‑1 (MQ=255)
                                                                                                         tgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGcc                                        >  2:540023/1‑151 (MQ=255)
                                                                                                          gtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCg                                       >  1:146198/1‑151 (MQ=255)
                                                                                                           tCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGACGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGTCATTCTGGTTGGTAAGGTAACGCCg                                       >  2:672173/1‑150 (MQ=255)
                                                                                                               gTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAgg                                  >  2:113247/1‑151 (MQ=255)
                                                                                                                tGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAgg                                  >  1:222179/1‑150 (MQ=255)
                                                                                                                    acCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGaa                              >  2:739166/1‑150 (MQ=255)
                                                                                                                                  ggAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc                >  1:861391/1‑150 (MQ=255)
                                                                                                                                       agaTCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCag              <  2:149557/147‑1 (MQ=255)
                                                                                                                                        gaTCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAgaagaa         <  1:736073/151‑1 (MQ=255)
                                                                                                                                         aTCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGaaa        <  1:381443/151‑1 (MQ=255)
                                                                                                                                             ccGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAAct    <  1:655880/151‑1 (MQ=255)
                                                                                                                                              cGCTGGCGTCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAGTCTGGTTGGTAAGGTAACGCCGAGAGGTGAAACTCAGCTGACCCCAGAAGAAAAActg   <  2:217121/151‑1 (MQ=255)
                                                                                                                                               gCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActgc  <  2:124256/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGC  >  NC_000913/4183655‑4183948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑