Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 590,711 T→A 6.5% W279R (TGG→AGG)  PP_0500 → dTDP‑4‑rhamnose reductase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947590,7110TA 6.5% 185.8 / 3.0 61W279R (TGG→AGG) PP_0500dTDP‑4‑rhamnose reductase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base A (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTCGACGA  >  NC_002947/590571‑590857
                                                                                                                                             |                                                                                                                                                   
ccggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCt                                                                                                                                             <  4:136588/149‑1 (MQ=255)
 cggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGgcacgcg                                                                                                                                             >  3:98639/1‑148 (MQ=255)
 cggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTg                                                                                                                                            <  3:339974/149‑1 (MQ=255)
       aTACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGc                                                                                                                                      <  6:330022/149‑1 (MQ=255)
        tACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGcc                                                                                                                                     <  3:77470/149‑1 (MQ=255)
         aCCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCt                                                                                                                                                     <  7:102655/132‑1 (MQ=255)
         aCCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCt                                                                                                                                                     >  8:102655/1‑132 (MQ=255)
          ccGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGGCACGCGGCTTgata                                                                                                                                    <  4:98639/148‑4 (MQ=255)
          ccGCAAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                                   >  5:184638/1‑149 (MQ=255)
           cGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCAcc                                                                                                                                  >  7:31847/1‑149 (MQ=255)
           cGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCAcc                                                                                                                                  >  7:37188/1‑149 (MQ=255)
              cAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACct                                                                                                                                 <  8:332309/147‑1 (MQ=255)
               aGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTa                                                                                                                              >  4:139601/1‑149 (MQ=255)
                gCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                             <  2:316253/149‑1 (MQ=255)
                gCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                             <  3:139601/149‑1 (MQ=255)
                     ccGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGAc                                                                                                                        >  1:170527/1‑149 (MQ=255)
                     ccGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGAc                                                                                                                        >  8:306515/1‑149 (MQ=255)
                      cGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCa                                                                                                                                    <  5:211079/136‑1 (MQ=255)
                      cGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCa                                                                                                                                    >  6:211079/1‑136 (MQ=255)
                       gTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTgg                                                                                                                          >  4:170492/1‑145 (MQ=255)
                       gTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTgg                                                                                                                          <  3:170492/145‑1 (MQ=255)
                        ttCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCgg                                                                                                                     >  4:190644/1‑149 (MQ=255)
                             gCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTa                                                                                                               <  5:349495/149‑1 (MQ=255)
                             gCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTa                                                                                                               <  6:124670/149‑1 (MQ=255)
                              cccccACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                                   >  6:186191/1‑129 (MQ=255)
                                 ccACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGc                                                                                                           <  2:170527/149‑1 (MQ=255)
                                  cACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGc                                                                                                           >  7:236551/1‑148 (MQ=255)
                                  cACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGc                                                                                                           >  5:224074/1‑148 (MQ=255)
                                    cgccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCAt                                                                                                        <  5:129357/149‑1 (MQ=255)
                                    cgccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCAt                                                                                                        <  3:190644/149‑1 (MQ=255)
                                    cgccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCAt                                                                                                        <  8:236551/149‑1 (MQ=255)
                                            gCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACg                                                                                                >  3:69146/1‑149 (MQ=255)
                                             cccATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGc                                                                                               <  8:31847/149‑1 (MQ=255)
                                             cccATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGc                                                                                               <  8:37188/149‑1 (MQ=255)
                                                           caggATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCa                                                                                 >  2:266550/3‑149 (MQ=255)
                                                                tGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACcggc                                                                            <  6:224074/149‑1 (MQ=255)
                                                                               acaGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                             >  5:23392/1‑149 (MQ=255)
                                                                                caGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                             <  1:266550/149‑1 (MQ=255)
                                                                                    aCGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTc                                                        <  4:69146/149‑1 (MQ=255)
                                                                                     cgcgGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTcc                                                       >  6:320932/1‑149 (MQ=255)
                                                                                             tGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACcg                                                                              <  6:169651/118‑1 (MQ=255)
                                                                                             tGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACcg                                                                              >  5:169651/1‑118 (MQ=255)
                                                                                                   gCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                             >  3:49166/1‑129 (MQ=255)
                                                                                                   gCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                             <  4:49166/129‑1 (MQ=255)
                                                                                                        aaaaTCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCgttgt                                    <  5:320932/149‑1 (MQ=255)
                                                                                                            tCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGa                                >  4:284159/1‑149 (MQ=255)
                                                                                                               ttCATACCTTCGGTATCAAGCCGCGTGCCGGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAt                                                              <  6:248709/116‑1 (MQ=255)
                                                                                                                 cATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACaaa                            <  6:184638/148‑1 (MQ=255)
                                                                                                                        tCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCt                                                  <  4:39861/119‑1 (MQ=255)
                                                                                                                        tCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCt                                                  >  3:39861/1‑119 (MQ=255)
                                                                                                                        tCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGc                    >  6:363587/1‑149 (MQ=255)
                                                                                                                          ggTATCTAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTGCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTCGACAAAGGCTACGcc                   <  6:192725/148‑1 (MQ=255)
                                                                                                                            tATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTa                >  5:29205/1‑149 (MQ=255)
                                                                                                                            tATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTa                >  1:272088/1‑149 (MQ=255)
                                                                                                                              tCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTCATCACCGCCGGTGCCGGCTTCATTGTCTCCCACCTGTGCCATGCGTTGTTGGAc                               >  4:50798/1‑132 (MQ=255)
                                                                                                                              tCAAGCCGCGTGCCTGGCGCGCTGGCTTGCAACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGAc                               <  3:50798/132‑1 (MQ=255)
                                                                                                                              gcAAGCCAGCGCGCCAGGCACGCGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACaaa                            <  6:237632/134‑1 (MQ=255)
                                                                                                                              gcAAGCCAGCGCGCCAGGCACGCGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACaaa                            >  5:237632/2‑135 (MQ=255)
                                                                                                                                  gCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCAtt          >  7:171571/1‑149 (MQ=255)
                                                                                                                                        gTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAgg                          >  6:5109/1‑128 (MQ=255)
                                                                                                                                        gTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAgg                          <  5:5109/128‑1 (MQ=255)
                                                                                                                                         tGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTcgac    <  5:363587/149‑1 (MQ=255)
                                                                                                                                           ccTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTcgacga  >  4:61764/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCGCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTCGACGA  >  NC_002947/590571‑590857

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: