Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,592:1 +C 100% coding (960/993 nt) rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,5881.C98.0% 368.7 / 0.1 102coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (52/48);  total (53/49)
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.86e-01

Only 100 of 107 total aligned reads displayed.

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑CCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682
                                                                                                      |                                                                                               
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑g                                                                                                  <  2:1711956/100‑1 (MQ=255)
   cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑ccc                                                                                              >  2:1027416/1‑100 (MQ=255)
       ccTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑ccc‑cctg                                                                                         >  1:947002/1‑97 (MQ=255)
       ccTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑ccc‑cctg                                                                                         <  1:1783009/100‑4 (MQ=255)
         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑ccc‑cctgcg                                                                                       >  1:1822470/1‑95 (MQ=255)
             aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTg                                                                                    >  1:419345/1‑100 (MQ=255)
             aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTg                                                                                    <  1:2251076/100‑1 (MQ=255)
              ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGc                                                                                   <  1:1135172/100‑1 (MQ=255)
               tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCa                                                                                  <  2:1143092/100‑1 (MQ=255)
               tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCa                                                                                  >  2:1444328/1‑100 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAg                                                                                 <  2:946997/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAg                                                                                 <  1:1529951/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAg                                                                                 <  1:2044993/100‑1 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAgg                                                                                >  2:1539489/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAgg                                                                                >  1:300091/1‑100 (MQ=255)
                   nacaGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAt                                                                              >  1:2070564/2‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAt                                                                              >  1:889537/1‑100 (MQ=255)
                    acaGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAtt                                                                             >  1:510358/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAttt                                                                            >  2:48961/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAttt                                                                            >  2:1625018/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAttt                                                                            <  1:1861167/100‑1 (MQ=255)
                     caGAAAACACCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCCTACTTACTCGCGGAACAGCGCTTCGATCTTCACC‑CCC‑CTGCGTCTCAAGGAtt                                                                             >  2:1401325/1‑99 (MQ=255)
                      aGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTc                                                                           >  2:1677781/1‑100 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTcg                                                                          <  1:905573/100‑1 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTcg                                                                          >  2:556336/1‑100 (MQ=255)
                        aaaaGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTcgc                                                                         <  2:1422561/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTcgcg                                                                        <  2:1117526/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTcgcg                                                                        <  2:907793/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTcgcg                                                                        <  2:845989/100‑1 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCgcgc                                                                       <  1:1250886/100‑1 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCaa                                                                     >  2:1858668/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:74603/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:674382/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   <  1:1318661/100‑1 (MQ=255)
                               cAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTCCTCGCGGAACAGCGCTTCGACATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACg                                                                  >  1:22207/1‑100 (MQ=255)
                                   cTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  1:917071/100‑1 (MQ=255)
                                     cGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  2:1723460/100‑1 (MQ=255)
                                         cgaggcTGGCCTTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:484173/95‑1 (MQ=255)
                                          gagaCTGGCC‑TTTCGTACAGATGCTTACTTCCTAGTTGAACTGCGCTTCGATGTCCA‑G‑CCCTCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  2:1515870/100‑1 (MQ=255)
                                           agaCTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  1:441204/1‑100 (MQ=255)
                                           agaCTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  1:347969/1‑100 (MQ=255)
                                           agaCTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      <  2:286187/100‑1 (MQ=255)
                                            gaCTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:1860309/1‑100 (MQ=255)
                                              cTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  1:1528896/100‑1 (MQ=255)
                                              cTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  2:723588/1‑100 (MQ=255)
                                                ggCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  2:2319991/1‑100 (MQ=255)
                                                ggCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:1424657/1‑100 (MQ=255)
                                                ggCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:989944/1‑100 (MQ=255)
                                                 gCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                <  1:1481749/100‑1 (MQ=255)
                                                   c‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  1:1539499/100‑1 (MQ=255)
                                                      ttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATc                                            <  2:2188048/100‑1 (MQ=255)
                                                       tCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCt                                           >  1:2051652/1‑100 (MQ=255)
                                                        cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                          <  1:1042791/100‑1 (MQ=255)
                                                          gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:984520/100‑1 (MQ=255)
                                                           aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       >  1:995826/1‑100 (MQ=255)
                                                            cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  1:1258826/1‑100 (MQ=255)
                                                            cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGACGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:157033/1‑100 (MQ=255)
                                                              gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  1:603202/100‑1 (MQ=255)
                                                                     acttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                             >  1:1812025/1‑100 (MQ=255)
                                                                      cttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  1:1843135/1‑100 (MQ=255)
                                                                      cttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  2:792970/1‑100 (MQ=255)
                                                                      cttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  2:1469356/1‑100 (MQ=255)
                                                                      cttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  2:2320462/1‑100 (MQ=255)
                                                                       ttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           <  2:2102471/100‑1 (MQ=255)
                                                                       tgcctCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           <  1:1827236/97‑1 (MQ=255)
                                                                        tactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  1:1880912/100‑1 (MQ=255)
                                                                         actCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  2:233808/100‑1 (MQ=255)
                                                                          ctCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:1439707/99‑1 (MQ=255)
                                                                          ctCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:24320/1‑100 (MQ=255)
                                                                          ctCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:2095614/1‑100 (MQ=255)
                                                                          ctCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATCTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         >  2:1112392/1‑99 (MQ=255)
                                                                           tCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  2:951730/100‑1 (MQ=255)
                                                                           tCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       >  1:969061/1‑100 (MQ=255)
                                                                            cgcgGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGGACACGTTCACGGGTGa                      >  2:252383/1‑100 (MQ=255)
                                                                            cgcgGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:1447872/1‑100 (MQ=255)
                                                                              cgGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:1760775/100‑1 (MQ=255)
                                                                                gAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  <  2:2066628/100‑1 (MQ=255)
                                                                                gAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGCGAGGcc                  >  1:1981401/1‑100 (MQ=255)
                                                                                   cAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               <  2:953662/100‑1 (MQ=255)
                                                                                    aGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  1:610037/1‑100 (MQ=255)
                                                                                      cgcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            >  1:558476/1‑100 (MQ=255)
                                                                                      cgcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            >  2:1380886/1‑100 (MQ=255)
                                                                                       gcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           >  2:1692441/1‑100 (MQ=255)
                                                                                       gcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  1:1219733/100‑1 (MQ=255)
                                                                                       gcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:889533/100‑1 (MQ=255)
                                                                                       gcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  1:876496/100‑1 (MQ=255)
                                                                                       gcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTCAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:1258821/100‑1 (MQ=255)
                                                                                        cTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1930839/100‑1 (MQ=255)
                                                                                        cTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:116154/100‑1 (MQ=255)
                                                                                          tCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  2:1604875/1‑100 (MQ=255)
                                                                                           cGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  2:1894203/1‑100 (MQ=255)
                                                                                           cGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:1497218/100‑1 (MQ=255)
                                                                                           cGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  2:347967/100‑1 (MQ=255)
                                                                                            gATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  2:419342/100‑1 (MQ=255)
                                                                                            gATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  1:562844/1‑100 (MQ=255)
                                                                                             atatTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAGACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCt     >  1:564258/1‑100 (MQ=255)
                                                                                              tatTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:1444336/100‑1 (MQ=255)
                                                                                              tatTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    >  2:35103/1‑100 (MQ=255)
                                                                                                ttCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  1:1318320/1‑100 (MQ=255)
                                                                                                ttCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:194553/1‑100 (MQ=255)
                                                                                                      |                                                                                               
CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCC‑TTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑CCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

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