Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,440,212 T→A E273V (GAA→GTA)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,2120TA100.0% 248.2 / NA 65E273V (GAA→GTA) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/35);  total (30/35)

CAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTC  >  NC_000913/3440068‑3440354
                                                                                                                                                |                                                                                                                                              
cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGCAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCtt                                                                                                                                          >  1:692902/1‑151 (MQ=255)
         aGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGt                                                                                                                                 <  1:568091/151‑1 (MQ=255)
          ggCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGt                                                                                                                                 <  1:572104/150‑1 (MQ=255)
           gcgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCATTTa                                                                                                                               <  2:851600/151‑1 (MQ=255)
           gcgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTa                                                                                                                               <  2:434403/151‑1 (MQ=255)
           gcgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTa                                                                                                                               <  1:439039/151‑1 (MQ=255)
            cgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAg                                                                                                                              <  2:784603/151‑1 (MQ=255)
            cgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAg                                                                                                                              <  2:757132/151‑1 (MQ=255)
              cATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGc                                                                                                                             >  1:859217/1‑150 (MQ=255)
              cATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCa                                                                                                                            <  2:980761/151‑1 (MQ=255)
                 gCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCaga                                                                                                                          <  1:438334/150‑1 (MQ=255)
                 gCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCag                                                                                                                           >  1:655756/1‑149 (MQ=255)
                            gTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCaa                                                                                                              <  2:480349/151‑1 (MQ=255)
                                     aaGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTcc                                                                                                          >  1:175636/1‑146 (MQ=255)
                                       gCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGt                                                                                                   <  1:831474/151‑1 (MQ=255)
                                         cAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGt                                                                                                   <  2:488141/149‑1 (MQ=255)
                                          aGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCaa                                                                                                <  1:840285/151‑1 (MQ=255)
                                            cACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACa                                                                                              <  2:285219/151‑1 (MQ=255)
                                               gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAggg                                                                                           <  1:135846/151‑1 (MQ=255)
                                                 cTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                                                          >  2:571633/1‑150 (MQ=255)
                                                   ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgcgc                                                                                        <  2:482948/150‑1 (MQ=255)
                                                   ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgcgc                                                                                        >  1:941185/1‑150 (MQ=255)
                                                   ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgcgc                                                                                        >  1:44145/1‑150 (MQ=255)
                                                     aaTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCgcag                                                                                      <  1:39159/150‑1 (MQ=255)
                                                        ctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAgg                                                                                  >  2:10504/1‑151 (MQ=255)
                                                          cAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggat                                                                                >  2:745012/1‑151 (MQ=255)
                                                          cAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAgga                                                                                 <  1:631998/150‑1 (MQ=255)
                                                           aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                                                               <  2:118721/151‑1 (MQ=255)
                                                           aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                                                               <  1:865159/151‑1 (MQ=255)
                                                             tAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgg                                                                             >  2:599154/1‑151 (MQ=255)
                                                              aaGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgg                                                                             <  2:125757/150‑1 (MQ=255)
                                                               agagATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggat                                                                           <  2:665679/151‑1 (MQ=255)
                                                                gagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                                                          >  1:857605/1‑151 (MQ=255)
                                                                 agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGgatcggatcg                                                                         >  1:115868/1‑151 (MQ=255)
                                                                    ttttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAAc                                                                      >  1:347450/1‑151 (MQ=255)
                                                                      ttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAAct                                                                     <  1:761381/150‑1 (MQ=255)
                                                                           ccAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                                                               >  2:305194/1‑151 (MQ=255)
                                                                             aaGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtc                                                             >  1:126347/1‑151 (MQ=255)
                                                                             aaGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtc                                                             <  1:479223/151‑1 (MQ=255)
                                                                               ggatAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTAATGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGCCAACCGGGCGCAGCAGGAACGGAGCGAACTCTGGTTtc                                                             >  2:38234/4‑149 (MQ=255)
                                                                               ggTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtctc                                                           >  2:140719/1‑151 (MQ=255)
                                                                                gTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTctct                                                          <  1:820899/151‑1 (MQ=255)
                                                                                  tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtctt                                                         <  1:332231/150‑1 (MQ=255)
                                                                                   aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                                                       >  2:859715/1‑151 (MQ=255)
                                                                                       gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTc                                                    >  2:331504/1‑150 (MQ=255)
                                                                                          ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACtt                                                >  2:780876/1‑151 (MQ=255)
                                                                                             aGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTc                                               >  1:176050/1‑149 (MQ=255)
                                                                                               gAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGc                                           <  2:675580/151‑1 (MQ=255)
                                                                                                  cTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg                                         >  1:850250/1‑150 (MQ=255)
                                                                                                    cAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGa                                        <  2:467748/149‑1 (MQ=255)
                                                                                                      aCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGt                                     >  1:166164/1‑150 (MQ=255)
                                                                                                         tCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACAt                                 >  2:411941/1‑151 (MQ=255)
                                                                                                             tACGCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACg                             >  1:767157/1‑151 (MQ=255)
                                                                                                                gCTGTACCAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTa                           >  2:885786/1‑150 (MQ=255)
                                                                                                                      ccAGATCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCaa                     <  2:887193/150‑1 (MQ=255)
                                                                                                                           tCACCGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaa                 >  2:279063/1‑149 (MQ=255)
                                                                                                                              ccGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCtt            <  2:838730/151‑1 (MQ=255)
                                                                                                                               cGATATAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTc           >  1:240619/1‑151 (MQ=255)
                                                                                                                                 atatAGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTcc          >  2:672515/1‑150 (MQ=255)
                                                                                                                                  tataGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTcc          >  2:697106/1‑149 (MQ=255)
                                                                                                                                   ataGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgt       <  2:447650/151‑1 (MQ=255)
                                                                                                                                   ataGTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAg        <  2:475237/150‑1 (MQ=255)
                                                                                                                                      gTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgttgt    <  2:151685/151‑1 (MQ=255)
                                                                                                                                      gTGGATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgttgt    <  2:220145/151‑1 (MQ=255)
                                                                                                                                        ggATAGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTc  <  1:426133/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTC  >  NC_000913/3440068‑3440354

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: