Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,763,638 G→T 8.2% K338N (AAG→AAT) ‡ 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,6380GT 8.2% 185.9 / 6.0 61K338N (AAG→AAT) ‡PP_5057M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base T (2/3);  total (28/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.37e-01

CCCGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACG  >  NC_002947/5763490‑5763780
                                                                                                                                                    |                                                                                                                                              
cccGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAg                                                                                                                                                >  5:117515/1‑149 (MQ=255)
  cGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTg                                                                                                                                              <  1:247937/149‑1 (MQ=255)
   gTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTgg                                                                                                                                             <  3:181761/149‑1 (MQ=255)
      ccaccCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGa                                                                                                                                           <  6:117515/148‑1 (MQ=255)
           cTCGGGCCGATGGTTTCCAGTGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGt                                                                                                                                     >  1:85773/1‑149 (MQ=255)
            tCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGtt                                                                                                                                    >  5:209626/1‑149 (MQ=255)
              gggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTAt                                                                                                                                  <  5:140167/149‑1 (MQ=255)
              gggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTAt                                                                                                                                  >  8:42138/1‑149 (MQ=255)
              gggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggcgtc                                                                                                                                   <  6:191215/148‑5 (MQ=255)
              gggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGAGGCCCGCGGTAGCGACGCCCGTGCCAAGGGGGACGGGGTGAt                                                                                                                                  >  6:257113/1‑149 (MQ=255)
                gCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCTCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGa                                                                                                                                <  8:245927/149‑1 (MQ=255)
                  cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGAt                                                                                                                               >  3:83380/1‑148 (MQ=255)
                  cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggcgtc                                                                                                                                   >  1:53572/1‑140 (MQ=255)
                  cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggcgtc                                                                                                                                   <  2:53572/144‑5 (MQ=255)
                    aTGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAg                                                                                                                            >  7:243972/1‑149 (MQ=255)
                     tGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAgc                                                                                                                           >  8:261122/1‑149 (MQ=255)
                          tCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCaa                                                                                                                      >  2:223708/1‑149 (MQ=255)
                                  tGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCa                                                                                                              >  7:285302/1‑149 (MQ=255)
                                    gcgcCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCAc                                                                                                             <  8:79190/148‑1 (MQ=255)
                                        cAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAgg                                                                                                        <  3:188063/149‑1 (MQ=255)
                                        cAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggcgtc                                                                                                                                   <  1:192151/122‑5 (MQ=255)
                                        cAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggcgtc                                                                                                                                   >  2:192151/1‑118 (MQ=255)
                                           cTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAc                                                                                                     <  6:209626/149‑1 (MQ=255)
                                            tACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGt                                                                                                       >  5:137252/1‑146 (MQ=255)
                                            tACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGt                                                                                                       <  6:137252/146‑1 (MQ=255)
                                                    cgcATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACgggg                                                                                                                                      <  2:308401/107‑1 (MQ=255)
                                                    cgcATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACgggg                                                                                                                                      >  1:308401/1‑107 (MQ=255)
                                                        ttttCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCa                                                                                                                  <  8:16571/123‑1 (MQ=255)
                                                        ttttCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCa                                                                                                                  >  7:16571/1‑123 (MQ=255)
                                                        ttttCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGGTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCAcgg                                                                                        >  7:122982/1‑149 (MQ=255)
                                                           tCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCAcggcg                                                                                      <  3:97632/148‑1 (MQ=255)
                                                                tgcGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCgtg                                                                                >  4:157284/1‑149 (MQ=255)
                                                                  cgcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCAGGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                              <  8:243972/149‑1 (MQ=255)
                                                                  cgcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                              <  3:157284/149‑1 (MQ=255)
                                                                         cAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTcg                                                                       >  3:247636/1‑149 (MQ=255)
                                                                          aaaCTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTcgc                                                                      >  1:145987/1‑149 (MQ=255)
                                                                          aaaCTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTcgc                                                                      >  1:191243/1‑149 (MQ=255)
                                                                                 cttGGCCAGTCAATGGGCGATTGCTGGCTCGTTTCGGTGATGCGCGCGGTAGCGACGCCCGTGCGAAGTGGGACGGGGTTATGATCAGCGCCAAGCCAGGCAGCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTgg                                                               <  5:175150/147‑1 (MQ=255)
                                                                                     gCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTgc                                                           >  3:128780/1‑149 (MQ=255)
                                                                                       cAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTgcgc                                                         >  2:43107/1‑149 (MQ=255)
                                                                                        aGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGcgcg                                                        >  1:96493/1‑149 (MQ=255)
                                                                                              tGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCg                                                  <  2:191243/149‑1 (MQ=255)
                                                                                                gTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCAGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCggg                                                <  6:153932/149‑1 (MQ=255)
                                                                                                        cTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTgg                                         <  1:223708/148‑1 (MQ=255)
                                                                                                          ggCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCa                                      >  2:130871/1‑149 (MQ=255)
                                                                                                           gCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTc                                       <  8:122982/147‑1 (MQ=255)
                                                                                                             aggTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGAGTGGTTGCGCGGTGCCGGGCTTCTGGTCATtc                                   <  5:257113/147‑1 (MQ=255)
                                                                                                             aCGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATtc                                   <  7:42138/149‑1 (MQ=255)
                                                                                                              cGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATtct                                  <  4:83380/149‑1 (MQ=255)
                                                                                                               gTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATtctc                                 <  7:261122/149‑1 (MQ=255)
                                                                                                                tttCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCg                                <  2:85773/149‑1 (MQ=255)
                                                                                                                   cGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGAc                              <  4:247636/148‑1 (MQ=255)
                                                                                                                        aTGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                         <  8:285302/148‑1 (MQ=255)
                                                                                                                         tGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGc                        <  2:145987/148‑1 (MQ=255)
                                                                                                                            ccGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACg                    >  3:113432/1‑149 (MQ=255)
                                                                                                                            ccGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACg                    <  4:113432/149‑1 (MQ=255)
                                                                                                                                     cGACGCCCGTGCCAAGTGGGCCGGGGTTATGATCAGCGCCAACCGAGTCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGCTATTAGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGAGCATGGCACCGGATACCt             >  8:229846/1‑147 (MQ=255)
                                                                                                                                           ccGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACAGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCg                                                  <  2:232640/104‑1 (MQ=255)
                                                                                                                                           ccGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACAGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCg                                                  >  1:232640/1‑104 (MQ=255)
                                                                                                                                           ccGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTACACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTg      >  5:198370/1‑148 (MQ=255)
                                                                                                                                            cGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTa    >  6:147117/1‑149 (MQ=255)
                                                                                                                                              tGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACg  <  1:130871/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCCGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACG  >  NC_002947/5763490‑5763780

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: