Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,608,463 T→C 5.5% 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 5.5% 220.8 / 2.4 74D36G (GAC→GGC) ‡rsuA16S rRNA pseudouridine(516) synthase
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

GAGAACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGTTCTT  >  NC_002947/1608318‑1608609
                                                                                                                                                 |                                                                                                                                                  
gagaACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGc                                                                                                                                                 >  2:296018/1‑149 (MQ=255)
 agaACGGTACGGGGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGGGGGCTTGTACAGCATCAGGTAGCGGGCCGGTCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCt                                                                                                                                                <  7:249615/149‑1 (MQ=255)
    acggtacggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGACCCGCTGAc                                                                                                                                             <  7:271546/149‑1 (MQ=255)
    acggtacggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCACGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCGGAATTCACGCACGTCGGTCAGCGACCCGCTGAc                                                                                                                                             <  8:188661/149‑1 (MQ=255)
          cggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCt                                                                                                                                        <  3:248746/148‑1 (MQ=255)
           ggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTgg                                                                                                                                      <  3:260950/149‑1 (MQ=255)
           ggtGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTgg                                                                                                                                      >  2:124729/1‑149 (MQ=255)
             tgtTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGc                                                                                                                                     >  7:100245/1‑148 (MQ=255)
                 gCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATc                                                                                                                                 <  3:188995/148‑1 (MQ=255)
                  cGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATcc                                                                                                                                <  7:50390/148‑1 (MQ=255)
                     aTCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTc                                                                                                                            <  2:200761/149‑1 (MQ=255)
                            gtggCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCAcgcg                                                                                                                     >  7:39784/1‑149 (MQ=255)
                             tggCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACgcgg                                                                                                                    >  5:76146/1‑149 (MQ=255)
                                  cTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTg                                                                                                               >  7:51646/1‑149 (MQ=255)
                                   ttACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTgc                                                                                                              >  5:59456/1‑149 (MQ=255)
                                     aCGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAACGGGTc                                                                                                                                                   >  4:69022/1‑110 (MQ=255)
                                     aCGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAACGGGTc                                                                                                                                                   <  3:69022/110‑1 (MQ=255)
                                      cGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGcc                                                                                                           <  1:296018/149‑1 (MQ=255)
                                      cGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGcc                                                                                                           <  3:341984/149‑1 (MQ=255)
                                       gCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACgcggc                                                                                                                   <  2:18312/140‑1 (MQ=255)
                                       gCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACgcggc                                                                                                                   >  1:18312/1‑140 (MQ=255)
                                         aGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAg                                                                                                         >  5:232261/1‑148 (MQ=255)
                                         aGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAg                                                                                                         >  7:180696/1‑148 (MQ=255)
                                         aGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAg                                                                                                         >  6:288896/1‑148 (MQ=255)
                                            cAGTGGTCTTGTCAAGGATCAGTTAGCGGTCCGGGCTGCCTGCCTTCAGGAGTTGCTCGTCCAGCTCTAGGCGGCTGAATTGACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATc                                                                                                     <  8:232650/149‑1 (MQ=255)
                                            cAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCAt                                                                                                      >  3:213732/1‑148 (MQ=255)
                                            cAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGGCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATc                                                                                                     <  2:125794/149‑1 (MQ=255)
                                            cAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGGTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATc                                                                                                     <  3:109271/149‑1 (MQ=255)
                                               tGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAgg                                                                                                  <  6:76146/149‑1 (MQ=255)
                                                 ggCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCg                                                                                                <  4:105474/149‑1 (MQ=255)
                                                  gCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCg                                                                                                >  5:197940/1‑148 (MQ=255)
                                                    ttGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAAc                                                                                             >  4:314106/1‑149 (MQ=255)
                                                      gTACAGCATCAGGTAGCGGTCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACg                                                                                            <  4:126586/148‑1 (MQ=255)
                                                            cATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGc                                                                                     >  1:350918/1‑149 (MQ=255)
                                                                aGGTAGCGGGCCGGGGGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGGCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGtt                                                                                 <  1:136422/149‑1 (MQ=255)
                                                                 tgctgcaggccggccGCCCGGCCCGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGATCTGGCGGTTg                                                                                >  6:152862/8‑149 (MQ=255)
                                                                 ggTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTg                                                                                <  7:256259/149‑1 (MQ=255)
                                                                 ggTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTg                                                                                >  6:327832/1‑149 (MQ=255)
                                                                  gTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGt                                                                               <  8:100245/149‑1 (MQ=255)
                                                                   tAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTa                                                                              <  4:302348/149‑1 (MQ=255)
                                                                   tAGCGGGCCGCGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTa                                                                              <  8:298734/149‑1 (MQ=255)
                                                                          ccgggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCt                                                                           <  7:136213/145‑1 (MQ=255)
                                                                          ccgggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCt                                                                           >  8:136213/1‑145 (MQ=255)
                                                                            gggcGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAg                                                                     <  8:51646/149‑1 (MQ=255)
                                                                                ggccggccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTg                                                                 >  7:177699/1‑149 (MQ=255)
                                                                                     gccTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGGCGGTCAGCGGGTCGCTGACTACCTGGCCATCCCCTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAg                                                            <  5:298363/149‑1 (MQ=255)
                                                                                      ccTGCAGCAGTTGCTCGTCCAGCTCTACGCGACTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAg                                                            >  8:162530/1‑148 (MQ=255)
                                                                                         gcagcaGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAgg                                                           <  4:156732/146‑1 (MQ=255)
                                                                                         gcagcaGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAgg                                                           >  3:156732/1‑146 (MQ=255)
                                                                                         gcagcaGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGaaa                                                        >  5:255815/1‑149 (MQ=255)
                                                                                           agcagTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACg                                                      <  5:288896/149‑1 (MQ=255)
                                                                                              agTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGt                                                    >  4:5445/1‑148 (MQ=255)
                                                                                                      gTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGgcgc                                            <  6:161874/148‑1 (MQ=255)
                                                                                                       tCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCAt                                          <  6:232261/149‑1 (MQ=255)
                                                                                                       tCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCAt                                          >  8:215895/1‑149 (MQ=255)
                                                                                                           gCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGc                                       >  7:308382/1‑148 (MQ=255)
                                                                                                             tctACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAAc                                    >  6:250440/1‑149 (MQ=255)
                                                                                                               tACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCATGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAAc                                    <  8:177699/147‑1 (MQ=255)
                                                                                                                       tGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTggg                          >  1:151801/1‑149 (MQ=255)
                                                                                                                          aTTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAgg                                                                    <  4:93795/104‑1 (MQ=255)
                                                                                                                          aTTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAgg                                                                    >  3:93795/1‑104 (MQ=255)
                                                                                                                             cacgcacgTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCg                     <  6:59456/148‑1 (MQ=255)
                                                                                                                              acgcgcgTCCGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGACCGCATCAGGCGAACGTTCTGGCGGTGGTAGCTGTGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGAGCAGTGGGCAGCTGa                   >  4:197717/1‑149 (MQ=255)
                                                                                                                              acgcacgTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGa                   >  1:85627/1‑149 (MQ=255)
                                                                                                                               cgcacgTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGAt                  <  1:124729/149‑1 (MQ=255)
                                                                                                                                   aggtagTCAGCGACCCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGAc              >  5:67491/6‑149 (MQ=255)
                                                                                                                                       agTCAGCGACCCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAg          >  7:57940/2‑149 (MQ=255)
                                                                                                                                         tCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCt        >  3:330402/1‑149 (MQ=255)
                                                                                                                                          cAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTg       <  6:255815/149‑1 (MQ=255)
                                                                                                                                            gCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGtt     >  2:80156/1‑149 (MQ=255)
                                                                                                                                              gggTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttct   <  4:213732/149‑1 (MQ=255)
                                                                                                                                              gggTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttct   >  3:138679/1‑149 (MQ=255)
                                                                                                                                               ggTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttctt  <  7:162530/149‑1 (MQ=255)
                                                                                                                                                gTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGttctt  >  5:170040/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GAGAACGGTACGGTGTTGCGGATCGTGGGTGGCGCTTACGCAGCCAGTGGGCTTGTACAGCATCAGGTAGCGGGCCGGGCGGCCGGCCTGCAGCAGTTGCTCGTCCAGCTCTACGCGGCTGAATTCACGCACGTCGGTCAGCGGGTCGCTGACTACCTGGCCATCCACTCGCACGCGGCGTTGCGCCAGCATCAGGCGAACGTTCTGGCGGTTGTAGCTGGGCAGGTTGGCGAGGAAACGGTCAAGGCGCATGGCAACGATCAGTGGGCAGCGGATGGACGGAGCTGTTCTT  >  NC_002947/1608318‑1608609

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: