Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,822,869 G→A 18.2% R157H (CGC→CAC)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,822,8690GA18.2% 346.2 / 58.4 157R157H (CGC→CAC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (70/57);  new base A (18/11);  total (88/69)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 175 total aligned reads displayed.

TGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCT  >  NC_000913/3822770‑3822968
                                                                                                   |                                                                                                   
tGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTcgccg                                                                                                     <  2:3549064/100‑1 (MQ=255)
 ggCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTcgccgc                                                                                                    <  1:1727949/100‑1 (MQ=255)
  gCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTcgccgc                                                                                                    >  2:3008014/1‑99 (MQ=255)
  gCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCcaca                                                                                                   >  1:121829/1‑100 (MQ=255)
  gCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCa                                                                                                   >  1:83664/1‑100 (MQ=255)
  gCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCa                                                                                                   >  1:2079537/1‑100 (MQ=255)
       ggTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGcc                                                                                              <  2:3655362/100‑1 (MQ=255)
        gTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGccc                                                                                             >  1:3374402/1‑100 (MQ=255)
        gTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGccc                                                                                             >  1:2515029/1‑100 (MQ=255)
         tccatgagatCCGCGTCATCCTCATCATACTTGCCCACCTTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCg                                                                                            <  2:3654138/92‑1 (MQ=255)
         tGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCACATCGCCCg                                                                                            >  2:3033177/1‑100 (MQ=255)
         tGCAGGATATCCGCGACATCATCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGccc                                                                                             >  2:911798/1‑99 (MQ=255)
          gCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCTACAACATGCGCACGCTGGGCTCACTTAGACCGGACAAACGTCGCCACATCGCCCGt                                                                                           >  1:2119746/1‑100 (MQ=255)
          gCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGcc                                                                                              >  2:352577/1‑97 (MQ=255)
          gCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGt                                                                                           >  2:233222/1‑100 (MQ=255)
           cAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCACATCGCCCGTg                                                                                          <  2:1879838/100‑1 (MQ=255)
            aGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATACGCACTCTAGACTCACTACGCCCGGACAAACGTCGCCACATCGCCCg                                                                                            >  1:572973/1‑97 (MQ=255)
             ggATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGACTCAATTCGCCCCAACAAAAGTCGCCACATCTcccataaa                                                                                        >  1:3733734/1‑95 (MQ=255)
               atatCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAAc                                                                                      <  1:1350312/100‑1 (MQ=255)
                  tCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAActct                                                                                   >  1:2102319/1‑100 (MQ=255)
                   ccGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACtctc                                                                                  >  2:1423157/1‑100 (MQ=255)
                    cgcgTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCg                                                                                 <  1:3513212/100‑1 (MQ=255)
                        tCATCCTCATCAAACTTGCCGACCGTACCTACTACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAAt                                                                             <  1:975013/100‑1 (MQ=255)
                          aTCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAAttt                                                                           >  2:2077923/1‑100 (MQ=255)
                               catcaAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGc                                                                      <  2:1322076/100‑1 (MQ=255)
                               catcaAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGc                                                                      >  2:1062241/1‑100 (MQ=255)
                                  caAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCg                                                                   >  2:2905228/1‑100 (MQ=255)
                                   aaaCTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGc                                                                  >  1:153947/1‑100 (MQ=255)
                                    aaCTTGCCGACCGTACCCACAACATGAGAACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGTCAGTGAAACTCTCGAAATTTATAGCCCGc                                                                  >  2:2110375/1‑99 (MQ=255)
                                     aCTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTg                                                                >  2:722412/1‑100 (MQ=255)
                                       ttGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGgc                                                              <  1:2318478/100‑1 (MQ=255)
                                         gCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCACGCTAgcg                                                             >  1:927626/1‑99 (MQ=255)
                                            gACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCAc                                                          <  2:1182084/99‑1 (MQ=255)
                                             aCCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCg                                                        >  2:3434885/1‑100 (MQ=255)
                                             aCCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCg                                                        >  2:36756/1‑100 (MQ=255)
                                                 tACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTa                                                    <  2:1974075/100‑1 (MQ=255)
                                                 tACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTa                                                    <  1:3434924/100‑1 (MQ=255)
                                                   cccACAACATGCGCACGCTGGGCTCCCTTCGCCCGGACAAACGTCGCCACATCGCCCGTGCAACTCTCGAAATTTATAGCCCGCTGGCGCACAGTTTAgg                                                  >  1:1434220/1‑100 (MQ=255)
                                                   cccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAgn                                                  <  1:3811812/100‑2 (MQ=255)
                                                   cccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAgg                                                  >  1:71998/1‑100 (MQ=255)
                                                   cccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAgg                                                  <  2:1932764/100‑1 (MQ=255)
                                                   cccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAgg                                                  <  1:1912226/100‑1 (MQ=255)
                                                   cccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAgg                                                  <  2:203529/100‑1 (MQ=255)
                                                   cccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAgg                                                  >  2:2004425/1‑100 (MQ=255)
                                                    ccACAACATGCGCACGCTGGGCTTACTTCGTCTGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAgg                                                  <  1:286124/99‑1 (MQ=255)
                                                    ccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGt                                                 >  1:2370780/1‑100 (MQ=255)
                                                     cacaACATGCTCACGCTGGGCTCACTTCGCTCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTa                                                <  1:674334/100‑1 (MQ=255)
                                                     cacaACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGTATTTTATAGCCCGCTGGCGCACCGTTTAGGTa                                                <  2:2484001/100‑1 (MQ=255)
                                                     cacaACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCATGAAACTCTAGAAATTTATAGCCCGCTGGCGCAACGTTTAGGTa                                                >  1:2282010/1‑100 (MQ=255)
                                                     cacaACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTa                                                >  1:1024665/1‑100 (MQ=255)
                                                      acaacaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTAt                                               >  1:1884758/1‑100 (MQ=255)
                                                       caacaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATc                                              >  2:1773665/1‑100 (MQ=255)
                                                         acaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATcca                                            <  1:3033204/100‑1 (MQ=255)
                                                           aTGCGCACGCTGGGCTCACTTCGCCCGGATAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATccacc                                          <  1:1383028/100‑1 (MQ=255)
                                                           aTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATccacn                                          <  2:3307502/100‑2 (MQ=255)
                                                           aTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATccacc                                          <  1:2002843/100‑1 (MQ=255)
                                                              cgcACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACcaca                                       >  2:76053/1‑100 (MQ=255)
                                                                cACGCTTCGCTCACTTCTCCCGGACAAACGTCGCCGCATCGCCCGTTAAACTCTCGAAATTTATAGCCCGCTTGCGCACCGTTTAGGTATccacc                                          <  2:3759644/95‑1 (MQ=255)
                                                                 aCGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAATTTTATAGCCCTCTGGCGCACCGTTTAGGTATCCACCACATTa                                    <  1:2644177/100‑1 (MQ=255)
                                                                  cGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTaa                                   >  1:3067701/1‑100 (MQ=255)
                                                                  cGCTGGGCTCACTTCGCCCGGACAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTaa                                   >  1:1617449/1‑100 (MQ=255)
                                                                   gCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTaaa                                  <  1:2583367/100‑1 (MQ=255)
                                                                    cTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTaaaa                                 <  1:251863/100‑1 (MQ=255)
                                                                     tGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCAGATAGCGCACCGTTTAGGTATCCACCACATTAAAAc                                >  1:2456636/1‑100 (MQ=255)
                                                                      gggCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAAcc                               >  2:1166154/1‑100 (MQ=255)
                                                                      gggCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAAc                                <  2:2322303/99‑1 (MQ=255)
                                                                       ggCTCACTGCGCCCGGACGAGCGTCTCCTGATCGCCCGTGAAGCTCTTGAAATGTATATCCCGCTGGCGCATCGTTCAGGTATCCACCACATTaaggcc                               <  2:2269057/99‑5 (MQ=255)
                                                                        gCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGa                             <  1:3762048/100‑1 (MQ=255)
                                                                        gCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTATAGAAATTTATAGCCCGATGGCGCACCGTTTAGGTATCCACCACATTAAAAACGa                             >  1:2335000/1‑100 (MQ=255)
                                                                         cTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGaa                            >  2:3703845/1‑100 (MQ=255)
                                                                          tCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAAc                           <  1:2876643/100‑1 (MQ=255)
                                                                           cACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACt                          >  2:2559983/1‑100 (MQ=255)
                                                                           cACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACt                          >  1:2350755/1‑100 (MQ=255)
                                                                             cTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCg                        >  2:226888/1‑100 (MQ=255)
                                                                             cTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCg                        <  2:3374364/100‑1 (MQ=255)
                                                                              ttCGCCCGGACAAACGTCGCCGCATCGCCTGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGa                       <  1:3661705/100‑1 (MQ=255)
                                                                              ttCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTAGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGa                       >  1:3317635/1‑100 (MQ=255)
                                                                               tCGCCCGGACAAACGTCGCCGCATCACCCGTGATTCTCTCGAAATTTATATCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACGCGaa                      <  2:3400302/100‑1 (MQ=255)
                                                                                cGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAag                     >  2:841404/1‑100 (MQ=255)
                                                                                cGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAag                     >  1:50714/1‑100 (MQ=255)
                                                                                    cGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCt                 >  2:3107922/1‑100 (MQ=255)
                                                                                     ggACAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTg                <  1:47369/100‑1 (MQ=255)
                                                                                       aCAAACGTCGCCGTATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTggg              <  1:892202/100‑1 (MQ=255)
                                                                                        cAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGt             >  1:536241/1‑100 (MQ=255)
                                                                                        cAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGt             >  1:2705217/1‑100 (MQ=255)
                                                                                        cAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTggg              >  2:3081936/1‑99 (MQ=255)
                                                                                        cAAACGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGAACCGATTAAGTAACCACCACATTAAAACCGAACTCGAAGAGCt                 >  1:43328/1‑96 (MQ=255)
                                                                                         aaaCGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGtt            <  2:1545354/100‑1 (MQ=255)
                                                                                          aaCGTCGCCACATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGttt           >  2:1021030/1‑100 (MQ=255)
                                                                                            cGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTg         >  1:2933812/1‑100 (MQ=255)
                                                                                               cgtcgcATCGACAGTGAAAATCTCGAAATTAATAGCCCGCTGGCACACCGTTTAGGTATCCACAACATTAAAACCGAACTCGAAGAGATggg              >  2:1426108/4‑92 (MQ=255)
                                                                                               cgccgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGa        >  2:584213/1‑98 (MQ=255)
                                                                                               cgccgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAgg      <  2:2932318/100‑1 (MQ=255)
                                                                                               cgccgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAgg      <  2:2932613/100‑1 (MQ=255)
                                                                                                gccgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGTCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGgc     >  2:442059/1‑100 (MQ=255)
                                                                                                gccgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGgc     >  1:3465324/1‑100 (MQ=255)
                                                                                                gccgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTAGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGgc     <  2:2804218/100‑1 (MQ=255)
                                                                                                gccgcAGCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAgg      <  2:1224985/99‑1 (MQ=255)
                                                                                                gccacATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAgg      <  2:1851968/99‑1 (MQ=255)
                                                                                                  cgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAAgcgc   >  2:556924/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCT  >  NC_000913/3822770‑3822968

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