Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,000,136 A→G 6.2% Q274R (CAG→CGG) ‡ PP_3528 → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,000,1360AG 6.2% 192.3 / 2.3 64Q274R (CAG→CGG) ‡PP_3528ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (2/2);  total (33/31)
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.45e-01

CTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTACC  >  NC_002947/3999990‑4000282
                                                                                                                                                   |                                                                                                                                                  
cTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAg                                                                                                                                                   <  4:203583/148‑1 (MQ=255)
 tACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGt                                                                                                                                                 >  7:29359/1‑149 (MQ=255)
     tGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTg                                                                                                                                             >  7:219607/1‑149 (MQ=255)
      gTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTg                                                                                                                                             >  3:238793/1‑148 (MQ=255)
         cacaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAg                                                                                                                                         >  4:20989/1‑149 (MQ=255)
           cagccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                       <  5:19880/149‑1 (MQ=255)
           cagccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                       <  2:181531/149‑1 (MQ=255)
            agccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTgg                                                                                                                                      >  5:141649/1‑149 (MQ=255)
            agccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTgg                                                                                                                                      <  2:171526/149‑1 (MQ=255)
              ccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGcc                                                                                                                                    <  1:45674/149‑1 (MQ=255)
                agcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCAGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAg                                                                                                                                  >  6:138887/1‑149 (MQ=255)
                 gcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGc                                                                                                                                 >  2:261080/1‑149 (MQ=255)
                  cTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCg                                                                                                                                <  8:29359/149‑1 (MQ=255)
                       ccGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACg                                                                                                                            <  8:125550/148‑1 (MQ=255)
                               cACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTa                                                                                                                   <  3:105566/149‑1 (MQ=255)
                                           ctgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACACCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGc                                                                                                       <  4:308928/149‑1 (MQ=255)
                                                aaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCa                                                                                                  <  4:238793/149‑1 (MQ=255)
                                                aaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCa                                                                                                  <  5:193070/149‑1 (MQ=255)
                                                    ttcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCAcccc                                                                                              >  1:136040/1‑149 (MQ=255)
                                                     tcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCAAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCg                                                                                             >  3:282161/1‑149 (MQ=255)
                                                            cGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCg                                                                                      >  3:86708/1‑149 (MQ=255)
                                                            cGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCg                                                                                      >  7:166424/1‑149 (MQ=255)
                                                               gCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTg                                                                                   <  2:136040/149‑1 (MQ=255)
                                                                  gggCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGcc                                                                                 >  7:65372/1‑148 (MQ=255)
                                                                     cAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGaac                                                                             <  3:20989/149‑1 (MQ=255)
                                                                     cAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGaac                                                                             <  4:86708/149‑1 (MQ=255)
                                                                        ggTCTGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACACCGCAGATCGTTGCCGaacaa                                                                           >  6:64252/1‑148 (MQ=255)
                                                                            aGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAg                                                                        >  2:60926/1‑147 (MQ=255)
                                                                                aGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctg                                                                  <  1:60926/149‑1 (MQ=255)
                                                                                    ggCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctggct                                                               <  5:138887/148‑1 (MQ=255)
                                                                                    ggCCAACGCCAAGCCTGCCAAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctggct                                                               <  4:282161/148‑1 (MQ=255)
                                                                                      ccAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCa                                                                                           <  6:94450/118‑1 (MQ=255)
                                                                                      ccAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCa                                                                                           >  5:94450/1‑118 (MQ=255)
                                                                                       cAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGaca                                                           >  6:199439/1‑149 (MQ=255)
                                                                                           gcaaaGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCt                                                        >  1:228729/4‑149 (MQ=255)
                                                                                              aaGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTc                                                      >  2:259755/1‑147 (MQ=255)
                                                                                                     ccGAAGTGGCCGCGCTGCAGGCACCGCAGATGGGGATCAGCACCCAGGTGCTGGAAGAGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCAGCAGATCGTTGCCGAACAACAGAACCTGGCTGACACCTTCCAggccctg                                             >  8:223797/1‑149 (MQ=255)
                                                                                                        aaGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgggc                                          >  7:98733/1‑149 (MQ=255)
                                                                                                                      ctgGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCaa                            <  8:166424/149‑1 (MQ=255)
                                                                                                                        ggCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          >  8:32911/1‑149 (MQ=255)
                                                                                                                         gCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          <  8:98733/148‑1 (MQ=255)
                                                                                                                         gCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          >  8:170856/1‑148 (MQ=255)
                                                                                                                         gCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          >  2:227415/1‑148 (MQ=255)
                                                                                                                         gCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCATGATTCCACGCAAGc                          >  7:237256/1‑148 (MQ=255)
                                                                                                                           aCCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTa                        >  3:174384/1‑148 (MQ=255)
                                                                                                                               cAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGt                   >  1:114127/1‑149 (MQ=255)
                                                                                                                                aGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTg                  <  7:170856/149‑1 (MQ=255)
                                                                                                                                 gCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgg                                                  >  6:188613/1‑116 (MQ=255)
                                                                                                                                 gCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgg                                                  <  5:188613/116‑1 (MQ=255)
                                                                                                                                  cTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTg                  <  8:65372/147‑1 (MQ=255)
                                                                                                                                    ccacttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTAc                       <  5:306774/134‑1 (MQ=255)
                                                                                                                                    ccacttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTAc                       >  6:306774/8‑141 (MQ=255)
                                                                                                                                      gggATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCg              <  5:199439/148‑1 (MQ=255)
                                                                                                                                      acttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGa                                                             >  5:235830/6‑101 (MQ=255)
                                                                                                                                      acttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGa                                                             <  6:235830/96‑1 (MQ=255)
                                                                                                                                          tCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCg         <  1:227415/149‑1 (MQ=255)
                                                                                                                                           cAGGACCCAGGTGCTTGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCg         <  7:223797/148‑1 (MQ=255)
                                                                                                                                           cAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCaa                            <  6:83243/129‑1 (MQ=255)
                                                                                                                                           cAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCaa                            >  5:83243/1‑129 (MQ=255)
                                                                                                                                            aGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCgg        <  2:228729/148‑1 (MQ=255)
                                                                                                                                                cccAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTAc   >  1:57561/1‑149 (MQ=255)
                                                                                                                                                 ccAGGTGCTGGACGTGGCCAGCTAACGGCGCAACTACAACGCCGTCGATATCACCCCGCAGATCGTTGCCTAAATACAGAAGCTGGCTGACACCTTCCAgggcctgggc                                          <  1:47601/109‑1 (MQ=255)
                                                                                                                                                 ccAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgggc                                          >  2:47601/1‑109 (MQ=255)
                                                                                                                                                 ccAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTAcc  >  1:219122/1‑149 (MQ=255)
                                                                                                                                                   aGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTc      <  5:318834/143‑1 (MQ=255)
                                                                                                                                                   aGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTc      >  6:318834/1‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTACC  >  NC_002947/3999990‑4000282

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: