Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,250,104 Δ1 bp coding (551/921 nt) malM → maltose regulon periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,250,1020G.98.4% 266.7 / ‑2.0 63coding (549/921 nt)malMmaltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (31/31);  total (32/31)
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.59e-01

GCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCC  >  NC_000913/4249968‑4250236
                                                                                                                                      |                                                                                                                                      
gCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAggcgtc                                                                                                                                   >  2:886574/1‑136 (MQ=255)
  aGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCgg                                                                                                                                >  2:834652/1‑140 (MQ=255)
   gATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGt                                                                                                                               >  2:457933/1‑140 (MQ=255)
   gATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGt                                                                                                                               <  2:471195/140‑1 (MQ=255)
       ggCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACt                                                                                                                           <  2:89469/140‑1 (MQ=255)
       ggCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACt                                                                                                                           <  1:1110623/140‑1 (MQ=255)
       ggCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACt                                                                                                                           >  2:122951/1‑140 (MQ=255)
       ggCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACt                                                                                                                           <  2:684874/140‑1 (MQ=255)
         cTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGGAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCg                                                                                                                         >  2:1129333/1‑140 (MQ=255)
         cTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTGTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCg                                                                                                                         >  2:1268822/1‑140 (MQ=255)
            gAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATc                                                                                                                      <  2:1398180/140‑1 (MQ=255)
                 cGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGGAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGa                                                                                                                 <  1:1129333/140‑1 (MQ=255)
                 cGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGa                                                                                                                 <  1:122951/140‑1 (MQ=255)
                   ttATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCg                                                                                                                                 <  2:1408656/122‑1 (MQ=255)
                   ttATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCg                                                                                                                                 >  1:1408656/1‑122 (MQ=255)
                       gcgcCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATccc                                                                                                           <  2:628259/140‑1 (MQ=255)
                       gcgcCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATccc                                                                                                           <  2:304952/140‑1 (MQ=255)
                          ccTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGa                                                                                                        >  2:550571/1‑140 (MQ=255)
                              acacCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCg                                                                                                    >  1:558289/1‑140 (MQ=255)
                                  cGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTg                                                                                                >  2:1208222/1‑140 (MQ=255)
                                  cGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTg                                                                                                <  2:255174/140‑1 (MQ=255)
                                    gCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCt                                                                                              <  1:191609/140‑1 (MQ=255)
                                    gCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCt                                                                                              <  2:343427/140‑1 (MQ=255)
                                     cGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTc                                                                                             <  1:886574/140‑1 (MQ=255)
                                     cGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGAATATCCCCGATCCGGTTGCTc                                                                                             >  2:548256/1‑140 (MQ=255)
                                          gggCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCAt                                                                                        <  1:39159/140‑1 (MQ=255)
                                             cagcaAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATacc                                                                                     >  1:297583/1‑140 (MQ=255)
                                               gcaAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATaccac                                                                                   >  1:404687/1‑140 (MQ=255)
                                                       tttATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCt                                                                           <  2:297583/140‑1 (MQ=255)
                                                             ttCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactga                                                                     <  2:715987/140‑1 (MQ=255)
                                                                 ggTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactgaaact                                                                 >  2:1131314/1‑140 (MQ=255)
                                                                 ggTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactgaaact                                                                 <  1:1131314/140‑1 (MQ=255)
                                                                 ggTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactgaaact                                                                 >  1:206363/1‑140 (MQ=255)
                                                                     tttACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactgaaactgaa                                                              >  2:1103831/1‑140 (MQ=255)
                                                                            cGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGaaaa                                                      >  2:588004/1‑140 (MQ=255)
                                                                            cGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGaaaa                                                      >  1:1087634/1‑140 (MQ=255)
                                                                             ggAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAAc                                                     <  1:834652/140‑1 (MQ=255)
                                                                                   aGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActc                                               <  2:206363/140‑1 (MQ=255)
                                                                                      tctcCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActccag                                            <  1:550571/140‑1 (MQ=255)
                                                                                        tcCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagct                                          <  1:1268822/140‑1 (MQ=255)
                                                                                        tcCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagct                                          >  2:419749/1‑140 (MQ=255)
                                                                                         ccAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagctc                                         >  2:589355/1‑140 (MQ=255)
                                                                                           agcagACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagctcca                                       <  2:404687/140‑1 (MQ=255)
                                                                                            gcagACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATATCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagctccag                                      <  2:1373240/140‑1 (MQ=255)
                                                                                               gacgacCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCgt                                   >  1:1088372/1‑140 (MQ=255)
                                                                                                  gacCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGtt                                >  1:274256/1‑140 (MQ=255)
                                                                                                    cccAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTgg                              >  1:330112/1‑140 (MQ=255)
                                                                                                     ccAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGAATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGt                             <  1:548256/140‑1 (MQ=255)
                                                                                                       aaCTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAg                           >  1:1174687/1‑140 (MQ=255)
                                                                                                         cTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGa                         >  2:1273436/1‑140 (MQ=255)
                                                                                                         cTGCTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGa                         >  2:94055/1‑140 (MQ=255)
                                                                                                            cTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCt                      >  2:423654/1‑140 (MQ=255)
                                                                                                            cTCGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCt                      >  2:810018/1‑140 (MQ=255)
                                                                                                              cGACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCttt                    >  2:1045041/1‑140 (MQ=255)
                                                                                                               gACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTa                   <  1:423654/140‑1 (MQ=255)
                                                                                                               gACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTa                   <  1:589355/140‑1 (MQ=255)
                                                                                                               gACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTa                   <  1:1045041/140‑1 (MQ=255)
                                                                                                               gACCCGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTACAGCGTGTTGGTAGGACCTTTa                   <  1:457933/140‑1 (MQ=255)
                                                                                                                 cccGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTAtt                 >  1:14924/1‑140 (MQ=255)
                                                                                                                 cccGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTAtt                 >  1:129544/1‑140 (MQ=255)
                                                                                                                  ccGGCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTAttt                <  1:810018/140‑1 (MQ=255)
                                                                                                                    ggCTAAAGCCTATGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTgg              <  2:558289/140‑1 (MQ=255)
                                                                                                                                tGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAACCGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTcc  >  2:283099/1‑140 (MQ=255)
                                                                                                                                tGCCAAGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTcc  >  2:301720/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCC  >  NC_000913/4249968‑4250236

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