Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,822,5540TC22.8% 405.7 / 121.0 219V52A (GTT→GCT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (85/84);  new base C (25/25);  total (110/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 220 total aligned reads displayed.

AAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGA  >  NC_000913/3822455‑3822650
                                                                                                   |                                                                                                
  aTCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCGCAAGAGGGGCAAACACGATCAAGCGGTGAACCCTATATCCCGCACCCGGTAGCGGTTg                                                                                                >  2:1767715/1‑100 (MQ=255)
   tCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGc                                                                                               >  2:1116785/1‑100 (MQ=255)
    cAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGcc                                                                                              <  2:3092238/100‑1 (MQ=255)
    cAAGCGTCTGCGGCAGGCGTATCACGTTGCACGTGATGCTCAGGAGGTGCAAAAACGTTCAAGCGGTGAACCATATATCACGCACCCGGTAGCGGTTGcc                                                                                              >  2:1894583/1‑100 (MQ=255)
      aGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTg                                                                                            >  1:3271189/1‑100 (MQ=255)
              cGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACTCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTgg                                                                                    <  1:167314/100‑1 (MQ=255)
              cGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTgg                                                                                    <  1:2853733/100‑1 (MQ=255)
              cGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTgg                                                                                    <  2:1359193/100‑1 (MQ=255)
                gCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGcc                                                                                  <  2:2588591/100‑1 (MQ=255)
                gCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGcc                                                                                  <  2:899789/100‑1 (MQ=255)
                gCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGcc                                                                                  >  1:3761639/1‑100 (MQ=255)
                gCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGcc                                                                                  <  2:720575/100‑1 (MQ=255)
                 cAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCg                                                                                 <  2:1339753/100‑1 (MQ=255)
                   ggCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCgag                                                                               >  1:4152831/1‑100 (MQ=255)
                   ggCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTCCATTCTGGCCgcg                                                                               >  2:4152028/1‑98 (MQ=255)
                   ggCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTCGCCTGCATTCTGGCCgag                                                                               <  2:5038667/100‑1 (MQ=255)
                   ggCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCgag                                                                               <  1:63953/100‑1 (MQ=255)
                     cGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGAt                                                                             >  2:3470332/1‑100 (MQ=255)
                     cGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGAt                                                                             <  1:1943588/100‑1 (MQ=255)
                       tATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGa                                                                           >  2:2779462/1‑100 (MQ=255)
                        aTCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGaa                                                                          >  1:1021559/1‑100 (MQ=255)
                        aTCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGaa                                                                          >  2:3360670/1‑100 (MQ=255)
                          ctcGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCGATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAAc                                                                        >  2:4605370/1‑100 (MQ=255)
                            cGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTc                                                                      >  2:1262600/1‑100 (MQ=255)
                             gTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCg                                                                     <  1:3398284/100‑1 (MQ=255)
                                 cACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTa                                                                 >  2:2155993/1‑100 (MQ=255)
                                 cACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTa                                                                 >  1:4705470/1‑100 (MQ=255)
                                 cACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTa                                                                 >  2:3286865/1‑100 (MQ=255)
                                  aCGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTAt                                                                <  2:3245948/100‑1 (MQ=255)
                                   cGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATg                                                               >  1:3792973/1‑100 (MQ=255)
                                    gTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGa                                                              <  1:1793967/100‑1 (MQ=255)
                                    gTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGa                                                              >  2:3664204/1‑100 (MQ=255)
                                    gTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGa                                                              >  1:4585126/1‑100 (MQ=255)
                                     tGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGaa                                                             <  1:599781/100‑1 (MQ=255)
                                     tGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGaa                                                             >  2:1884189/1‑100 (MQ=255)
                                      gATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGaaa                                                            <  1:4261260/100‑1 (MQ=255)
                                        tGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACg                                                          >  1:1046704/1‑100 (MQ=255)
                                         gCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGc                                                         <  1:4199818/100‑1 (MQ=255)
                                           tCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTg                                                       <  1:924191/100‑1 (MQ=255)
                                            cACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGa                                                      >  1:307768/1‑100 (MQ=255)
                                             aCGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGACACTCGACTATGAAACGCTGAt                                                     >  2:4949047/1‑100 (MQ=255)
                                             aCGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGAt                                                     <  2:3674965/100‑1 (MQ=255)
                                              cGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATg                                                    <  1:1705938/100‑1 (MQ=255)
                                              cGAAGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATg                                                    >  2:4660406/1‑100 (MQ=255)
                                               gAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATgg                                                   >  2:4899301/1‑100 (MQ=255)
                                                aGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATggc                                                  >  1:347225/1‑100 (MQ=255)
                                                  gggCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATggcgg                                                <  2:1069785/100‑1 (MQ=255)
                                                  gggCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATggcgg                                                <  1:2650687/100‑1 (MQ=255)
                                                   ggCGAATACGTTCAAGCGGTGAACACTAGAACACGCACCCAGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATggcggc                                               <  2:2900513/100‑1 (MQ=255)
                                                   ggCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATggcggc                                               >  1:977698/1‑100 (MQ=255)
                                                   ggCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATggcggc                                               >  1:4044057/1‑100 (MQ=255)
                                                    gCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGgcggcg                                              <  1:855597/100‑1 (MQ=255)
                                                     cAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGgcgc                                             >  2:4485215/1‑100 (MQ=255)
                                                     cAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGgcgc                                             >  1:1462858/1‑100 (MQ=255)
                                                      aaaCACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCAACTATGAAACGCTGATGGCGGCgct                                            <  1:2096854/100‑1 (MQ=255)
                                                       aaCACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCgctg                                           >  1:2195123/1‑100 (MQ=255)
                                                        acacGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCgctgc                                          <  1:4954379/100‑1 (MQ=255)
                                                         cacGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCggtgct                                         <  2:4297666/100‑1 (MQ=255)
                                                          acGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGctgctg                                        <  1:4846263/100‑1 (MQ=255)
                                                          acGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGctgctg                                        >  2:3972313/1‑100 (MQ=255)
                                                           cGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCtgctgc                                       >  1:3235196/1‑100 (MQ=255)
                                                            gTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCa                                      <  1:226026/100‑1 (MQ=255)
                                                             ttCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCAt                                     >  1:4770136/1‑100 (MQ=255)
                                                               cAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGa                                   >  1:5125015/1‑100 (MQ=255)
                                                               cAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGa                                   >  2:756014/1‑100 (MQ=255)
                                                                aaGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGAc                                  <  1:831365/100‑1 (MQ=255)
                                                                aaGCGGTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGAc                                  >  2:3049947/1‑100 (MQ=255)
                                                                  gCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGt                                >  2:523372/1‑100 (MQ=255)
                                                                  gCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGt                                <  1:3105901/100‑1 (MQ=255)
                                                                   cGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTg                               <  1:3386430/100‑1 (MQ=255)
                                                                    ggTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGa                              >  2:2424874/1‑100 (MQ=255)
                                                                    ggTGAACCCTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGa                              >  1:2908380/1‑100 (MQ=255)
                                                                     gTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGCCTATGAAACGCTGATGGCGGCGCTGCTGCATGACATGAt                             <  2:3904144/100‑1 (MQ=255)
                                                                        aaCCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGa                          <  2:3445710/100‑1 (MQ=255)
                                                                        aaCCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGa                          <  1:4145294/100‑1 (MQ=255)
                                                                         aCCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGaa                         <  1:1128960/100‑1 (MQ=255)
                                                                         aCCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGaa                         >  2:1156846/1‑100 (MQ=255)
                                                                           ccTATATCACGCACCCGGTAGCGGTTGCATGAATTCTGACCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGa                       <  2:17640/100‑1 (MQ=255)
                                                                            cTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGAt                      >  2:3772231/1‑100 (MQ=255)
                                                                            cTATATCACGCACCCGGTAGCGGCTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGAt                      >  2:2661398/1‑100 (MQ=255)
                                                                               tatCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACt                   >  1:2365507/1‑100 (MQ=255)
                                                                               tatCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACt                   >  1:2531063/1‑100 (MQ=255)
                                                                                atCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTc                  >  2:2350225/1‑100 (MQ=255)
                                                                                 tCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAACCTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTcc                 <  2:1651412/100‑1 (MQ=255)
                                                                                  cACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTccc                >  2:836228/1‑100 (MQ=255)
                                                                                  cACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTccc                >  1:3926559/1‑100 (MQ=255)
                                                                                   aCGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGTTGCATGACGTGATTGAAGATACTCCCg               >  2:2505712/1‑100 (MQ=255)
                                                                                   aCGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCg               >  1:1531581/1‑100 (MQ=255)
                                                                                   aCGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCg               <  1:3560725/100‑1 (MQ=255)
                                                                                    cGCACCCGGTAGGGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGc              >  1:4433422/1‑100 (MQ=255)
                                                                                     gCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGcc             <  2:1046393/100‑1 (MQ=255)
                                                                                     gCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGcc             <  1:899856/100‑1 (MQ=255)
                                                                                      cACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCa            <  2:3158895/100‑1 (MQ=255)
                                                                                      cACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCa            >  1:2173299/1‑100 (MQ=255)
                                                                                      cACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCa            >  2:960888/1‑100 (MQ=255)
                                                                                       aCCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCAc           <  2:5045743/100‑1 (MQ=255)
                                                                                        cccGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCAcc          <  1:653113/100‑1 (MQ=255)
                                                                                        cccGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCAcc          <  1:2303563/100‑1 (MQ=255)
                                                                                            gTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTAcc      >  2:4405529/1‑100 (MQ=255)
                                                                                            gTAGCGGTTGCCTGCACTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGCCGTGATTGAAGATACTCCCGCCACCTAcc      <  2:589958/100‑1 (MQ=255)
                                                                                                   |                                                                                                
AAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGA  >  NC_000913/3822455‑3822650

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