Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,006,479 C→A 14.8% D1020E (GAC→GAA)  PP_2627 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,006,4790CA14.8% 159.1 / 18.5 61D1020E (GAC→GAA) PP_2627membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/20);  new base A (6/3);  total (38/23)
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.17e-01

AGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAGCCCTTG  >  NC_002947/3006333‑3006626
                                                                                                                                                  |                                                                                                                                                    
aGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACaa                                                                                                                                                    >  6:316309/1‑149 (MQ=255)
         cAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGc                                                                                                                                           >  4:129413/1‑149 (MQ=255)
            cAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCtt                                                                                                                                        <  5:316309/149‑1 (MQ=255)
               cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGcc                                                                                                                                          >  6:24378/1‑144 (MQ=255)
               cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGcc                                                                                                                                          <  5:24378/144‑1 (MQ=255)
                                 cccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGc                                                                                                                    >  6:18186/1‑148 (MQ=255)
                                   cTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAAt                                                                                                                 >  3:68047/1‑149 (MQ=255)
                                           tgAAGGCAACCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑tt                                                                                                        <  3:283133/149‑1 (MQ=255)
                                            gAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTc                                                                                                       >  2:52741/1‑149 (MQ=255)
                                               ggCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTg                                                                                                    <  5:18186/149‑1 (MQ=255)
                                                 cACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGcc                                                                                                  >  4:59940/1‑149 (MQ=255)
                                                  aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGcc                                                                                                  >  4:59226/1‑148 (MQ=255)
                                                  aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGcc                                                                                                  >  1:306026/1‑148 (MQ=255)
                                                     cGATGCTTCTTTGGATGACTTCGAGGCCTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGAACGGGCAATCCATGA‑TTCGTGCCGACg                                                                                              >  7:160145/1‑149 (MQ=255)
                                                         gCTTATTTGGATGACTTCGTGGCGTTCTTCGGGCCTAATGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGa                                                                                           <  2:306026/148‑1 (MQ=255)
                                                              tttGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGc                                                                                                   >  8:1223/1‑135 (MQ=255)
                                                              tttGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGc                                                                                                   <  7:1223/135‑1 (MQ=255)
                                                                tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCt                                                                                                             <  8:214450/124‑1 (MQ=255)
                                                                tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCt                                                                                                             >  7:214450/1‑124 (MQ=255)
                                                                tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGagc                                                                                     <  2:144181/147‑1 (MQ=255)
                                                                  gATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCt                                                                                 >  5:145723/1‑149 (MQ=255)
                                                                  gATGACTTCGAGGCGTTCTTCCGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCt                                                                                 >  8:101279/1‑149 (MQ=255)
                                                                   aTGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTc                                                                                >  5:239769/1‑149 (MQ=255)
                                                                    tGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTc                                                                                >  3:124672/1‑148 (MQ=255)
                                                                        ttCGATGCGTTCTTCGGGCCGAAAGGTCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTGCTGAAGTACAACCTTGATGCCTTGCAGTCCGCTCAGCACTGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCt                                                                           <  8:160145/149‑1 (MQ=255)
                                                                        ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCt                                                                           >  3:299500/1‑149 (MQ=255)
                                                                            aGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCg                                                                       <  4:124672/149‑1 (MQ=255)
                                                                               cGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGa                                                                     <  1:52741/149‑1 (MQ=255)
                                                                                gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                                   >  2:104690/1‑149 (MQ=255)
                                                                                    ttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGcgcggag                                                               >  6:206904/1‑146 (MQ=255)
                                                                                    ttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCAt                                                               >  2:177565/1‑149 (MQ=255)
                                                                                       gggCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATcc                                                             <  6:239769/148‑1 (MQ=255)
                                                                                                  gcgcCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACtttct                                                 <  4:299500/149‑1 (MQ=255)
                                                                                                            cagTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCg                                       >  4:243222/1‑149 (MQ=255)
                                                                                                                  aaTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                                   >  3:199796/1‑115 (MQ=255)
                                                                                                                  aaTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                                   <  4:199796/115‑1 (MQ=255)
                                                                                                                  aaTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGc                                   >  2:193093/1‑147 (MQ=255)
                                                                                                                   aTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAg                                <  5:206904/149‑1 (MQ=255)
                                                                                                                    tGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCg                                                                       >  2:40861/1‑109 (MQ=255)
                                                                                                                    tGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCg                                                                       <  1:40861/109‑1 (MQ=255)
                                                                                                                      aTCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCt                              <  6:110140/148‑1 (MQ=255)
                                                                                                                              ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCg                                                                               <  5:174878/91‑1 (MQ=38)
                                                                                                                              ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCg                                                                               >  6:174878/1‑91 (MQ=38)
                                                                                                                              ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg                       >  2:68595/1‑147 (MQ=255)
                                                                                                                               tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg                    >  8:110628/1‑149 (MQ=255)
                                                                                                                               tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg                    >  1:307809/1‑149 (MQ=255)
                                                                                                                               tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttctt                                                >  5:218507/1‑121 (MQ=255)
                                                                                                                               tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttctt                                                <  6:218507/121‑1 (MQ=255)
                                                                                                                               tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg                    >  3:114634/1‑149 (MQ=255)
                                                                                                                               tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGGGATCGAGCAGCTCGCGCTGGCGGAGCGCATCAGCGAAACCTTCTGCGCTCAGCGTGGCAAGCTAAGTGTGCAg                    >  8:283981/1‑149 (MQ=255)
                                                                                                                                cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt                   >  1:191833/1‑149 (MQ=255)
                                                                                                                                cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt                   >  1:178306/1‑149 (MQ=255)
                                                                                                                                 aaGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt                   >  7:96686/1‑148 (MQ=255)
                                                                                                                                   ggTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCa                <  4:114634/149‑1 (MQ=255)
                                                                                                                                        tCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATc           >  3:290315/1‑149 (MQ=255)
                                                                                                                                         ccTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGCAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATc           >  8:209220/1‑148 (MQ=255)
                                                                                                                                           tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCa                     <  3:253525/136‑1 (MQ=255)
                                                                                                                                           tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGCCAAGCTAAGTGTGCa                     >  4:253525/1‑136 (MQ=255)
                                                                                                                                             aaGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCACTGA‑TTAGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAAACTTCTTCGATCAGCGTCGCAAGCTAAGTGTGCAGTACAGCATCGAGcc      >  2:115025/1‑149 (MQ=255)
                                                                                                                                              aGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAGccc     <  7:101279/149‑1 (MQ=255)
                                                                                                                                                 aCAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAGCCCTTg  <  2:191833/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAGCCCTTG  >  NC_002947/3006333‑3006626

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: ‑