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.C94.7% 454.8 / 16.5 133coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (2/5);  new base C (70/56);  total (72/61)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682
                                                                                                    |                                                                                                
naTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                   >  2:3027992/2‑100 (MQ=255)
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                   >  2:1527358/1‑100 (MQ=255)
    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCACGGAACAGCGCTTCGCTATTCAG‑c‑c‑cc                                                                                            >  2:1530507/1‑100 (MQ=255)
     gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c‑ccc                                                                                           <  2:3333545/100‑2 (MQ=255)
       ccTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c‑ccctg                                                                                         <  2:2944995/100‑4 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c‑ccctgc                                                                                        >  2:549758/1‑96 (MQ=255)
         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c‑ccctgcg                                                                                       <  2:2241258/100‑6 (MQ=255)
            tATTGTGCACAGAAAAGGCTAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGttt                                                                                     >  2:1610151/1‑100 (MQ=255)
              ttctgCACAGAAAAGGCCAGCCTGCCTTGAGACTGGCCTTTCTGACACATGCTTACTTACTCGCGGAACAGCTCTTCGATATTCAGCC‑C‑CCTGCGttt                                                                                     <  1:2924638/95‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAg                                                                                 <  2:3401524/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAg                                                                                 <  1:1498824/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAg                                                                                 >  2:2720939/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAgg                                                                                >  1:349384/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAgg                                                                                <  2:1781788/100‑1 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGTCCTTTCTGACAGATGCTTACTTACTCGCGGAACCGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGTAGGa                                                                               <  1:248517/100‑1 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGa                                                                               >  2:928636/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGa                                                                               >  2:3521224/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAt                                                                              >  2:1343230/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAt                                                                              >  1:2098016/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGACCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAt                                                                              >  1:2499894/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAttt                                                                            >  1:3286398/1‑100 (MQ=255)
                      aGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTc                                                                           >  1:2006999/1‑100 (MQ=255)
                        aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTcgc                                                                         <  2:2485810/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTcgcg                                                                        <  2:282767/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTcgcg                                                                        <  1:3274155/100‑1 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCgcgc                                                                       <  2:1526658/100‑1 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  2:655551/1‑100 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCACGATTTCGCGCaa                                                                     >  1:1881153/1‑100 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGATTCGATATTCCGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:3624062/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:2606293/1‑100 (MQ=255)
                                 gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:2827238/1‑100 (MQ=255)
                                   cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  1:355025/100‑1 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             >  1:855360/1‑100 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  2:924341/100‑1 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  1:2154051/100‑1 (MQ=255)
                                    tCGCGTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  1:1002386/100‑1 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  2:1918564/100‑1 (MQ=255)
                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:2993013/100‑1 (MQ=255)
                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  1:495844/1‑100 (MQ=255)
                                          gagaCTGGCCTTTCTGACAGATGCTTACTTGCTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  1:667077/100‑1 (MQ=255)
                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  2:3447626/1‑100 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  2:2549879/100‑1 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  1:3447600/100‑1 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:1410697/100‑1 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  >  1:1683666/1‑100 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTGCTCGCGGAACAGCGCTTCGATATTCAG‑C‑CACGTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:1135550/99‑1 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:1477351/1‑100 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:623430/1‑100 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  2:2403249/1‑100 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  2:604356/1‑100 (MQ=255)
                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:2834567/1‑100 (MQ=255)
                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               <  2:753851/100‑1 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  2:2429905/1‑100 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  2:1056924/100‑1 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  1:621381/1‑100 (MQ=255)
                                                        tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  1:228107/1‑100 (MQ=255)
                                                        tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         <  1:3098618/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGTAACAGCGCGTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:3786726/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:3538773/1‑100 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:3375991/1‑100 (MQ=255)
                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  1:3587058/100‑1 (MQ=255)
                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:2294267/1‑100 (MQ=255)
                                                            aGATGCTTATTTACTTGTGGAACAGCGTTTCGTTATTCAG‑CTC‑CCTTCGTTTGCAGGATTTCGCTCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaa                                     <  1:1547589/100‑1 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  2:1090965/100‑1 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  1:3282174/100‑1 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   <  1:1108064/100‑1 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   <  1:1900968/100‑1 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   >  1:3645020/1‑100 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   >  1:3648996/1‑100 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  <  2:235163/100‑1 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  2:1134952/1‑100 (MQ=255)
                                                                gCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                 >  2:1084957/1‑100 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  2:1324750/1‑100 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  2:2902651/1‑100 (MQ=255)
                                                                 attacttacTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTCCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:159655/2‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  2:1979514/100‑1 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  1:479762/100‑1 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:130719/1‑100 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:2680455/100‑1 (MQ=255)
                                                                              ggAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   >  1:3653912/1‑100 (MQ=255)
                                                                                aaCAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCACAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa                 <  1:876044/100‑1 (MQ=255)
                                                                                 aCAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  2:568001/1‑100 (MQ=255)
                                                                                  cAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               >  2:421651/1‑100 (MQ=255)
                                                                                   aGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              <  2:159655/100‑1 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:3785039/1‑100 (MQ=255)
                                                                                     cgcTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCACTTTc            >  1:926706/1‑100 (MQ=255)
                                                                                      gcTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:825838/100‑1 (MQ=255)
                                                                                      gcTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:168242/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:3710103/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1860053/100‑1 (MQ=255)
                                                                                        ttCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAACTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2786733/99‑1 (MQ=255)
                                                                                         tCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        <  2:3486773/100‑1 (MQ=255)
                                                                                         tCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  1:1497152/1‑100 (MQ=255)
                                                                                         tCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  2:3120232/1‑100 (MQ=255)
                                                                                          cGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCTCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:2432352/100‑1 (MQ=255)
                                                                                            atatTCAG‑C‑CACCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCTTTCAACATGAATCTGGCGAAGACGTTCACGGGTGAAGACAATTTCACGACCt     >  1:3013868/1‑100 (MQ=255)
                                                                                             tatTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:1563714/100‑1 (MQ=255)
                                                                                             tatTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  2:3286422/100‑1 (MQ=255)
                                                                                               ttCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:3735005/1‑100 (MQ=255)
                                                                                               ttCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:1935676/1‑100 (MQ=255)
                                                                                                    |                                                                                                
CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

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