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.C97.8% 476.1 / 3.1 135coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (63/68);  total (65/70)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682
                                                                                                     |                                                                                               
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                  <  1:1407587/100‑1 (MQ=255)
  tCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGcc                                                                                               <  2:3212092/100‑1 (MQ=255)
   cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGTCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGcc                                                                                               >  2:1387916/1‑99 (MQ=255)
   cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccc                                                                                              >  2:1953376/1‑100 (MQ=255)
    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGcccc                                                                                            >  2:1530766/1‑100 (MQ=255)
     gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccc                                                                                           >  1:1150736/1‑99 (MQ=255)
     gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccc                                                                                           >  1:2772784/1‑99 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccctgc                                                                                        >  1:2271530/1‑96 (MQ=255)
           ttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGtt                                                                                      <  2:1642606/100‑1 (MQ=255)
            tATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGttt                                                                                     >  1:581959/1‑100 (MQ=255)
             aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTg                                                                                    >  2:48037/1‑100 (MQ=255)
               tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCa                                                                                  >  1:2938420/1‑100 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAg                                                                                 <  2:1300863/100‑1 (MQ=255)
                 tgCACAGAAAAAACAAGCCTCGCTTGAGACCGGCCTTTCTGACAGACGATTAATTACTCACGAAAAAGCGCAACAATAATCAGCCCCCTGCGTTTGCAgg                                                                                >  2:1760746/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAgg                                                                                >  1:795547/1‑99 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGa                                                                               >  2:2118269/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGa                                                                               >  2:3182890/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCATCCCCCTGCGTTTGCAgg                                                                                >  1:474746/1‑98 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAt                                                                              <  1:2884204/100‑1 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAt                                                                              <  2:585636/100‑1 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAttt                                                                            >  1:684129/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGAAtt                                                                             >  1:942389/1‑99 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTcgcg                                                                        <  2:595572/100‑1 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  2:1173760/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   <  1:260907/100‑1 (MQ=255)
                                aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACgg                                                                 >  2:52733/1‑100 (MQ=255)
                                  ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               >  2:2541840/1‑100 (MQ=255)
                                   cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCATACGgcgc                                                              <  1:1709558/100‑1 (MQ=255)
                                   cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  1:1147376/100‑1 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  1:2090331/100‑1 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:1956453/100‑1 (MQ=255)
                                       cTTGAGACTGGCCTTTCTTACATACGCTTACTTTCTCGCGGAACTGCGCTTCGATATTCAGCCCTCTGCGTTTGCATGATTTCGCCCAAACGGCGCAgg                                                           <  2:2009089/99‑1 (MQ=255)
                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAggaga                                                        >  2:1348025/1‑97 (MQ=255)
                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:3108667/1‑100 (MQ=255)
                                          gagaCTGTCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  1:2606636/100‑1 (MQ=255)
                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCTCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      <  2:2071259/100‑1 (MQ=255)
                                            gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGACTTCa                                                     <  2:1388033/100‑1 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:2392795/100‑1 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                <  2:1586329/100‑1 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  2:116092/1‑100 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                <  1:1797585/100‑1 (MQ=255)
                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  1:2458749/1‑100 (MQ=255)
                                                     ttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATc                                            <  2:774228/100‑1 (MQ=255)
                                                       cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                          <  1:1749815/100‑1 (MQ=255)
                                                        tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         <  1:2452089/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  1:2696880/1‑100 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:1270810/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:2386733/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:3108695/100‑1 (MQ=255)
                                                          aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  1:669462/100‑1 (MQ=255)
                                                          aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:1171324/100‑1 (MQ=255)
                                                          aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  1:3198849/100‑1 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  1:839877/1‑100 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  <  1:2337618/100‑1 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  2:2355690/1‑100 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                 <  2:698096/99‑1 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  2:3055925/1‑100 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                <  2:148657/100‑1 (MQ=255)
                                                                   tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              <  2:280933/100‑1 (MQ=255)
                                                                   tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              <  1:1387924/100‑1 (MQ=255)
                                                                    ncttactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                             >  1:1221993/2‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          >  1:2647044/1‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  2:2591539/100‑1 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          >  2:2423478/1‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          >  1:1384627/1‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  2:1577409/100‑1 (MQ=255)
                                                                         ctCGCGGAACAGCGGTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGCACACGTTCACGGGt                        >  2:882418/1‑100 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:437886/1‑100 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:1673750/1‑100 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:34055/1‑100 (MQ=255)
                                                                          tCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  1:2118284/100‑1 (MQ=255)
                                                                          tCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       >  1:142879/1‑100 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:1318594/1‑100 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:505343/1‑100 (MQ=255)
                                                                            gcgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     <  1:806741/100‑1 (MQ=255)
                                                                            gcgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     >  1:386005/1‑100 (MQ=255)
                                                                               gAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  2:1698371/1‑100 (MQ=255)
                                                                                 aCAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACATTCAAGAAAGAGGCCaa                >  1:810399/1‑100 (MQ=255)
                                                                                  cAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               <  1:1398494/100‑1 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  2:865723/1‑100 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             <  2:2345642/100‑1 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  2:1819831/1‑100 (MQ=255)
                                                                                      gcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  1:34055/100‑1 (MQ=255)
                                                                                      gcTTAGATATTCAGTCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCATCCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:2856552/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1157700/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:505341/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1602686/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1173766/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2647917/100‑1 (MQ=255)
                                                                                        ttCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg         <  1:1953388/100‑1 (MQ=255)
                                                                                        ttCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg         <  1:533563/100‑1 (MQ=255)
                                                                                         tCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAcag        >  2:1157694/1‑98 (MQ=255)
                                                                                          cGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  1:1507613/1‑100 (MQ=255)
                                                                                          cGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  2:2746511/100‑1 (MQ=255)
                                                                                           gATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        <  2:1063179/98‑1 (MQ=255)
                                                                                           gATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  2:3138502/1‑100 (MQ=255)
                                                                                           gATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  1:411900/100‑1 (MQ=255)
                                                                                           gATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCAAATATCACGAcc      >  1:836627/1‑100 (MQ=255)
                                                                                            atatTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCt     >  1:176785/1‑100 (MQ=255)
                                                                                               ttCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  2:2973700/100‑1 (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 < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: