Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,554 C→A 10.4% R49L (CGC→CTC)  PP_1191 ← S4 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,366,5540CA10.4% 164.1 / 9.5 58R49L (CGC→CTC) PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (22/30);  new base A (3/3);  total (25/33)
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.32e-02

GGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCG  >  NC_002947/1366406‑1366694
                                                                                                                                                     |                                                                                                                                            
ggCGCGCACGCCTTCGCGGTGGGCCTCTCTGAGGATGCCCATGTTCATGCTGGCGCGGGCATCTTCAATGTCAATGACGGCTGTCTGTTTCAGCAGCAGTGTCACCGGATACAGGGTTTCGGCACCGGCGCCCGGCAGCAGCTCG‑Acgc                                                                                                                                              >  2:257443/1‑149 (MQ=255)
ggCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGCGTTTCGGCCCGGGCGCCCGGCAGCCGCTCG‑Acgc                                                                                                                                              >  6:245641/1‑149 (MQ=255)
    cgcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCgctg                                                                                                                                          <  5:181923/149‑1 (MQ=255)
    cgcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCgctg                                                                                                                                          <  1:90919/149‑1 (MQ=255)
                cggcggGCCGCGCTGAGGCGGCCCATCTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCtt                                                                                                                              <  8:203246/149‑1 (MQ=255)
                cgGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCtt                                                                                                                              <  8:194034/149‑1 (MQ=255)
                 ggTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTg                                                                                                                             <  6:108724/149‑1 (MQ=255)
                   tGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGaa                                                                                                                           <  7:26375/149‑1 (MQ=255)
                   tGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGaa                                                                                                                           <  2:164434/149‑1 (MQ=255)
                       ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGc                                                                                                                       >  6:329897/1‑149 (MQ=255)
                       ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGc                                                                                                                       >  8:235529/1‑149 (MQ=255)
                       ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGc                                                                                                                       >  4:264290/1‑149 (MQ=255)
                         tCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCgg                                                                                                                     >  3:228276/1‑149 (MQ=255)
                           gCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCgg                                                                                                                     <  4:8284/147‑1 (MQ=255)
                             tGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGgctgc                                                                                                                 >  2:208960/1‑149 (MQ=255)
                                ggCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGAACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTc                                                                                                               <  1:208960/149‑1 (MQ=255)
                                   ttcccATTTTCATGCTGTCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGacc                                                                                                           <  7:198883/147‑1 (MQ=255)
                                              aTGCTGGCGCGGGCTTCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATc                                                                                                <  2:105194/149‑1 (MQ=255)
                                                    gcgcGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCATGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCg                                                                                           <  6:162591/148‑1 (MQ=255)
                                                    gcgcGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGt                                                                                          <  3:264290/149‑1 (MQ=255)
                                                      gcgGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTAG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTa                                                                                        >  4:81137/1‑149 (MQ=255)
                                                                 cACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGcctc                                                                               <  7:235529/147‑1 (MQ=255)
                                                                         aTGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTa                                                                     >  1:170897/1‑149 (MQ=255)
                                                                         aTGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTa                                                                     >  2:135052/1‑149 (MQ=255)
                                                                          tGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTAc                                                                    <  4:228276/149‑1 (MQ=255)
                                                                               gctggctgGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTc                                                               >  8:292398/1‑149 (MQ=255)
                                                                                 tggctggTTCAGCAGAAGGGTCACCGGCTCCAGGGTTTCGGCACTGGCTCCAGGCAGCAGCTCT‑AAGCGCTGATGATCAACCTTGAACTGTGTCTGCTCGACCACAACGCCATCTACCGTTACCCAGCCTCCTTCAATGTACTGCTctg                                                             <  1:257443/149‑3 (MQ=255)
                                                                                  ggctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGCGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCCCCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGc                                                            >  8:6872/1‑149 (MQ=255)
                                                                                   gctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGcc                                                           >  6:273383/1‑149 (MQ=255)
                                                                                   gctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGc                                                            >  8:143035/1‑148 (MQ=255)
                                                                                       gTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTcgc                                                       <  2:106733/149‑1 (MQ=255)
                                                                                           agcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggcg                                                   <  5:245641/149‑1 (MQ=255)
                                                                                            gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                  <  2:170897/149‑1 (MQ=255)
                                                                                            gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACACAGCCTCCTTCACTGTACTGCTCGGCCTCGcggcgg                                                  >  8:70475/1‑149 (MQ=255)
                                                                                            gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGACGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACACCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                  <  7:298235/149‑1 (MQ=255)
                                                                                              agcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                <  5:329897/149‑1 (MQ=255)
                                                                                                   ggTCACCGGCTCCAGGGTTTAGGCACGGGCCCACGGCAGCAGCTAG‑ACGCGCTGGTGCTCAACCATGAACTGCGGCTGCTCGACCCCCACGCCATCGACCGTTAACCAGCCTCCTTCAATTTACTGCTCGGACTCGCGGCGCGagcacc                                           >  2:194465/1‑147 (MQ=255)
                                                                                                       aCCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTCAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAg                                       <  7:6872/149‑1 (MQ=255)
                                                                                                         cGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAgct                                     >  1:94965/1‑149 (MQ=255)
                                                                                                               cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTAGCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTc                                 <  7:292398/147‑1 (MQ=255)
                                                                                                               cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                <  4:130867/132‑1 (MQ=255)
                                                                                                               cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                >  3:130867/1‑132 (MQ=255)
                                                                                                               cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCTCGGCTGCTcc                                <  8:20115/148‑2 (MQ=255)
                                                                                                               cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGa                               >  1:193632/1‑149 (MQ=255)
                                                                                                                aGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGGCCTGCGCGAt                              >  3:1369/1‑149 (MQ=255)
                                                                                                                   gTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCTCGGCTGCTcccgccg                           <  3:81137/149‑7 (MQ=255)
                                                                                                                           aCGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTg                    >  8:92609/1‑148 (MQ=255)
                                                                                                                                 gcCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGc             <  5:273383/149‑1 (MQ=255)
                                                                                                                                  cccGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCg            <  8:273677/149‑1 (MQ=255)
                                                                                                                                    cGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGa          <  2:193632/149‑1 (MQ=255)
                                                                                                                                     ggCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGAc         >  5:123555/1‑149 (MQ=255)
                                                                                                                                      gcagcaGCGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAAt                                                                        >  7:36448/1‑85 (MQ=39)
                                                                                                                                      gcagcaGCGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAAt                                                                        <  8:36448/85‑1 (MQ=39)
                                                                                                                                      gcagcaGCGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACg        >  2:287336/1‑149 (MQ=255)
                                                                                                                                        agcagcGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACggg      >  5:150433/1‑149 (MQ=255)
                                                                                                                                         gcagcTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGt     >  1:233835/1‑149 (MQ=255)
                                                                                                                                         gcagcGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGt     <  1:287336/149‑1 (MQ=255)
                                                                                                                                           agcGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTc   <  6:150433/149‑1 (MQ=255)
                                                                                                                                            gcTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCg  >  5:233314/1‑149 (MQ=255)
                                                                                                                                            gcTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCg  <  1:135052/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCG  >  NC_002947/1366406‑1366694

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