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.

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682
                                                                                                    |                                                                                                
naTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                   >  2:3027992/2‑100 (MQ=255)
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                   >  2:1527358/1‑100 (MQ=255)
    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCACGGAACAGCGCTTCGCTATTCAGcccc                                                                                            >  2:1530507/1‑100 (MQ=255)
     gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccc                                                                                           <  2:3333545/100‑2 (MQ=255)
       ccTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccctg                                                                                         <  2:2944995/100‑4 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccctgc                                                                                        >  2:549758/1‑96 (MQ=255)
         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccctgcg                                                                                       <  2:2241258/100‑6 (MQ=255)
            tATTGTGCACAGAAAAGGCTAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGttt                                                                                     >  2:1610151/1‑100 (MQ=255)
              ttctgCACAGAAAAGGCCAGCCTGCCTTGAGACTGGCCTTTCTGACACATGCTTACTTACTCGCGGAACAGCTCTTCGATATTCAGCCCCCTGCGttt                                                                                     <  1:2924638/95‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAg                                                                                 <  2:3401524/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAg                                                                                 <  1:1498824/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAg                                                                                 >  2:2720939/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAgg                                                                                >  1:349384/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAgg                                                                                <  2:1781788/100‑1 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGTCCTTTCTGACAGATGCTTACTTACTCGCGGAACCGCGCTTCGATATTCAGCCCCCTGCGTTTGTAGGa                                                                               <  1:248517/100‑1 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGa                                                                               >  2:928636/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGa                                                                               >  2:3521224/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAt                                                                              >  2:1343230/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAt                                                                              >  1:2098016/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGACCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAt                                                                              >  1:2499894/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAttt                                                                            >  1:3286398/1‑100 (MQ=255)
                      aGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTc                                                                           >  1:2006999/1‑100 (MQ=255)
                        aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTcgc                                                                         <  2:2485810/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTcgcg                                                                        <  2:282767/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTcgcg                                                                        <  1:3274155/100‑1 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCgcgc                                                                       <  2:1526658/100‑1 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  2:655551/1‑100 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCACGATTTCGCGCaa                                                                     >  1:1881153/1‑100 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGATTCGATATTCCGCCCCCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:3624062/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:2606293/1‑100 (MQ=255)
                                 gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:2827238/1‑100 (MQ=255)
                                   cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  1:355025/100‑1 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             >  1:855360/1‑100 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  2:924341/100‑1 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  1:2154051/100‑1 (MQ=255)
                                    tCGCGTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  1:1002386/100‑1 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  2:1918564/100‑1 (MQ=255)
                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:2993013/100‑1 (MQ=255)
                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  1:495844/1‑100 (MQ=255)
                                          gagaCTGGCCTTTCTGACAGATGCTTACTTGCTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  1:667077/100‑1 (MQ=255)
                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  2:3447626/1‑100 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  2:2549879/100‑1 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  1:3447600/100‑1 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:1410697/100‑1 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  >  1:1683666/1‑100 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTGCTCGCGGAACAGCGCTTCGATATTCAGCCACGTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:1135550/99‑1 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:1477351/1‑100 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:623430/1‑100 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  2:2403249/1‑100 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  2:604356/1‑100 (MQ=255)
                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:2834567/1‑100 (MQ=255)
                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               <  2:753851/100‑1 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  2:2429905/1‑100 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  2:1056924/100‑1 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  1:621381/1‑100 (MQ=255)
                                                        tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  1:228107/1‑100 (MQ=255)
                                                        tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         <  1:3098618/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGTAACAGCGCGTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:3786726/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:3538773/1‑100 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:3375991/1‑100 (MQ=255)
                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  1:3587058/100‑1 (MQ=255)
                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:2294267/1‑100 (MQ=255)
                                                            aGATGCTTATTTACTTGTGGAACAGCGTTTCGTTATTCAGCTCCCTTCGTTTGCAGGATTTCGCTCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaa                                     <  1:1547589/100‑1 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  2:1090965/100‑1 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  1:3282174/100‑1 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   <  1:1108064/100‑1 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   <  1:1900968/100‑1 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   >  1:3645020/1‑100 (MQ=255)
                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   >  1:3648996/1‑100 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  <  2:235163/100‑1 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  2:1134952/1‑100 (MQ=255)
                                                                gCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                 >  2:1084957/1‑100 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  2:1324750/1‑100 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  2:2902651/1‑100 (MQ=255)
                                                                 attacttacTCGCGGAACAGCGCTTCGATATTCAGCCCCCTCCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:159655/2‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  2:1979514/100‑1 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  1:479762/100‑1 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:130719/1‑100 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:2680455/100‑1 (MQ=255)
                                                                              ggAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   >  1:3653912/1‑100 (MQ=255)
                                                                                aaCAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCACAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa                 <  1:876044/100‑1 (MQ=255)
                                                                                 aCAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  2:568001/1‑100 (MQ=255)
                                                                                  cAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               >  2:421651/1‑100 (MQ=255)
                                                                                   aGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              <  2:159655/100‑1 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:3785039/1‑100 (MQ=255)
                                                                                     cgcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCACTTTc            >  1:926706/1‑100 (MQ=255)
                                                                                      gcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:825838/100‑1 (MQ=255)
                                                                                      gcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:168242/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:3710103/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1860053/100‑1 (MQ=255)
                                                                                        ttCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAACTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2786733/99‑1 (MQ=255)
                                                                                         tCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        <  2:3486773/100‑1 (MQ=255)
                                                                                         tCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  1:1497152/1‑100 (MQ=255)
                                                                                         tCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  2:3120232/1‑100 (MQ=255)
                                                                                          cGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCTCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:2432352/100‑1 (MQ=255)
                                                                                            atatTCAGCCACCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCTTTCAACATGAATCTGGCGAAGACGTTCACGGGTGAAGACAATTTCACGACCt     >  1:3013868/1‑100 (MQ=255)
                                                                                             tatTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:1563714/100‑1 (MQ=255)
                                                                                             tatTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  2:3286422/100‑1 (MQ=255)
                                                                                               ttCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:3735005/1‑100 (MQ=255)
                                                                                               ttCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:1935676/1‑100 (MQ=255)
                                                                                                    |                                                                                                
CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  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: