Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,944 G→T 5.1% noncoding (123/140 nt) ryjA ← novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,277,9440GT 5.1% 959.1 / 26.8 257noncoding (123/140 nt)ryjAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (111/133);  new base T (6/7);  total (117/140)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 264 total aligned reads displayed.

GTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTT  >  NC_000913/4277796‑4278092
                                                                                                                                                    |                                                                                                                                                    
gTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAgg                                                                                                                                                     <  2:4251981/150‑1 (MQ=255)
 tGGCGCTGCGTGACGAACTGGACGGATGTTTTGGTTGTGGCTGCCTTTCGCGCTGTGATTGCCCTTTGCGTAACCCGGGCGACCGCTTAGTAGAATAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAACGCGCCACAAggg                                                                                                                                                    <  1:3255186/150‑1 (MQ=255)
 tGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAggg                                                                                                                                                    >  2:4089454/1‑150 (MQ=255)
   gcgcTGCGTGACTAAGTGGACGGTTGTTTTGGTTGGGGCTGCCTTTCGCGCTGTGATTGCCCGTTGCGTAACCGGGGCGACCGCTTCGGATAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgcg                                                                                                                                                  <  1:655301/150‑1 (MQ=255)
   gcgcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgcg                                                                                                                                                  >  2:1506477/1‑150 (MQ=255)
   gcgcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgcg                                                                                                                                                  >  2:1290119/1‑150 (MQ=255)
    cgcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgcgc                                                                                                                                                 >  2:2298929/1‑150 (MQ=255)
     gcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt                                                                                                                                                >  2:3844980/1‑150 (MQ=255)
     gcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt                                                                                                                                                >  2:288109/1‑150 (MQ=255)
      cTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCtt                                                                                                                                               >  1:716871/1‑150 (MQ=255)
        gCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTa                                                                                                                                             <  2:4139137/150‑1 (MQ=255)
         cGTGACGAACTGGACGGTTGTTTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAg                                                                                                                                            <  1:853392/150‑1 (MQ=255)
         cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATGGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAg                                                                                                                                            <  1:2752818/150‑1 (MQ=255)
          gTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgt                                                                                                                                           <  2:1028489/150‑1 (MQ=255)
          gTGACGAACTGGACGGAGGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgt                                                                                                                                           >  2:4181616/1‑150 (MQ=255)
            gACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAATAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgttt                                                                                                                                         <  1:305202/150‑1 (MQ=255)
             aCGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttg                                                                                                                                        >  1:1944928/1‑150 (MQ=255)
              ctaaCTGGACGGTTTTTTTGGTTGTGGCGGCCTTTCGCGCTGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAATATCAACAAAACTAACGCGCCCTTGTGGCGCTTTAgttttg                                                                                                                                       <  1:1487177/148‑3 (MQ=255)
               gAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgttttgt                                                                                                                                      <  1:596159/150‑4 (MQ=255)
               gAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttgtt                                                                                                                                      >  1:2688541/1‑150 (MQ=255)
               gAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttgtt                                                                                                                                      >  2:2680529/1‑150 (MQ=255)
                 aCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGtttt                                                                                                                                    >  1:453263/1‑150 (MQ=255)
                           tGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttg                                                                                                                          <  2:1686511/150‑1 (MQ=255)
                            gTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgt                                                                                                                         >  1:1008392/1‑150 (MQ=255)
                              aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtct                                                                                                                       <  2:2900171/150‑1 (MQ=255)
                              aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtct                                                                                                                       >  2:3016510/1‑150 (MQ=255)
                               ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtctt                                                                                                                      <  2:3370778/150‑1 (MQ=255)
                                tGGTTGTGGCTGCCTTTCGCGCAGTGAGTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGATGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt                                                                                                                     <  1:4062525/150‑1 (MQ=255)
                                  gTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtct                                                                                                                   >  1:2967669/1‑150 (MQ=255)
                                   ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctc                                                                                                                  <  2:336802/150‑1 (MQ=255)
                                   ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctc                                                                                                                  >  2:932188/1‑150 (MQ=255)
                                    tgtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctct                                                                                                                 >  2:3661158/1‑150 (MQ=255)
                                      tgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTctctct                                                                                                               <  2:1439190/150‑1 (MQ=255)
                                       ggCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTa                                                                                                              <  2:3889143/150‑1 (MQ=255)
                                        gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTAt                                                                                                             >  2:44244/1‑150 (MQ=255)
                                          tGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATcc                                                                                                           <  1:538234/150‑1 (MQ=255)
                                              tttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                                                                                       <  1:3843870/150‑1 (MQ=255)
                                              tttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                                                                                       <  2:332134/150‑1 (MQ=255)
                                              tttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                                                                                       <  2:3898895/150‑1 (MQ=255)
                                               ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTg                                                                                                      <  1:1195309/150‑1 (MQ=255)
                                                  gcgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTa                                                                                                   >  1:19149/1‑150 (MQ=255)
                                                      aGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTAcaca                                                                                               <  2:2061514/150‑1 (MQ=255)
                                                            tGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAggg                                                                                         >  1:1364179/1‑150 (MQ=255)
                                                            tGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAggg                                                                                         >  1:3112969/1‑150 (MQ=255)
                                                               ccGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGttt                                                                                      <  1:392456/150‑1 (MQ=255)
                                                                      gTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGa                                                                               <  2:4041559/150‑1 (MQ=255)
                                                                        aaCCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTGGCGTGGGTTTTCCGGTCTTTGGCTTTCTCGCTCGCCCGCTGGTCCACAGGAGGGTTTCCCCCGACg                                                                             >  2:3083636/1‑150 (MQ=255)
                                                                            cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                         <  2:3206187/150‑1 (MQ=255)
                                                                            cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                         <  1:2035727/150‑1 (MQ=255)
                                                                            cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                         <  2:2552883/150‑1 (MQ=255)
                                                                             gggCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAac                                                                        >  1:1774602/1‑150 (MQ=255)
                                                                             gggCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAac                                                                        >  2:3776101/1‑150 (MQ=255)
                                                                             gggCGACCGCATAGGAGAAGAAGGAACCGGAGAACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGTCGCTTTAGTTTGTTTTCCGGTCTTTGACTTTCTCTCAAACCCGCATGTACACACGAGGGATTCCCCCGACGACAac                                                                        >  2:849641/1‑150 (MQ=255)
                                                                                cGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacaca                                                                     >  1:400044/1‑150 (MQ=255)
                                                                                 gACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacacac                                                                    <  1:4135504/150‑1 (MQ=255)
                                                                                 gACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacacac                                                                    <  2:194471/150‑1 (MQ=255)
                                                                                   ccGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCt                                                                  >  1:1082951/1‑150 (MQ=255)
                                                                                    cGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTc                                                                 >  2:681902/1‑150 (MQ=255)
                                                                                    cGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCGATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTc                                                                 <  1:1882988/150‑1 (MQ=255)
                                                                                      cTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCAt                                                               >  1:3291712/1‑150 (MQ=255)
                                                                                       ttAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCAtt                                                              >  2:2431693/1‑150 (MQ=255)
                                                                                       ttAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTGTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCAtt                                                              >  2:2431714/1‑150 (MQ=255)
                                                                                          ggAGAAGAAGGTACCGGCGCACGCTTTCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGa                                                           <  1:1591350/150‑1 (MQ=255)
                                                                                          ggAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGa                                                           >  1:2513137/1‑150 (MQ=255)
                                                                                           gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                          <  1:279525/150‑1 (MQ=255)
                                                                                           gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                          <  1:1290119/150‑1 (MQ=255)
                                                                                              aagaagGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTCGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCAc                                                       >  2:275287/1‑150 (MQ=255)
                                                                                                gaagGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgt                                                     >  2:4294393/1‑150 (MQ=255)
                                                                                                gaagGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgt                                                     >  2:2150837/1‑150 (MQ=255)
                                                                                                gaagGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgt                                                     >  1:1229453/1‑150 (MQ=255)
                                                                                                   ggTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtggt                                                  <  2:104932/150‑1 (MQ=255)
                                                                                                   ggTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtggt                                                  <  2:3733352/150‑1 (MQ=255)
                                                                                                    gTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtggtg                                                 <  1:4181616/150‑1 (MQ=255)
                                                                                                      aCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGa                                               <  1:315590/150‑1 (MQ=255)
                                                                                                      aCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGa                                               <  2:1871223/150‑1 (MQ=255)
                                                                                                       ccGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                              <  2:1097932/150‑1 (MQ=255)
                                                                                                             cACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCg                                        >  1:4249982/1‑150 (MQ=255)
                                                                                                                 cTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTg                                    <  1:3016510/150‑1 (MQ=255)
                                                                                                                   ttcTGGAAGATGAACTAAACTAAAGCGCCTCAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgt                                  <  2:2446275/148‑1 (MQ=255)
                                                                                                                   tGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgt                                  <  1:3827973/150‑1 (MQ=255)
                                                                                                                     cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGATTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                <  2:2363649/150‑1 (MQ=255)
                                                                                                                      tGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGGTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATGCGGCTGGTAGACTGGTGGGTTTCCCCCGTCGTCAAGACACCTCATTCGAGCACGTGGTGGAGGTTGGGGTTGGTGtt                               >  2:816251/1‑150 (MQ=255)
                                                                                                                       ggAAGATGAACAAAACTAAAGCGCCACACGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                              >  2:3976222/1‑150 (MQ=255)
                                                                                                                        gAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGa                             <  2:1163923/150‑1 (MQ=255)
                                                                                                                                aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                     >  2:472340/1‑150 (MQ=255)
                                                                                                                                 cAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAt                    <  2:3407611/150‑1 (MQ=255)
                                                                                                                                  aaaaCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                   >  1:3670980/1‑150 (MQ=255)
                                                                                                                                  aaaaCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                   >  1:3450169/1‑150 (MQ=255)
                                                                                                                                     aCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa                >  2:2245055/1‑150 (MQ=255)
                                                                                                                                     aCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATACCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTACGGTTGGGGGTGATGCTTTAATTGTa                >  1:276254/1‑150 (MQ=255)
                                                                                                                                       tAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTTTATg              >  1:4283917/1‑150 (MQ=255)
                                                                                                                                       tAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATg              >  1:2879952/1‑150 (MQ=255)
                                                                                                                                       tAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTACGGTCTTTTTCTTTATATATATAACGCTGGTACACATGAGGGTTTCCCCCGACGTCAACACAACTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTTTTGATTCTTTAATTGTATg              >  1:99/1‑150 (MQ=255)
                                                                                                                                          aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa           <  1:1694265/150‑1 (MQ=255)
                                                                                                                                          aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa           >  1:3384345/1‑150 (MQ=255)
                                                                                                                                          aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa           >  1:2619995/1‑150 (MQ=255)
                                                                                                                                          aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTGTCTCTCGATCCCGCTGGGACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa           <  1:3484436/150‑1 (MQ=255)
                                                                                                                                                 caAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGt    <  2:1121715/150‑1 (MQ=255)
                                                                                                                                                   aGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGttt  <  1:932188/150‑1 (MQ=255)
                                                                                                                                                   aGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGttt  <  1:232312/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTT  >  NC_000913/4277796‑4278092

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