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.7% 288.5 / ‑4.0 78D853G (GAC→GGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (40/37);  total (41/37)
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.28e-01

CATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTG  >  NC_000913/4183656‑4183947
                                                                                                                                                  |                                                                                                                                                 
cATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCAt                                                                                                                                                <  2:606676/150‑1 (MQ=255)
      gTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGaa                                                                                                                                          >  1:1027506/1‑150 (MQ=255)
      gTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAAc                                                                                                                                         <  1:488102/151‑1 (MQ=255)
       tGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAAc                                                                                                                                         >  2:671947/1‑150 (MQ=255)
               ttACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGa                                                                                                                                 >  2:908027/1‑150 (MQ=255)
                tACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGa                                                                                                                                 >  1:1147421/1‑149 (MQ=255)
                 aCAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCTGAACGTGGGTGaa                                                                                                                                <  2:38678/149‑1 (MQ=255)
                    aCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCt                                                                                                                             >  1:803803/1‑149 (MQ=255)
                     cTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTgcg                                                                                                                          >  1:247046/1‑151 (MQ=255)
                     cTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTgcg                                                                                                                          >  1:713754/1‑151 (MQ=255)
                       tCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGct                                                                                                                        >  2:252626/1‑151 (MQ=255)
                         gAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGctct                                                                                                                      >  1:103838/1‑151 (MQ=255)
                             aCTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCaa                                                                                                                  <  2:916734/151‑1 (MQ=255)
                             aCTCCATCCTCGTATCCGAGCGGGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCaa                                                                                                                  <  2:1091241/151‑1 (MQ=255)
                               tCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAc                                                                                                                >  2:994797/1‑151 (MQ=255)
                                ccATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAc                                                                                                                <  2:495750/150‑1 (MQ=255)
                                ccATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACt                                                                                                               <  2:575157/151‑1 (MQ=255)
                                 cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTg                                                                                                              >  1:88925/1‑151 (MQ=255)
                                 cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTg                                                                                                              >  1:474565/1‑151 (MQ=255)
                                 cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTg                                                                                                              <  2:415223/151‑1 (MQ=255)
                                  aTCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTgg                                                                                                             >  1:621333/1‑151 (MQ=255)
                                     cTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCAGCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGGAGCTGCGCTCTCCAAACTGGa                                                                                                            <  2:1026105/149‑1 (MQ=255)
                                      tCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATg                                                                                                          >  1:843926/1‑150 (MQ=255)
                                       cGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGa                                                                                                         >  1:1038386/1‑150 (MQ=255)
                                        gTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAAt                                                                                                       >  1:685826/1‑151 (MQ=255)
                                         tATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATc                                                                                                      >  1:786774/1‑151 (MQ=255)
                                            ccGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCgg                                                                                                   >  2:322383/1‑151 (MQ=255)
                                             cGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCgg                                                                                                   >  2:567242/1‑150 (MQ=255)
                                             cGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGt                                                                                                  <  1:700905/151‑1 (MQ=255)
                                                 cGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCg                                                                                              <  1:841745/151‑1 (MQ=255)
                                                   tgttgtTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGtt                                                                                            >  2:706206/1‑151 (MQ=255)
                                                        ttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAt                                                                                       >  1:410218/1‑151 (MQ=255)
                                                         tCAGGAAGACCGTTTCACCACCATCCACATCCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAtt                                                                                      >  1:256623/1‑151 (MQ=255)
                                                          cAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTg                                                                                     >  1:857754/1‑151 (MQ=255)
                                                          cAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTg                                                                                     <  1:381738/151‑1 (MQ=255)
                                                          cAGGAAGACCGTTTAACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTg                                                                                     >  2:354462/1‑151 (MQ=255)
                                                              aaGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGc                                                                                 >  1:544245/1‑151 (MQ=255)
                                                               aGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCg                                                                                <  2:96975/151‑1 (MQ=255)
                                                                gACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCgg                                                                               <  1:388694/151‑1 (MQ=255)
                                                                 aCCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCgg                                                                               >  2:322216/1‑150 (MQ=255)
                                                                 aCCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                              <  1:400096/151‑1 (MQ=255)
                                                                  ccGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                              <  2:399878/150‑1 (MQ=255)
                                                                   cGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                              >  2:166506/1‑149 (MQ=255)
                                                                    gTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGt                                                                           >  1:1020006/1‑151 (MQ=255)
                                                                      ttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAGGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTg                                                                          <  1:117342/150‑1 (MQ=255)
                                                                         accaccATGCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCCGCGTCCGGAAAAGATCACCGCTGGCATCCCGAACGGGGGTGAAGCTGCGCCCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                                                                      >  2:1071261/1‑151 (MQ=255)
                                                                         accaccATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                                                                      <  1:109073/151‑1 (MQ=255)
                                                                          ccaccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg                                                                     <  1:786380/151‑1 (MQ=255)
                                                                             ccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTgg                                                                  <  2:781969/151‑1 (MQ=255)
                                                                              caTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGc                                                                 <  1:772043/151‑1 (MQ=255)
                                                                                        cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtg                                                        <  2:523820/150‑1 (MQ=255)
                                                                                        cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtg                                                        <  2:999821/150‑1 (MQ=255)
                                                                                        cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtg                                                        <  2:1161192/150‑1 (MQ=255)
                                                                                              cTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttggt                                                  >  1:14168/1‑150 (MQ=255)
                                                                                               tGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggta                                                 >  1:246228/1‑150 (MQ=255)
                                                                                                 gCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGCGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGAGACCGGTGGCGACATTCTGGTTggtaagg                                              >  1:787733/1‑151 (MQ=255)
                                                                                                    tgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggta                                            >  1:318004/1‑150 (MQ=255)
                                                                                                         gtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGcc                                       <  1:943807/150‑1 (MQ=255)
                                                                                                            ccGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaaa                                   <  1:1162227/151‑1 (MQ=255)
                                                                                                             cGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaa                                    >  2:512715/1‑149 (MQ=255)
                                                                                                                 acacCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTg                               >  2:404037/1‑150 (MQ=255)
                                                                                                                 acacCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGa                              <  2:432681/151‑1 (MQ=255)
                                                                                                                  cacCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGaa                             <  2:89830/151‑1 (MQ=255)
                                                                                                                   acCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGTTATCGTTTACATTGGAGCGGAGGTGACCGGAGGCGACATTCTCGGTGGTAAGGCAACGCCGAAAGGTGaa                             >  2:7028/1‑150 (MQ=255)
                                                                                                                     cAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACt                          >  1:719981/1‑151 (MQ=255)
                                                                                                                      aaGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTc                         <  1:274906/151‑1 (MQ=255)
                                                                                                                           gggTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCt                     <  1:938282/150‑1 (MQ=255)
                                                                                                                           gggTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTg                    <  2:1139018/151‑1 (MQ=255)
                                                                                                                                 ggAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc               >  1:282979/1‑150 (MQ=255)
                                                                                                                                    agagATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCagaa           <  1:34856/151‑1 (MQ=255)
                                                                                                                                      agaTCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCagaaga         <  1:1091415/151‑1 (MQ=255)
                                                                                                                                      agaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGCATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCaga            >  1:496761/1‑148 (MQ=255)
                                                                                                                                        aTCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAgaagaa        >  2:248636/1‑150 (MQ=255)
                                                                                                                                          cACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGaaaaa     <  1:385448/151‑1 (MQ=255)
                                                                                                                                            ccGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAAct   <  2:870349/151‑1 (MQ=255)
                                                                                                                                            ccGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAAct   <  1:740918/151‑1 (MQ=255)
                                                                                                                                             cGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActg  >  2:1029322/1‑151 (MQ=255)
                                                                                                                                             cGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActg  <  2:910726/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTG  >  NC_000913/4183656‑4183947

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