Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,763,642 G→C 6.6% D340H (GAC→CAC)  PP_5057 → M23/M37 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,763,6420GC 6.6% 340.4 / 10.7 106D340H (GAC→CAC) PP_5057M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (49/50);  new base C (3/4);  total (52/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 109 total aligned reads displayed.

CCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCA  >  NC_002947/5763496‑5763790
                                                                                                                                                  |                                                                                                                                                    
ccaccCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGa                                                                                                                                                     <  2:252516/148‑1 (MQ=255)
 caccCTCGGGCCGATGGTTTCCAGCGCTGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGGCAGTCTATGGTCGATTGCTGGCACGTGTCGGTGATGCCCGCGGTGGCGACGCGCGTGCCAAGGGGGACg                                                                                                                                                   >  8:330250/1‑149 (MQ=255)
 caccCTCGGGCCGATGGTTTCCAGCCATGGCGCCAAATACGGCGGCGCATTTTCTGATGCGCGTGGAAAACTTCCATGGCAAGTCAATGGTCGATTGCTGCCACGTTTCGGTGATGCCCGCTGTAGCGACGCCCGTGCCAAGTGGCACg                                                                                                                                                   >  2:374021/1‑149 (MQ=255)
   cccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACggg                                                                                                                                                 >  6:314876/1‑149 (MQ=255)
   cccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACggg                                                                                                                                                 <  3:98222/149‑1 (MQ=255)
    ccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACgggg                                                                                                                                                <  4:252197/149‑1 (MQ=255)
    ccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggc                                                                                                                                                <  6:78241/149‑2 (MQ=255)
     cTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGTAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGACCCACTTGGCACgggcg                                                                                                                                               <  5:44345/149‑3 (MQ=255)
      tCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGt                                                                                                                                               >  6:252248/1‑148 (MQ=255)
      tCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGt                                                                                                                                               >  5:353980/1‑148 (MQ=255)
         ggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATg                                                                                                                                           <  2:380998/149‑1 (MQ=255)
         ggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggcgtcgc                                                                                                                                           <  3:139357/149‑7 (MQ=255)
          gCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGa                                                                                                                                          >  2:266412/1‑149 (MQ=255)
            cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGAGGCCCGTGCCAAGTGGGACGGGGTTATGATc                                                                                                                                        <  1:38934/149‑1 (MQ=255)
            cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGAt                                                                                                                                         >  8:225442/1‑148 (MQ=255)
            cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATc                                                                                                                                        <  7:88884/149‑1 (MQ=255)
                 gTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAgcgc                                                                                                                                   >  1:281422/1‑149 (MQ=255)
                  tttCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGTCGGGGTTATGATCAgcgc                                                                                                                                   >  3:97771/1‑148 (MQ=255)
                    tCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCaa                                                                                                                                >  8:373218/1‑149 (MQ=255)
                     ccAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAc                                                                                                                               >  3:27714/1‑149 (MQ=255)
                       aGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAccc                                                                                                                             >  7:218412/1‑149 (MQ=255)
                       aGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAccc                                                                                                                             >  7:212307/1‑149 (MQ=255)
                       aGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAcc                                                                                                                              >  5:169945/1‑148 (MQ=255)
                        gCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAcc                                                                                                                              >  8:378734/1‑147 (MQ=255)
                        gCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCa                                                                                                                            <  5:271487/149‑1 (MQ=255)
                           aTGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCAAgg                                                                                                                          >  3:47922/1‑148 (MQ=255)
                                    aCTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTa                                                                                                                >  5:124792/1‑149 (MQ=255)
                                    aCTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTa                                                                                                                <  8:218412/149‑1 (MQ=255)
                                    aCTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTa                                                                                                                <  8:212307/149‑1 (MQ=255)
                                     cTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCCGGCACCCAGGTAc                                                                                                               <  4:349565/149‑1 (MQ=255)
                                      tACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcg                                                                                                              >  5:57509/1‑149 (MQ=255)
                                        cggcggCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCAGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgt                                                                                                             <  4:47922/148‑1 (MQ=255)
                                         ggcggcGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgc                                                                                                           >  8:119588/1‑149 (MQ=255)
                                           cggcgCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccg                                                                                                         >  8:242165/1‑149 (MQ=255)
                                           cggcgCATTTTCTGCTGCGCGGGGCAAACTTCAATGTCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccg                                                                                                         <  7:302020/149‑1 (MQ=255)
                                            ggcgCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccgt                                                                                                        <  5:252248/149‑1 (MQ=255)
                                            ggcgCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccgt                                                                                                        >  5:124970/1‑149 (MQ=255)
                                               gcATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCa                                                                                                     >  1:57822/1‑149 (MQ=255)
                                               gcATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGCTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCa                                                                                                     <  7:395239/149‑1 (MQ=255)
                                                  ttttCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCAcgg                                                                                                  <  5:314876/149‑1 (MQ=255)
                                                            cgcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCgtgt                                                                                         <  1:266412/148‑1 (MQ=255)
                                                            cgcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                                        <  5:364189/149‑1 (MQ=255)
                                                             gcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTgg                                                                                       <  8:64008/149‑1 (MQ=255)
                                                             gcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTgg                                                                                       <  4:247862/149‑1 (MQ=255)
                                                              cgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTgg                                                                                       >  6:216378/1‑148 (MQ=255)
                                                               ggggCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTa                                                                                     >  7:350841/1‑149 (MQ=255)
                                                                gggCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTAt                                                                                    >  1:413051/1‑149 (MQ=255)
                                                                      aCTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCgg                                                                              <  5:216378/149‑1 (MQ=255)
                                                                        ttCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGAc                                                                            <  2:281422/149‑1 (MQ=255)
                                                                          ccATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGAc                                                                            <  7:242165/147‑1 (MQ=255)
                                                                           cATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTgg                                                                         <  6:124970/149‑1 (MQ=255)
                                                                           cATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTg                                                                          <  6:432657/148‑1 (MQ=255)
                                                                            aTGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCTGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGt                                                                        >  3:360967/1‑149 (MQ=255)
                                                                               gCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCTGTAGCGACGCCCGTGCCATGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGATTGGTTgc                                                                     <  1:374021/149‑1 (MQ=255)
                                                                               gCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                                        <  8:316462/130‑1 (MQ=255)
                                                                               gCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                                        >  7:316462/1‑130 (MQ=255)
                                                                                  aGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGcgcg                                                                  <  7:373218/149‑1 (MQ=255)
                                                                                  aGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGcgcg                                                                  <  2:57822/149‑1 (MQ=255)
                                                                                  aGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGcgcg                                                                  <  2:64118/149‑1 (MQ=255)
                                                                                     cAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTg                                                               >  6:439980/1‑149 (MQ=255)
                                                                                      aaTGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTg                                                               <  2:413051/148‑1 (MQ=255)
                                                                                      aaTGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGc                                                              >  8:404786/1‑149 (MQ=255)
                                                                                         ggTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCgg                                                           >  5:243729/1‑149 (MQ=255)
                                                                                           tCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGc                                                         <  7:119588/149‑1 (MQ=255)
                                                                                             gATTGCTGGCACGTTTCGGTGATGCCCGCGGGAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCtt                                                       >  1:267650/1‑149 (MQ=255)
                                                                                              aTTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                                      >  7:144915/1‑149 (MQ=255)
                                                                                                   tGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                                      <  2:267178/144‑1 (MQ=255)
                                                                                                   tGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                                      >  1:267178/1‑144 (MQ=255)
                                                                                                   tGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGt                                                  >  5:387549/1‑148 (MQ=255)
                                                                                                     gCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCa                                                <  1:57757/148‑1 (MQ=255)
                                                                                                      cACGTTTCGGTGATGCCCGCGGTAGCGACGGCCGTGCCCAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCAtt                                              <  3:314808/149‑1 (MQ=255)
                                                                                                       aCGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGATATCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATtc                                             >  1:363399/1‑149 (MQ=255)
                                                                                                         gTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                                      <  2:105647/138‑1 (MQ=255)
                                                                                                         gTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                                      >  1:105647/1‑138 (MQ=255)
                                                                                                          tttCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATtctc                                           >  4:410995/1‑148 (MQ=255)
                                                                                                            tCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGAc                                        >  3:169257/1‑149 (MQ=255)
                                                                                                                tGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATg                                    <  6:169945/149‑1 (MQ=255)
                                                                                                                 gATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                                   <  6:243729/149‑1 (MQ=255)
                                                                                                                 gATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAAGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                                   >  2:341199/1‑149 (MQ=255)
                                                                                                                  aTGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCATGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                                   <  1:341199/148‑1 (MQ=255)
                                                                                                                    gCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTc                                                 <  8:80000/132‑1 (MQ=255)
                                                                                                                    gCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTc                                                 >  7:80000/1‑132 (MQ=255)
                                                                                                                    gCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCaa                                >  4:275182/1‑149 (MQ=255)
                                                                                                                     cccGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                                   <  1:70136/145‑1 (MQ=255)
                                                                                                                     cccGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                                   >  2:70136/1‑145 (MQ=255)
                                                                                                                       cgcgGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACgg                             >  7:87568/1‑149 (MQ=255)
                                                                                                                        gcgGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGt                            <  7:378734/149‑1 (MQ=255)
                                                                                                                         cgGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTTTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGtt                           <  4:27714/149‑1 (MQ=255)
                                                                                                                           gTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGTGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTAc                         <  7:143321/149‑1 (MQ=255)
                                                                                                                               cGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCt                       >  3:101348/1‑147 (MQ=255)
                                                                                                                               ataaccccGTCCCACTTGGCACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCt                       >  6:393121/6‑147 (MQ=255)
                                                                                                                                taaccccGTCCCACTTGGCACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAg                    >  6:222303/5‑149 (MQ=255)
                                                                                                                                gACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAg                    >  2:318966/1‑149 (MQ=255)
                                                                                                                                 aCGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGc                   >  6:36973/1‑149 (MQ=255)
                                                                                                                                 aCGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGc                   <  6:57509/149‑1 (MQ=255)
                                                                                                                                   gCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCt                 <  6:124792/149‑1 (MQ=255)
                                                                                                                                          ccAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGc          <  4:169257/149‑1 (MQ=255)
                                                                                                                                                  gACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAAc    <  3:74948/147‑1 (MQ=255)
                                                                                                                                                  gACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAAc    >  4:74948/1‑147 (MQ=255)
                                                                                                                                                  gACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCa  <  5:439980/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCA  >  NC_002947/5763496‑5763790

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: