Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,608,462 G→C 6.1% D36E (GAC→GAG) ‡ rsuA ← 16S rRNA pseudouridine(516) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,608,4620GC 6.1% 210.9 / 2.5 65D36E (GAC→GAG) ‡rsuA16S rRNA pseudouridine(516) synthase
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base C (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GGTCGAGAACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGTTCT  >  NC_002947/1608314‑1608608
                                                                                                                                                    |                                                                                                                                                  
ggTCGAGAACGGTATGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCggg                                                                                                                                                    >  4:220107/1‑149 (MQ=255)
   cGAGAACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCg                                                                                                                                                 >  1:218574/1‑149 (MQ=255)
    gagaACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGc                                                                                                                                                >  2:395453/1‑149 (MQ=255)
        acggtacggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGACCCGCTGAc                                                                                                                                            <  1:8734/149‑1 (MQ=255)
              cggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCt                                                                                                                                       <  1:334496/148‑1 (MQ=255)
              cggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTg                                                                                                                                      <  4:136828/149‑1 (MQ=255)
              cggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTg                                                                                                                                      <  4:277900/149‑1 (MQ=255)
               ggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTgg                                                                                                                                     >  6:260332/1‑149 (MQ=255)
                   ttGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCAt                                                                                                                                 >  8:135618/1‑149 (MQ=255)
                     gCGGATCGTGGGTGGCGCTTACGCAGCCAGAGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATcc                                                                                                                               >  5:383128/1‑149 (MQ=255)
                       ggATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCAc                                                                                                                             >  1:105368/1‑149 (MQ=255)
                              gtgtggCGCTTACGCAGCCAGTTGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCAcg                                                                                                                      <  6:383128/147‑1 (MQ=255)
                                gtggCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCAcgcg                                                                                                                    >  6:235642/1‑149 (MQ=255)
                                 tggCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACgcgg                                                                                                                   <  5:138302/149‑1 (MQ=255)
                                     gcTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGtt                                                                                                               >  2:250423/1‑149 (MQ=255)
                                         aCGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTgcg                                                                                                            <  2:218574/148‑1 (MQ=255)
                                          cGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTgcgc                                                                                                           <  3:220107/148‑1 (MQ=255)
                                           gCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTgcgc                                                                                                           >  8:51966/1‑147 (MQ=255)
                                             aGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGCCTCCCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGc                                                                                                       >  4:303188/1‑149 (MQ=255)
                                               ccAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCa                                                                                                      >  5:180499/1‑148 (MQ=255)
                                                     ggCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGc                                                                                                <  8:234598/148‑1 (MQ=255)
                                                       cTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGaa                                                                                             <  7:135618/149‑1 (MQ=255)
                                                          gTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAAc                                                                                            <  4:256952/147‑1 (MQ=255)
                                                          gTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGACCCGCTGAccg                                                                                                                                          >  4:192202/1‑99 (MQ=255)
                                                          gTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGACCCGCTGAccg                                                                                                                                          <  3:192202/101‑3 (MQ=255)
                                                           tACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGtt                                                                                         >  4:238004/1‑149 (MQ=255)
                                                                    aGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGtt                                                                                <  4:348460/149‑1 (MQ=255)
                                                                     ggTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTg                                                                               <  3:39818/149‑1 (MQ=255)
                                                                           gggccgggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTg                                                                         >  3:336578/1‑149 (MQ=255)
                                                                           gggccgggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTg                                                                         >  7:128513/1‑149 (MQ=255)
                                                                           gggccgggcGGCCGGCCTGCAGAAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTg                                                                         >  7:248052/1‑149 (MQ=255)
                                                                            ggccgggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTgg                                                                        <  5:235642/149‑1 (MQ=255)
                                                                             gccgggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTggg                                                                       >  8:217197/1‑149 (MQ=255)
                                                                                gggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAg                                                                    <  8:208376/149‑1 (MQ=255)
                                                                                        ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  2:105368/149‑1 (MQ=255)
                                                                                        ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  1:250423/149‑1 (MQ=255)
                                                                                        ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  5:260332/149‑1 (MQ=255)
                                                                                        ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  2:349125/149‑1 (MQ=255)
                                                                                        ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  8:128513/149‑1 (MQ=255)
                                                                                            tGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGa                                                         >  8:27660/1‑148 (MQ=255)
                                                                                                   gTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTc                                                  >  7:359998/1‑148 (MQ=255)
                                                                                                   gTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTc                                                  >  3:145680/1‑148 (MQ=255)
                                                                                                    ttGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCaa                                                <  8:248052/149‑1 (MQ=255)
                                                                                                    ttGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCaa                                                <  6:180499/149‑1 (MQ=255)
                                                                                                      gCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAgg                                              >  1:247996/1‑149 (MQ=255)
                                                                                                      gCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAgg                                              >  3:228630/1‑149 (MQ=255)
                                                                                                      gCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAg                                               <  5:149682/148‑1 (MQ=255)
                                                                                                          gTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCa                                          >  1:108499/1‑149 (MQ=255)
                                                                                                            ccAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  4:174792/129‑1 (MQ=255)
                                                                                                            ccAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            >  3:174792/1‑129 (MQ=255)
                                                                                                              aGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATc                                                                                                                                <  3:318789/59‑1 (MQ=255)
                                                                                                              aGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATc                                                                                                                                >  4:318789/1‑59 (MQ=255)
                                                                                                                ctctACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCaa                                    >  8:351875/1‑149 (MQ=255)
                                                                                                                     cgcgGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACg                                  <  2:247996/146‑1 (MQ=255)
                                                                                                                            gAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTggg                         >  5:344912/1‑148 (MQ=255)
                                                                                                                                   cgcacgTCGGTCAGCGGGTCGCTGACTACCTGTCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAACGCGCATGGCAACGATCAGTGGGCAGCGGAt                 <  4:145680/149‑1 (MQ=255)
                                                                                                                                    gcacgTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGATCTCGCGGTAGTAGCTGCGCAGGTTCGCGAGGAAACGGTCAAGCCGCAAGCCAACGAACAGTGGGCAGCGCATg                >  2:203561/1‑149 (MQ=255)
                                                                                                                                        gTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTgg                          <  1:176635/135‑1 (MQ=255)
                                                                                                                                        gTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTgg                          >  2:176635/1‑135 (MQ=255)
                                                                                                                                         gtagTCAGCGACCCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACg            >  1:1095/4‑148 (MQ=255)
                                                                                                                                           ggTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCCCGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAg         <  4:228630/149‑1 (MQ=255)
                                                                                                                                           ggTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAg         >  3:138629/1‑149 (MQ=255)
                                                                                                                                            gTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCCGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGc        >  4:375970/1‑149 (MQ=255)
                                                                                                                                                 cGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttc   >  8:148571/1‑149 (MQ=255)
                                                                                                                                                 cGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttc   <  7:217197/149‑1 (MQ=255)
                                                                                                                                                 cGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttc   >  7:200506/1‑149 (MQ=255)
                                                                                                                                                  gggTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttct  <  7:27660/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTCGAGAACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGTTCT  >  NC_002947/1608314‑1608608

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