Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,690,087 G→T H338N (CAT→AAT)  glrK ← sensory histidine kinase GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,690,0870GT96.8% 214.8 / 0.7 62H338N (CAT→AAT) glrKsensory histidine kinase GlrK
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (28/32);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATGAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGGCTGCTGCTATCAAGAATG  >  NC_000913/2689953‑2690223
                                                                                                                                      |                                                                                                                                        
ttCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAg                                                                                                                                     >  2:1003223/1‑140 (MQ=255)
 tCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGa                                                                                                                                    <  2:269812/140‑1 (MQ=255)
 tCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGa                                                                                                                                    >  1:566317/1‑140 (MQ=255)
   aGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGaaa                                                                                                                                  >  2:993262/1‑140 (MQ=255)
         gTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTg                                                                                                                            <  2:254707/140‑1 (MQ=255)
         gTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTg                                                                                                                            <  1:993262/140‑1 (MQ=255)
               cACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCa                                                                                                                               <  2:75409/131‑1 (MQ=255)
               cACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCa                                                                                                                               >  1:75409/1‑131 (MQ=255)
               cACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTcca                                                                                                                      >  2:910781/1‑140 (MQ=255)
                  cGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTccacca                                                                                                                   <  2:727221/140‑1 (MQ=255)
                   gCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAg                                                                                                                  <  2:691688/140‑1 (MQ=255)
                    cATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGtctc                                                                                                                        <  1:333525/133‑1 (MQ=255)
                    cATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGtctc                                                                                                                        >  2:333525/1‑133 (MQ=255)
                         gAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGc                                                                                                            <  2:519919/140‑1 (MQ=255)
                         gAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGc                                                                                                            <  1:1050079/140‑1 (MQ=255)
                           gTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTa                                                                                                          <  1:1003223/140‑1 (MQ=255)
                            tAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTaa                                                                                                         <  1:628134/140‑1 (MQ=255)
                              aaGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACt                                                                                                       >  1:850661/1‑140 (MQ=255)
                                   tATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACa                                                                                                  >  2:545070/1‑140 (MQ=255)
                                    aTCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACAt                                                                                                 >  2:262941/1‑140 (MQ=255)
                                    aTCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACAt                                                                                                 >  1:380298/1‑140 (MQ=255)
                                      ccAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATtc                                                                                               >  2:748917/1‑140 (MQ=255)
                                               tcatcaGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGATACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTcc                                                                                      >  2:90740/1‑140 (MQ=255)
                                               tcatcaGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTcc                                                                                      >  2:106541/1‑140 (MQ=255)
                                                    agcagcATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCg                                                                                  <  1:884522/139‑1 (MQ=255)
                                                    agcagcATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCg                                                                                  >  2:884522/1‑139 (MQ=255)
                                                          aTTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTc                                                                           <  1:106541/140‑1 (MQ=255)
                                                           ttGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTcc                                                                          >  1:68488/1‑140 (MQ=255)
                                                           ttGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTcc                                                                          >  1:664063/1‑140 (MQ=255)
                                                            tGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCg                                                                         <  2:380298/140‑1 (MQ=255)
                                                            tGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCg                                                                         <  1:262941/140‑1 (MQ=255)
                                                             ggCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGc                                                                        <  1:247189/140‑1 (MQ=255)
                                                                tCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATGAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTg                                                                     >  1:1018012/1‑140 (MQ=255)
                                                                  cGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGtt                                                                   <  1:704837/140‑1 (MQ=255)
                                                                  cGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGtt                                                                   <  1:80181/140‑1 (MQ=255)
                                                                   gCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGttt                                                                  >  2:911125/1‑140 (MQ=255)
                                                                    ccAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTa                                                                 <  2:163788/140‑1 (MQ=255)
                                                                    ccAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTa                                                                 <  2:206553/140‑1 (MQ=255)
                                                                       ggCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACgg                                                              <  1:401153/140‑1 (MQ=255)
                                                                        gCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGt                                                             >  2:1016549/1‑140 (MQ=255)
                                                                        gCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGt                                                             >  2:86885/1‑140 (MQ=255)
                                                                                ttGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCa                                                     >  2:746665/1‑140 (MQ=255)
                                                                                  gCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAg                                                   <  2:943385/140‑1 (MQ=255)
                                                                                      tGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGATACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAgt                                               <  1:90740/140‑1 (MQ=255)
                                                                                          aTCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAgttgtt                                           >  2:822238/1‑140 (MQ=255)
                                                                                             gACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGa                                        >  1:629225/1‑140 (MQ=255)
                                                                                                      aTGCATCATTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGttttt                               <  1:86885/140‑1 (MQ=255)
                                                                                                          atcatTTTAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCaa                           >  1:758812/1‑140 (MQ=255)
                                                                                                                 tAGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCgg                    <  1:822238/140‑1 (MQ=255)
                                                                                                                  aGCCCGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgc                   <  2:701502/140‑1 (MQ=255)
                                                                                                                    cccGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctg                 <  2:850661/140‑1 (MQ=255)
                                                                                                                    cccGTGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctg                 <  2:519461/140‑1 (MQ=255)
                                                                                                                      cGTGCGGGCAGGCTATTAGCAGAAACCCCTGTCCCCACCAGCGGTGCTAACTTAACATTCTCCAGTTCCCCCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgct               <  2:566317/140‑1 (MQ=255)
                                                                                                                        tGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgctgc             >  1:698708/1‑140 (MQ=255)
                                                                                                                        tGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgctgc             >  1:676329/1‑140 (MQ=255)
                                                                                                                        tGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgctgc             >  1:639251/1‑140 (MQ=255)
                                                                                                                        tGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgctgc             >  1:532734/1‑140 (MQ=255)
                                                                                                                        tGCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgctgc             >  2:134115/1‑140 (MQ=255)
                                                                                                                         gCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgctgct            <  2:532734/140‑1 (MQ=255)
                                                                                                                         gCGGGCAGGCTATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGgctgctgct            <  1:859784/140‑1 (MQ=255)
                                                                                                                                ggCTATGAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGGCTGCTGCTATCAAGa     <  2:1018012/140‑1 (MQ=255)
                                                                                                                                   tATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGGCTGCTGCTATCAAGAATg  <  2:758812/140‑1 (MQ=255)
                                                                                                                                   tATTAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGGCTGCTGCTATCAAGAATg  <  1:748917/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
TTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGAGATCGACGTCGGTATGCATCATTTTAGCCCGTGCGGGCAGGCTATGAGCAGAAACCACTGTCTCCACCAGCGGTGCTAACTCAACATTCTCCAGTTCCACCGCACTGTCCGCCTGTTTACGGTTGTAATCAAGCAGTTGTTCGATCAGTTTTTGCAAATTGCGGCTGCTGCTATCAAGAATG  >  NC_000913/2689953‑2690223

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