Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,085,461 A→G S151P (TCC→CCC)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,085,4610AG100.0% 265.0 / NA 72S151P (TCC→CCC) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/36);  total (36/36)

CCAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAACGC  >  NC_000913/3085326‑3085596
                                                                                                                                       |                                                                                                                                       
ccAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                      <  1:599/139‑1 (MQ=255)
ccAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                      <  2:594215/139‑1 (MQ=255)
     tAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCtt                                                                                                                                 <  1:1302459/139‑1 (MQ=255)
                        tctcGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGcc                                                                                                              >  1:81640/1‑139 (MQ=255)
                         ctcGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGccc                                                                                                             >  1:1201385/1‑139 (MQ=255)
                             ccAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaat                                                                                                         >  1:602251/1‑139 (MQ=255)
                                 ttAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatgaat                                                                                                     >  1:1067469/1‑139 (MQ=255)
                                 ttAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatgaat                                                                                                     >  2:228322/1‑139 (MQ=255)
                                          ggCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaa                                                                                                          <  1:589886/125‑1 (MQ=255)
                                          ggCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaa                                                                                                          >  2:589886/1‑125 (MQ=255)
                                          ggCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTc                                                                                            >  1:827143/1‑139 (MQ=255)
                                           gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                      >  2:1446867/1‑96 (MQ=255)
                                           gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                      <  1:1446867/96‑1 (MQ=255)
                                            cGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCaa                                                                                          >  2:1051122/1‑139 (MQ=255)
                                                     tCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcggc                                                                                 >  1:367580/1‑139 (MQ=255)
                                                       gcgGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGcggcgg                                                                               <  1:995422/139‑1 (MQ=255)
                                                         ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                             <  2:602251/139‑1 (MQ=255)
                                                         ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                             >  1:1170725/1‑139 (MQ=255)
                                                         ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCAGCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                             >  2:981203/1‑139 (MQ=255)
                                                         ggTCTTTATAACCGTTGCAGACGATGACGCCACGGGTCATGCCAGGATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGGCCCAGCGGTTCGCCCGAATTAATCAGGGACTCAATCACGCGGCGGTg                                                                             <  1:855423/139‑1 (MQ=255)
                                                          gTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGc                                                                                                               <  1:395806/104‑1 (MQ=255)
                                                          gTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGc                                                                                                               >  2:395806/1‑104 (MQ=255)
                                                             tttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCTGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTgg                                                                         <  2:1067469/139‑1 (MQ=255)
                                                              ttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  <  1:373431/129‑1 (MQ=255)
                                                              ttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  >  2:373431/1‑129 (MQ=255)
                                                              ttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                        <  2:1201385/139‑1 (MQ=255)
                                                                     cGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    <  2:484935/104‑1 (MQ=255)
                                                                     cGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    >  1:484935/1‑104 (MQ=255)
                                                                      gTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTg                                                                <  2:81640/139‑1 (MQ=255)
                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    >  2:1055873/1‑102 (MQ=255)
                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    <  1:1055873/102‑1 (MQ=255)
                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCAcg                                                                                     <  1:1126546/117‑1 (MQ=255)
                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCAcg                                                                                     >  2:1126546/1‑117 (MQ=255)
                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  >  2:238345/1‑120 (MQ=255)
                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  <  1:238345/120‑1 (MQ=255)
                                                                       ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                               <  1:228322/139‑1 (MQ=255)
                                                                         gCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATc                                                             >  1:1115854/1‑139 (MQ=255)
                                                                         gCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTAATGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATc                                                             >  1:926913/1‑139 (MQ=255)
                                                                           aGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACGTTGATCgg                                                           <  1:981203/139‑1 (MQ=255)
                                                                                 tGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                              >  2:1033933/1‑130 (MQ=255)
                                                                                 tGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                              <  1:1033933/130‑1 (MQ=255)
                                                                                 tGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAAc                                                     <  2:404021/139‑1 (MQ=255)
                                                                                 tGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAAc                                                     <  2:367580/139‑1 (MQ=255)
                                                                                  gACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAGACa                                                    >  2:1467624/1‑139 (MQ=255)
                                                                                         cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  <  2:557600/102‑1 (MQ=255)
                                                                                         cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  >  1:557600/1‑102 (MQ=255)
                                                                                            gTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAAt                                          <  2:827143/139‑1 (MQ=255)
                                                                                              cATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAGACAAGGAAGTAATCg                                        <  1:1467624/139‑1 (MQ=255)
                                                                                              cATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCg                                        >  2:498521/1‑139 (MQ=255)
                                                                                              cATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCg                                        <  2:1115854/139‑1 (MQ=255)
                                                                                                    aGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTta                                  <  1:498521/139‑1 (MQ=255)
                                                                                                      cATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTtata                                >  1:981510/1‑139 (MQ=255)
                                                                                                       aTGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGAGTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAg                               <  2:981510/139‑1 (MQ=255)
                                                                                                                tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  <  2:1376267/79‑1 (MQ=255)
                                                                                                                tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  >  1:1376267/1‑79 (MQ=255)
                                                                                                                    gCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTccc                  >  1:889722/1‑139 (MQ=255)
                                                                                                                    gCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTccc                  >  2:725244/1‑139 (MQ=255)
                                                                                                                     ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCgg                                                           <  2:175337/97‑1 (MQ=255)
                                                                                                                     ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCgg                                                           >  1:175337/1‑97 (MQ=255)
                                                                                                                     ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCt                 >  2:392437/1‑139 (MQ=255)
                                                                                                                         cAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCa             >  2:486571/1‑139 (MQ=255)
                                                                                                                          aaCTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  >  2:562960/1‑69 (MQ=255)
                                                                                                                          aaCTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                  <  1:562960/69‑1 (MQ=255)
                                                                                                                          aaCTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCAc            >  2:238776/1‑139 (MQ=255)
                                                                                                                               ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    <  1:1479466/46‑1 (MQ=255)
                                                                                                                               ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    >  2:1479466/1‑46 (MQ=255)
                                                                                                                               ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                             >  1:848586/1‑69 (MQ=255)
                                                                                                                               ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                             <  2:848586/69‑1 (MQ=255)
                                                                                                                                gCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    >  1:746015/1‑45 (MQ=255)
                                                                                                                                gCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATc                                                                                                    <  2:746015/45‑1 (MQ=255)
                                                                                                                                  tttGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAAc    <  2:1170725/139‑1 (MQ=255)
                                                                                                                                    tGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAAcgc  <  1:725244/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CCAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAACGC  >  NC_000913/3085326‑3085596

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