Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,922,483 G→A intergenic (‑43/+574) atpI ← / ← rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,922,4830GA100.0% 237.1 / NA 60intergenic (‑43/+574)atpI/rsmGATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/26);  total (34/26)

AGAAGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGC  >  NC_000913/3922335‑3922629
                                                                                                                                                    |                                                                                                                                                  
aGAAGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCAta                                                                                                                                                  <  1:33872/151‑1 (MQ=255)
   aGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCAtata                                                                                                                                                >  1:336580/1‑150 (MQ=255)
    ggggTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATaaa                                                                                                                                              <  1:613755/151‑1 (MQ=255)
             gAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActttctt                                                                                                                                      <  1:817576/150‑1 (MQ=255)
                 cTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAgg                                                                                                                                 >  2:30722/1‑151 (MQ=255)
                       aGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCtt                                                                                                                              >  1:562782/1‑148 (MQ=255)
                          aaTCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTc                                                                                                                        >  1:912016/1‑151 (MQ=255)
                          aaTCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTc                                                                                                                        <  1:796196/151‑1 (MQ=255)
                          aaTCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTc                                                                                                                        <  1:278359/151‑1 (MQ=255)
                            tCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCa                                                                                                                       >  2:847150/1‑150 (MQ=255)
                              cACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGc                                                                                                                    <  1:95322/151‑1 (MQ=255)
                                      aTCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCaaaa                                                                                                             >  2:488871/1‑150 (MQ=255)
                                       tCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCAACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCaaaaa                                                                                                            <  2:938075/150‑1 (MQ=255)
                                        caccacCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGAGATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAg                                                                                                           <  1:736181/150‑1 (MQ=255)
                                         accaccAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCg                                                                                                         <  1:795429/151‑1 (MQ=255)
                                           caccaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGt                                                                                                       >  2:446661/1‑151 (MQ=255)
                                           caccaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGt                                                                                                       >  2:345148/1‑151 (MQ=255)
                                           caccaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGt                                                                                                       <  1:61837/151‑1 (MQ=255)
                                             ccaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCa                                                                                                     <  1:79491/151‑1 (MQ=255)
                                                 tAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAtt                                                                                                 >  1:509295/1‑151 (MQ=255)
                                                      gAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACg                                                                                            >  1:769056/1‑151 (MQ=255)
                                                       aaCGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACgg                                                                                           <  1:146534/151‑1 (MQ=255)
                                                          gAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTgc                                                                                        >  2:256474/1‑151 (MQ=255)
                                                           aGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTgcg                                                                                       >  1:162198/1‑151 (MQ=255)
                                                                  gCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGt                                                                                <  2:825224/151‑1 (MQ=255)
                                                                   cTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTg                                                                               >  1:211828/1‑151 (MQ=255)
                                                                      ccGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGa                                                                              <  2:894092/149‑1 (MQ=255)
                                                                      ccGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGAtt                                                                            >  2:700852/1‑151 (MQ=255)
                                                                      ccGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGAtt                                                                            >  2:483732/1‑151 (MQ=255)
                                                                        gAGCAATGTCTCGACTCACGAGCGCCACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGGGTCATATAAATCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTc                                                                          >  2:626506/1‑151 (MQ=255)
                                                                        gAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTc                                                                          >  1:521935/1‑151 (MQ=255)
                                                                          gCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCa                                                                         >  2:433458/1‑150 (MQ=255)
                                                                           cAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCaa                                                                        >  1:950721/1‑150 (MQ=255)
                                                                           cAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCaa                                                                        >  2:375420/1‑150 (MQ=255)
                                                                                tttCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAAt                                                                   <  2:371416/150‑1 (MQ=255)
                                                                                    gACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAg                                                                >  2:149775/1‑149 (MQ=255)
                                                                                     aCTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGt                                                               >  1:843808/1‑149 (MQ=255)
                                                                                        cACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCa                                                          <  1:944468/151‑1 (MQ=255)
                                                                                         aCGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCa                                                          >  2:11319/1‑150 (MQ=255)
                                                                                           gAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaa                                                        <  2:608585/150‑1 (MQ=255)
                                                                                                acacAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGa                                                  <  1:553235/151‑1 (MQ=255)
                                                                                                 cacaGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGa                                                  <  2:222026/150‑1 (MQ=255)
                                                                                                  acaGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAAt                                                <  2:742439/151‑1 (MQ=255)
                                                                                                        aTCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGt                                          <  1:949331/151‑1 (MQ=255)
                                                                                                         tCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGtt                                         >  1:835344/1‑151 (MQ=255)
                                                                                                          cACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGtt                                         <  1:68057/150‑1 (MQ=255)
                                                                                                                       ccTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCa                           <  2:615783/151‑1 (MQ=255)
                                                                                                                          gCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTg                        >  2:41673/1‑151 (MQ=255)
                                                                                                                            tCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGc                       >  1:250886/1‑150 (MQ=255)
                                                                                                                             cccTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCAt                     >  2:882713/1‑151 (MQ=255)
                                                                                                                             cccTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCAt                     <  1:380358/151‑1 (MQ=255)
                                                                                                                                ttCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATAc                   >  2:139878/1‑150 (MQ=255)
                                                                                                                                   gAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAtt               >  1:755141/1‑151 (MQ=255)
                                                                                                                                       tATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCg            >  2:540822/1‑150 (MQ=255)
                                                                                                                                            cgcgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACtttt      >  1:785742/1‑151 (MQ=255)
                                                                                                                                            cgcgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACtttt      <  2:791928/151‑1 (MQ=255)
                                                                                                                                              cgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCt    <  2:813302/151‑1 (MQ=255)
                                                                                                                                               gtgtCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTg   >  1:95962/1‑151 (MQ=255)
                                                                                                                                               gtgtCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTg   >  2:470820/1‑151 (MQ=255)
                                                                                                                                                tgtCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCACTGCATACAATTGCGACTTTTCTGc  >  2:704808/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGAAGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGC  >  NC_000913/3922335‑3922629

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