Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,183,669 G→T 94.0% G809C (GGT→TGT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,6690GT94.0% 294.4 / 5.0 98G809C (GGT→TGT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (48/43);  total (52/46)
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.78e-01

AACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTG  >  NC_000913/4183522‑4183815
                                                                                                                                                   |                                                                                                                                                  
aaCCAGAACACCTGTATCAACCAGATGCCGGGGGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGt                                                                                                                                                  <  2:1057841/150‑1 (MQ=255)
 aCCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGtt                                                                                                                                                 <  1:704587/150‑1 (MQ=255)
 aCCAGAACACCTGTATCAACCAGATGCCGTGTGTGGCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGtt                                                                                                                                                 <  2:1469184/150‑1 (MQ=255)
 aCCAGAACACCTGTATCAACCAGATGCCGTGTGTGGCTCTGGGGGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGtt                                                                                                                                                 <  2:956398/150‑1 (MQ=255)
  ccAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTa                                                                                                                                                >  2:1343586/1‑150 (MQ=255)
  ccAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTa                                                                                                                                                >  2:38674/1‑150 (MQ=255)
  ccAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTa                                                                                                                                                >  2:1490048/1‑150 (MQ=255)
   cAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTAc                                                                                                                                               <  1:47490/150‑1 (MQ=255)
         acCTGTATCAACCAGAGGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTGCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTc                                                                                                                                         <  2:1057894/150‑1 (MQ=255)
          cctcaaacAACCAGAAACCGTCTGTCTCTCTGGCTGAACCGGTTGAACCTGGCGACGTCCTCGCAGACGGTCCGTCCACCGACCACGGTGAACTGGCGCATCCACACAACATGCGCGTAGCGTTCAAGCCGTGCAATAGTTACAACATCg                                                                                                                                        >  2:964576/8‑150 (MQ=255)
          ccTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGCCCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAAGTGTTACAACTTCg                                                                                                                                        >  1:1388079/1‑150 (MQ=255)
          ccTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCg                                                                                                                                        >  1:1382380/1‑150 (MQ=255)
               aTCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGAc                                                                                                                                   >  2:1633165/1‑150 (MQ=255)
                 cAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACTGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTc                                                                                                                                 >  2:772380/1‑150 (MQ=255)
                 cAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTc                                                                                                                                 >  2:1624270/1‑150 (MQ=255)
                  aaCCAGATGTCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTcc                                                                                                                                >  2:1387681/1‑150 (MQ=255)
                  aaCCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTcc                                                                                                                                >  1:555571/1‑150 (MQ=255)
                         tGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCTTCCACCGACCTCGGTGAACTGGCGCTTGGGCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCg                                                                                                                         <  2:1042214/150‑1 (MQ=255)
                            cGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGCTATTGTTACAACTTCGAAGACTCCATCCTCGTAt                                                                                                                      >  1:848460/1‑150 (MQ=255)
                             gtgtgtgtCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGTCGCTTGGTCAGACCATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTACCCAGACTCCATCCGCGTATc                                                                                                                     >  2:1152866/1‑150 (MQ=255)
                                 gtgtCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAg                                                                                                                 >  1:601446/1‑150 (MQ=255)
                                   gtCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCg                                                                                                               <  1:1577575/150‑1 (MQ=255)
                                     ctctGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCgtg                                                                                                             >  2:644245/1‑150 (MQ=255)
                                      tctGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTAGTATCCGAGCgtgt                                                                                                            >  2:1273242/1‑150 (MQ=255)
                                      tctGGGGTAACCGGTTGAACGTGGCGACGGGCTGGCTGACGGTCCGGCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTGCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCgtgt                                                                                                            <  1:310505/150‑1 (MQ=255)
                                         gggTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACATCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGtgttgt                                                                                                         <  1:745182/150‑1 (MQ=255)
                                          ggTGAACCGGTTGCACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGCACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTgttgtt                                                                                                        >  2:238081/1‑150 (MQ=255)
                                          ggTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTgttgtt                                                                                                        <  2:782329/150‑1 (MQ=255)
                                          ggTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTgttgtt                                                                                                        >  1:925719/1‑150 (MQ=255)
                                           gTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTc                                                                                                       >  1:1410109/1‑150 (MQ=255)
                                              aaCCGGTTGAACGTGGCGACGTGCGGGCAGACGGGACGTCCACCGACCTCGGTGAAATGGCGCGTGGGCAGAACATGCGCGCAGCGTTCATGCGGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAgg                                                                                                    <  1:919209/150‑1 (MQ=255)
                                                  ggTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGa                                                                                                >  1:384803/1‑150 (MQ=255)
                                                  ggTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGa                                                                                                >  2:760660/1‑150 (MQ=255)
                                                  ggTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGa                                                                                                >  1:391138/1‑150 (MQ=255)
                                                      gAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGt                                                                                            >  2:1733175/1‑150 (MQ=255)
                                                      gAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGt                                                                                            >  2:1709165/1‑150 (MQ=255)
                                                       aaCGTGGCGACGTGCTGGCAGACGGTCCGTCCACAGACCTCGGTGAACTGGCGCTTTGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGtt                                                                                           <  1:324198/150‑1 (MQ=255)
                                                         cGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTGCAGGAAGACCGTTTc                                                                                         >  1:1380259/1‑150 (MQ=255)
                                                         cGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGGCAGAACCTGCGCGGAGCGGTCATGCCGGGGAATTGTTACAACTTCGAAGACTCCATCCTCGAATCCGAGCGTGTTGTTCAGGAAGACCGTTTc                                                                                         >  1:276565/1‑150 (MQ=255)
                                                          ggggCGCCGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATCCGCGTAGCGTTCATGCCGTGTAATTGTTACAACTTCGACGACTCCATCCTCGGATCCGAGCGTGTGGTTCAGGAAGACCGTTTca                                                                                        <  2:97700/148‑1 (MQ=255)
                                                          gTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTca                                                                                        <  1:1084470/150‑1 (MQ=255)
                                                           tGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcac                                                                                       >  2:1063279/1‑150 (MQ=255)
                                                            ggCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcacc                                                                                      <  1:849619/150‑1 (MQ=255)
                                                            ggCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcacc                                                                                      <  2:1209614/150‑1 (MQ=255)
                                                            ggCGACGTGCTGGCAGACGGTCCGGCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcacc                                                                                      <  1:843947/150‑1 (MQ=255)
                                                              cGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcaccac                                                                                    <  2:1580596/150‑1 (MQ=255)
                                                               gACGTGCTGGCAGACGGTCCGTCCACCGACCGCGGGGTACTGGCGCTTGGTCAGAACATGCGCGAAGCGGTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGGGCTGTTCAGGAAGACCGTTTCaccacc                                                                                   <  1:774013/150‑1 (MQ=255)
                                                                  gTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATc                                                                                <  2:618797/150‑1 (MQ=255)
                                                                   tGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATcc                                                                               <  1:404498/150‑1 (MQ=255)
                                                                   tGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATcc                                                                               <  2:1154787/150‑1 (MQ=255)
                                                                     cTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCcac                                                                             >  1:622212/1‑150 (MQ=255)
                                                                     cTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCcac                                                                             <  2:973369/150‑1 (MQ=255)
                                                                        gCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACAtt                                                                          >  2:543914/1‑150 (MQ=255)
                                                                         cAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTc                                                                         >  2:1462792/1‑150 (MQ=255)
                                                                           gACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAg                                                                       >  1:1176934/1‑150 (MQ=255)
                                                                             cGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGa                                                                     <  1:1213274/150‑1 (MQ=255)
                                                                                  cgtccACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCCCCCCCATCCACAGTCAGGAACTgg                                                                >  2:834603/1‑150 (MQ=255)
                                                                                   gtccACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGc                                                               >  1:29103/1‑150 (MQ=255)
                                                                                   gtccACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGc                                                               <  1:1624270/150‑1 (MQ=255)
                                                                                      cACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCgtg                                                            <  2:1382380/150‑1 (MQ=255)
                                                                                       aCCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCgtgt                                                           >  2:983879/1‑150 (MQ=255)
                                                                                        ccGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGCCCGTTTCACCACCATCCACATTCAGGAACTGGCgtgtg                                                          >  2:643042/1‑150 (MQ=255)
                                                                                        ccGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCgtgtg                                                          <  1:1343586/150‑1 (MQ=255)
                                                                                        ccGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCgtgtg                                                          <  1:1490048/150‑1 (MQ=255)
                                                                                            ccTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGTACTGGCGTGTGTGTc                                                      >  2:39419/1‑150 (MQ=255)
                                                                                             cTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTcc                                                     <  2:1099606/150‑1 (MQ=255)
                                                                                                ggTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGt                                                  >  1:1020558/1‑150 (MQ=255)
                                                                                                      cTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACAcc                                            <  1:772380/150‑1 (MQ=255)
                                                                                                        ggCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCaa                                          >  1:1652096/1‑150 (MQ=255)
                                                                                                         gcggTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAg                                         <  2:1693942/150‑1 (MQ=255)
                                                                                                          cgcTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACTACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTAAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTACCGTGCCACCaatc                                        >  2:1004783/1‑148 (MQ=255)
                                                                                                            cTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTg                                      <  1:409327/150‑1 (MQ=255)
                                                                                                             ttGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTgg                                     <  2:1403675/150‑1 (MQ=255)
                                                                                                             ttGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTggg                                    <  2:656102/150‑1 (MQ=255)
                                                                                                              tGGTCATAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGGGTGTGTCCCGTGACACCAAGCTggg                                    <  2:311185/150‑1 (MQ=255)
                                                                                                              tGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTggg                                    <  1:644245/150‑1 (MQ=255)
                                                                                                                gTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCg                                                              <  2:398181/122‑1 (MQ=255)
                                                                                                                gTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCg                                                              >  1:398181/1‑122 (MQ=255)
                                                                                                                 tCAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTcc                                 <  1:863994/150‑1 (MQ=255)
                                                                                                                  cAGAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCg                                >  2:711346/1‑150 (MQ=255)
                                                                                                                    gAACATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCCCCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGa                              >  2:829826/1‑150 (MQ=255)
                                                                                                                      aCATGCGCGTAGCGTGCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCGCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAag                            <  1:73367/150‑1 (MQ=255)
                                                                                                                       cATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCTATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAaga                           >  1:506457/1‑150 (MQ=255)
                                                                                                                       cATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAaga                           >  1:778555/1‑150 (MQ=255)
                                                                                                                       cATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAaga                           >  1:1116580/1‑150 (MQ=255)
                                                                                                                       cATGCGCGTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACAACAAGCTGGGTCCGGAaga                           >  1:108014/1‑150 (MQ=255)
                                                                                                                         ttcgcgTAGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCCGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGGGTGTCCCGTGACACCAAGCTGGGTCCGGAAgaga                         <  2:225538/148‑1 (MQ=255)
                                                                                                                           cgcgTAGCGTTCAAGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGATCAGCAAGACCGTTTCACAACCATCCACATTCAGGAACTGGCGTGTCTGTCCCGTGACACCAACCTGGGTCCGGAAGAGAta                       >  2:976091/1‑149 (MQ=255)
                                                                                                                                aGCGTTCATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCCCCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGc                  >  1:963957/1‑150 (MQ=255)
                                                                                                                                 gCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCt                 <  2:678079/150‑1 (MQ=255)
                                                                                                                                  cGTTCATGCCCTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCCCCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGCTCCCCGcgg                >  1:405747/1‑148 (MQ=255)
                                                                                                                                      cATGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACAt            <  2:925719/150‑1 (MQ=255)
                                                                                                                                        tGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATcc          <  1:1562666/150‑1 (MQ=255)
                                                                                                                                        tGCCGTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATcc          <  2:1176934/150‑1 (MQ=255)
                                                                                                                                         gCCGTGGAATTGTTACAACTTCGAATACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGAGTGTGTGTCCCGTGACACCAAGCTGGGTCAGGAAGAGAGCACCGCTGACATacc         >  1:81793/1‑150 (MQ=255)
                                                                                                                                         gCCGTGGAATTGTAACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGATGACAGTTTCACCACGATCCACATTCAGGAGCTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGTAAGAGATCAACGCTGACATccc         >  1:263504/1‑150 (MQ=255)
                                                                                                                                            gTGGAATTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGaa      >  2:1455422/1‑150 (MQ=255)
                                                                                                                                                aaTTGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTg  <  2:967735/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTG  >  NC_000913/4183522‑4183815

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