Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,866,592 +C 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.C97.1% 374.9 / 3.9 104coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (57/44);  total (59/45)
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.99e-01

Only 100 of 109 total aligned reads displayed.

ATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑CCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866490‑2866682
                                                                                                    |                                                                                               
aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑c                                                                                                >  2:2706275/1‑100 (MQ=255)
aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑c                                                                                                <  2:515971/100‑1 (MQ=255)
    gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑ccc‑cc                                                                                           <  1:1041778/100‑2 (MQ=255)
    gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑ccc‑cc                                                                                           >  2:484836/1‑99 (MQ=255)
        ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑ccc‑cctgcg                                                                                       >  1:1225745/1‑95 (MQ=255)
          ttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGtt                                                                                      <  1:1738435/100‑1 (MQ=255)
            aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTg                                                                                    >  1:620696/1‑100 (MQ=255)
             ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGc                                                                                   <  1:1266067/100‑1 (MQ=255)
               gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAg                                                                                 <  1:348671/100‑1 (MQ=255)
               gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAg                                                                                 >  2:1031574/1‑100 (MQ=255)
               gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAg                                                                                 >  2:1279398/1‑100 (MQ=255)
               gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAg                                                                                 >  1:1998611/1‑100 (MQ=255)
                tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGTACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAgg                                                                                <  1:1942201/100‑1 (MQ=255)
                  cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAt                                                                              <  2:1313893/100‑1 (MQ=255)
                   acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAtt                                                                             >  1:334335/1‑100 (MQ=255)
                   acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAtt                                                                             >  1:941460/1‑100 (MQ=255)
                   acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAtt                                                                             >  1:127210/1‑100 (MQ=255)
                    caGAAAAGNCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCNGGAttt                                                                            >  1:2521215/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAttt                                                                            >  1:1409305/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAttt                                                                            >  2:2255135/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGAttt                                                                            >  1:1054107/1‑100 (MQ=255)
                      gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTc                                                                           >  2:830959/1‑99 (MQ=255)
                       aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTcgc                                                                         <  1:606610/100‑1 (MQ=255)
                         aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCgcgc                                                                       <  2:2519279/100‑1 (MQ=255)
                           ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCaa                                                                     <  1:1127235/100‑1 (MQ=255)
                             ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:1262634/1‑100 (MQ=255)
                             ccAGCCTCGCTTGAGACTGGCCTTTATGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCaaa                                                                    >  1:2721143/1‑99 (MQ=255)
                              cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACg                                                                  <  2:1998594/100‑1 (MQ=255)
                              cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACg                                                                  >  2:2261554/1‑100 (MQ=255)
                              cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACg                                                                  >  2:1176618/1‑100 (MQ=255)
                               aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACgg                                                                 <  1:828734/100‑1 (MQ=255)
                               aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACgg                                                                 >  1:2656402/1‑100 (MQ=255)
                                gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                <  1:2336835/100‑1 (MQ=255)
                                  cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              >  1:606827/1‑100 (MQ=255)
                                   tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  1:511680/100‑1 (MQ=255)
                                    cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             >  2:1327266/1‑99 (MQ=255)
                                    cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  1:704905/100‑1 (MQ=255)
                                     gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             >  2:2197146/1‑98 (MQ=255)
                                     gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  2:295406/100‑1 (MQ=255)
                                     gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  2:1077804/1‑100 (MQ=255)
                                      cTTGAGACTGGCCTTTCTGCCAGATGCTTACTTCCTCGCGGAACAGCGCTTCGATACTCACA‑CCC‑CCGCGTTCGCACGATCTCGCGCACACgccccacgc                                                          >  2:2630542/1‑92 (MQ=255)
                                         gcgaCTGGCCTTTCTGACAGATGCTTAGTTACTCGCGGAACAGCGCTGCGATATTGA‑G‑CCCACTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  2:231063/96‑1 (MQ=255)
                                         gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       >  2:1355155/1‑100 (MQ=255)
                                          agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      <  1:1279405/100‑1 (MQ=255)
                                           gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  2:2747697/1‑100 (MQ=255)
                                             cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  1:1344413/1‑100 (MQ=255)
                                               ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  2:2002656/1‑100 (MQ=255)
                                                 ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:1986922/1‑100 (MQ=255)
                                                 ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:1847175/1‑100 (MQ=255)
                                                  cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  2:1008366/1‑100 (MQ=255)
                                                      cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATc                                            <  2:1675119/98‑1 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:2606900/1‑100 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  1:2584309/1‑100 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:510779/1‑100 (MQ=255)
                                                        gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  2:2568125/100‑1 (MQ=255)
                                                        gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:1372315/1‑100 (MQ=255)
                                                        gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATCCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:305314/100‑1 (MQ=255)
                                                         aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:1752701/100‑1 (MQ=255)
                                                         aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:805305/100‑1 (MQ=255)
                                                          ttccgatctacttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  2:560920/92‑1 (MQ=255)
                                                          cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:711481/1‑100 (MQ=255)
                                                          cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  1:484840/100‑1 (MQ=255)
                                                          cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:1457787/1‑100 (MQ=255)
                                                           aGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaa                                     <  2:1840994/100‑1 (MQ=255)
                                                           aGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaa                                     >  1:1676618/1‑100 (MQ=255)
                                                            gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  2:2295746/100‑1 (MQ=255)
                                                            gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  2:1352721/1‑100 (MQ=255)
                                                            gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  2:199849/100‑1 (MQ=255)
                                                              tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  2:149133/1‑100 (MQ=255)
                                                               gCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                 <  2:1103025/100‑1 (MQ=255)
                                                                cttacttaCTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                <  1:2706294/100‑1 (MQ=255)
                                                                 ttacttacTCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGtt                               <  2:2562453/100‑1 (MQ=255)
                                                                  tacttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              <  1:2002673/100‑1 (MQ=255)
                                                                   acttactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaagatcggaa                             >  1:560924/1‑92 (MQ=255)
                                                                    cttactCGCGGAACAGCGCTTCGATATTCA‑G‑CCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           >  2:660355/1‑100 (MQ=255)
                                                                      tactCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          >  1:170641/1‑100 (MQ=255)
                                                                        ctCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCATGGGt                        >  2:700186/1‑100 (MQ=255)
                                                                        ctCGCGGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:1145636/1‑100 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCA‑GCCCC‑CTGTGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:2685165/1‑100 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAATACGTTCACGGGTg                       >  2:567469/1‑99 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:1221501/1‑100 (MQ=255)
                                                                            ctcttccgatcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:566071/90‑1 (MQ=255)
                                                                            cgGAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:410880/100‑1 (MQ=255)
                                                                             ggAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  1:2240562/100‑1 (MQ=255)
                                                                             ggAACAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   >  1:2075837/1‑100 (MQ=255)
                                                                                 cAGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               >  2:1075823/1‑100 (MQ=255)
                                                                                  aGCGCTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCATGGGTGAGGCCAAtt              <  1:700191/100‑1 (MQ=255)
                                                                                   gcgcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             <  2:1087137/100‑1 (MQ=255)
                                                                                   gcgcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  2:2009166/1‑100 (MQ=255)
                                                                                    cgcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            >  1:2415190/1‑100 (MQ=255)
                                                                                    cgcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            <  1:2179223/100‑1 (MQ=255)
                                                                                     gcTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTcn           <  2:1253169/100‑2 (MQ=255)
                                                                                      cTTCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  2:1124405/1‑100 (MQ=255)
                                                                                       ttCGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg         <  2:855034/100‑1 (MQ=255)
                                                                                         cGATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  2:2137534/1‑100 (MQ=255)
                                                                                          gATATTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  1:1077811/100‑1 (MQ=255)
                                                                                            tatTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:1221508/100‑1 (MQ=255)
                                                                                            tatTCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:1143270/100‑1 (MQ=255)
                                                                                              ttCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  1:808179/1‑100 (MQ=255)
                                                                                              ttCA‑GCCCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  1:1750069/1‑100 (MQ=255)
                                                                                                    |                                                                                               
ATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCA‑G‑CCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866490‑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: ‑