Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,608,463 T→C 6.0% D36G (GAC→GGC) ‡ 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,4630TC 6.0% 191.4 / 2.4 66D36G (GAC→GGC) ‡rsuA16S rRNA pseudouridine(516) synthase
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base C (2/2);  total (34/32)
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.88e-01

CGAGAACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGTTCT  >  NC_002947/1608317‑1608608
                                                                                                                                                  |                                                                                                                                                 
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:277900/149‑1 (MQ=255)
           cggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTg                                                                                                                                      <  4:136828/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:349125/149‑1 (MQ=255)
                                                                                     ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  5:260332/149‑1 (MQ=255)
                                                                                     ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  2:105368/149‑1 (MQ=255)
                                                                                     ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  1:250423/149‑1 (MQ=255)
                                                                                     ggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGa                                                            <  8:128513/149‑1 (MQ=255)
                                                                                         tGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGa                                                         >  8:27660/1‑148 (MQ=255)
                                                                                                gTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTc                                                  >  3:145680/1‑148 (MQ=255)
                                                                                                gTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTc                                                  >  7:359998/1‑148 (MQ=255)
                                                                                                 ttGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCaa                                                <  8:248052/149‑1 (MQ=255)
                                                                                                 ttGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCaa                                                <  6:180499/149‑1 (MQ=255)
                                                                                                   gCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAgg                                              >  3:228630/1‑149 (MQ=255)
                                                                                                   gCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAgg                                              >  1:247996/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                                                                                                                                >  4:318789/1‑59 (MQ=255)
                                                                                                           aGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATc                                                                                                                                <  3:318789/59‑1 (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                          >  2:176635/1‑135 (MQ=255)
                                                                                                                                     gTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTgg                          <  1:176635/135‑1 (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)
                                                                                                                                                  |                                                                                                                                                 
CGAGAACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGTTCT  >  NC_002947/1608317‑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: ‑