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.C90.3% 351.2 / 31.5 113coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (3/8);  new base C (51/49);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

ATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866490‑2866682
                                                                                                   |                                                                                                
aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCTGAACAGCGCTTCGATATTCAG‑c                                                                                                 <  1:54118/100‑1 (MQ=255)
  cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑c                                                                                              <  2:58652/100‑1 (MQ=255)
    gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccc                                                                                            >  2:1287956/1‑99 (MQ=255)
       cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgc                                                                                         >  1:2142861/1‑96 (MQ=255)
       cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgc                                                                                         >  2:2517279/1‑96 (MQ=255)
        ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgcg                                                                                        <  1:1271792/100‑6 (MQ=255)
        ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgcg                                                                                        >  2:2450144/1‑95 (MQ=255)
        ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccc                                                                                            >  1:666315/1‑95 (MQ=255)
          ttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGtt                                                                                      <  2:116067/100‑1 (MQ=255)
            aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGGCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTg                                                                                    <  2:131099/100‑1 (MQ=255)
            aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGtttn                                                                                    <  1:2077187/100‑2 (MQ=255)
            aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTg                                                                                    <  2:2585313/100‑1 (MQ=255)
             ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACAGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGc                                                                                   >  1:1926498/1‑100 (MQ=255)
              tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCa                                                                                  <  2:2682134/100‑1 (MQ=255)
              tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCa                                                                                  <  1:1100458/100‑1 (MQ=255)
              tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCa                                                                                  <  2:1633610/100‑1 (MQ=255)
               gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAg                                                                                 <  1:763402/100‑1 (MQ=255)
                 gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGa                                                                               >  2:1587284/1‑100 (MQ=255)
                   acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGATCAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGcgcgaaa                                                                             >  1:1582221/1‑94 (MQ=255)
                   acaGAAAAGGCCAGCCTCGCGTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGAtt                                                                             >  1:1780758/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGCTATTCAG‑CC‑CC‑TGCGTTTGCAGGATTTc                                                                           >  2:1371684/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGAttt                                                                            >  1:990094/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTGTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGAttt                                                                            <  1:1203169/100‑1 (MQ=255)
                      gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGCTTTcg                                                                          >  1:1884045/1‑100 (MQ=255)
                      gAAAAGGCCAGCCTCGCGTGAGACTGGCCTTTCTGACAGATGCTGACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTcg                                                                          <  2:2716112/100‑1 (MQ=255)
                        aatgtccAGCCTCGCTTGAAACTTGCCTTTCTGTCAGATGCTTACTTACTCTCGGAACAGCTCTTCTATATTCAG‑CCTCC‑TGCGTTTGCAGGATTTcg                                                                          <  1:2501631/93‑1 (MQ=255)
                         aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCgcgc                                                                       <  2:1114716/100‑1 (MQ=255)
                         aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCgcgc                                                                       <  1:1784930/100‑1 (MQ=255)
                         aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCgcgc                                                                       <  1:1968508/100‑1 (MQ=255)
                           ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:1133020/1‑100 (MQ=255)
                           ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCaa                                                                     >  2:1782064/1‑100 (MQ=255)
                           ggCCAGCCTCGCTTGAAACTGGCCTTTCTGACTGATGCTTACTTAATCGCGGAACTGCCCTTCGATATTCTG‑CC‑CCATGCGTTTGCAGGATTTCgcgc                                                                       <  1:2794954/98‑1 (MQ=255)
                             ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAAc                                                                   <  2:58218/100‑1 (MQ=255)
                               aGCCTCGCTTGAGAGTGGCCTTTCTGACAGATACTTACTTACTTGCGGAACAGCGCTTCGATATTCAGACC‑CC‑TGCATTTCGAGGATTTCGCGCAAACgg                                                                 >  1:1863235/1‑100 (MQ=255)
                               aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACgg                                                                 <  2:84992/100‑1 (MQ=255)
                                gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:1533701/1‑100 (MQ=255)
                                gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:1053413/1‑100 (MQ=255)
                                gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:918750/1‑100 (MQ=255)
                                 ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               <  2:1857938/100‑1 (MQ=255)
                                     gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:1492551/100‑1 (MQ=255)
                                      cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  1:590012/1‑100 (MQ=255)
                                      cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  1:747597/1‑100 (MQ=255)
                                       ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  2:2477501/100‑1 (MQ=255)
                                        tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcat                                                        >  2:202807/1‑98 (MQ=255)
                                        tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:1638632/1‑100 (MQ=255)
                                        tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:1638685/1‑100 (MQ=255)
                                        tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCCCAAACGGCGCAGGCCt                                                        >  1:2172669/1‑100 (MQ=255)
                                           gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  2:2142874/100‑1 (MQ=255)
                                            aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    >  2:2288626/1‑100 (MQ=255)
                                             cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  1:1623448/1‑100 (MQ=255)
                                               ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:1790370/1‑100 (MQ=255)
                                               ggCCCTGCTGACAGATGCTCACTTCCTCGCCGAGCAGCGCGTCGATATTCAG‑CC‑CCGTGCGTTTGCAGGATTTCTCGCAAATGGCGCAGGCCGTTAACCt                                                 <  1:1437117/100‑1 (MQ=255)
                                                gCCGTTCTGACAGATGCTGACGTGCTGGGGTAACGGCGCTTCGATATTCAG‑CCGCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    <  1:2602666/96‑1 (MQ=255)
                                                 ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  1:2629801/1‑100 (MQ=255)
                                                   tttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAt                                             <  2:220873/100‑1 (MQ=255)
                                                    ttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATc                                            >  1:1394253/1‑100 (MQ=255)
                                                     tCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCt                                           >  1:2230657/1‑100 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:2650945/1‑100 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCAATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:1552530/1‑100 (MQ=255)
                                                        gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  1:732970/1‑100 (MQ=255)
                                                        gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:559339/1‑100 (MQ=255)
                                                         aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:870335/100‑1 (MQ=255)
                                                              tGCTTACTTACTCGCGGACAAGCGCTTCGATATTCAGACC‑TC‑TGCGTTTGCAGGATTTCGCGCAAGCGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  1:1320310/98‑1 (MQ=255)
                                                                cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:2462080/1‑100 (MQ=255)
                                                                cttacttaCTCGCGGAACAGCGCTTCGATATTCAG‑CC‑CCATGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAAAACGt                                >  1:1082983/1‑100 (MQ=255)
                                                                   acttactCGCGGCACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCTCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                <  1:967260/97‑1 (MQ=255)
                                                                    cttactCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  2:264279/1‑100 (MQ=255)
                                                                     ttactCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           <  1:1053410/100‑1 (MQ=255)
                                                                      tactCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  1:318957/100‑1 (MQ=255)
                                                                        ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCCGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:249533/1‑100 (MQ=255)
                                                                        ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:2525542/1‑100 (MQ=255)
                                                                        ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:271128/1‑100 (MQ=255)
                                                                        ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:2311074/1‑100 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:642322/1‑100 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:220163/1‑100 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:921584/1‑100 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCAG‑CCTCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:740045/1‑100 (MQ=255)
                                                                           gcgGAACAGCGCTTCGATATTCAG‑CC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:1601785/100‑1 (MQ=255)
                                                                            cgGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgt                    <  1:642319/100‑2 (MQ=255)
                                                                            cgGAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:1447193/100‑1 (MQ=255)
                                                                            cgGAACAGCGCTTCGATATTCAG‑CCGCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:1322412/100‑1 (MQ=255)
                                                                            cgGAACAGCGCTTCGATATTCAG‑CC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  1:2522678/100‑1 (MQ=255)
                                                                              gAACAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  1:1980380/1‑100 (MQ=255)
                                                                               aaCAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa                 >  2:80447/1‑100 (MQ=255)
                                                                                aCAGCGCTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                <  2:2311089/100‑1 (MQ=255)
                                                                                   gcgcTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             <  1:918747/100‑1 (MQ=255)
                                                                                   gcgcTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:85049/1‑100 (MQ=255)
                                                                                    cgcTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            >  1:1997551/1‑100 (MQ=255)
                                                                                     gcTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:1830950/100‑1 (MQ=255)
                                                                                     gcTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           >  2:451088/1‑100 (MQ=255)
                                                                                     gcTTCGATATTCAG‑CC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:453429/100‑1 (MQ=255)
                                                                                      cTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:1943838/100‑1 (MQ=255)
                                                                                      cTTCGATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1453665/100‑1 (MQ=255)
                                                                                          gATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  2:1251611/1‑100 (MQ=255)
                                                                                          gATATTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  1:204805/1‑100 (MQ=255)
                                                                                           atatTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCATCCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCt     <  1:161083/100‑1 (MQ=255)
                                                                                            tatTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAGGGGTGAGGCCAATTTCACGACCt     <  1:355604/99‑1 (MQ=255)
                                                                                            tatTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  1:2442965/98‑1 (MQ=255)
                                                                                            tatTCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  2:2883557/100‑1 (MQ=255)
                                                                                              ttCAGCCC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:1538903/1‑100 (MQ=255)
                                                                                                   |                                                                                                
ATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CC‑CC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866490‑2866682

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