Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,000,135 C→T 6.6% Q274* (CAG→TAG) ‡ 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,1350CT 6.6% 172.9 / 2.3 60Q274* (CAG→TAG) ‡PP_3528ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base T (2/2);  total (30/30)
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.69e-01

CTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTACC  >  NC_002947/3999987‑4000282
                                                                                                                                                     |                                                                                                                                                   
cTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCAccc                                                                                                                                                     <  3:46642/149‑1 (MQ=255)
   cTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAg                                                                                                                                                   <  4:203583/148‑1 (MQ=255)
    tACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGt                                                                                                                                                 >  7:29359/1‑149 (MQ=255)
        tGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTg                                                                                                                                             >  7:219607/1‑149 (MQ=255)
         gTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTg                                                                                                                                             >  3:238793/1‑148 (MQ=255)
            cacaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAg                                                                                                                                         >  4:20989/1‑149 (MQ=255)
              cagccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                       <  5:19880/149‑1 (MQ=255)
              cagccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                       <  2:181531/149‑1 (MQ=255)
               agccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTgg                                                                                                                                      <  2:171526/149‑1 (MQ=255)
               agccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTgg                                                                                                                                      >  5:141649/1‑149 (MQ=255)
                 ccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGcc                                                                                                                                    <  1:45674/149‑1 (MQ=255)
                   agcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCAGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAg                                                                                                                                  >  6:138887/1‑149 (MQ=255)
                    gcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGc                                                                                                                                 >  2:261080/1‑149 (MQ=255)
                     cTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCg                                                                                                                                <  8:29359/149‑1 (MQ=255)
                          ccGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACg                                                                                                                            <  8:125550/148‑1 (MQ=255)
                                  cACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTa                                                                                                                   <  3:105566/149‑1 (MQ=255)
                                              ctgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACACCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGc                                                                                                       <  4:308928/149‑1 (MQ=255)
                                                   aaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCa                                                                                                  <  5:193070/149‑1 (MQ=255)
                                                   aaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCa                                                                                                  <  4:238793/149‑1 (MQ=255)
                                                       ttcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCAcccc                                                                                              >  1:136040/1‑149 (MQ=255)
                                                        tcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCAAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCg                                                                                             >  3:282161/1‑149 (MQ=255)
                                                               cGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCg                                                                                      >  7:166424/1‑149 (MQ=255)
                                                               cGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCg                                                                                      >  3:86708/1‑149 (MQ=255)
                                                                  gCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTg                                                                                   <  2:136040/149‑1 (MQ=255)
                                                                     gggCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGcc                                                                                 >  7:65372/1‑148 (MQ=255)
                                                                        cAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGaac                                                                             <  3:20989/149‑1 (MQ=255)
                                                                        cAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGaac                                                                             <  4:86708/149‑1 (MQ=255)
                                                                           ggTCTGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACACCGCAGATCGTTGCCGaacaa                                                                           >  6:64252/1‑148 (MQ=255)
                                                                               aGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAg                                                                        >  2:60926/1‑147 (MQ=255)
                                                                                   aGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctg                                                                  <  1:60926/149‑1 (MQ=255)
                                                                                       ggCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctggct                                                               <  5:138887/148‑1 (MQ=255)
                                                                                       ggCCAACGCCAAGCCTGCCAAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctggct                                                               <  4:282161/148‑1 (MQ=255)
                                                                                         ccAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCa                                                                                           <  6:94450/118‑1 (MQ=255)
                                                                                         ccAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCa                                                                                           >  5:94450/1‑118 (MQ=255)
                                                                                          cAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGaca                                                           >  6:199439/1‑149 (MQ=255)
                                                                                              gcaaaGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCt                                                        >  1:228729/4‑149 (MQ=255)
                                                                                                 aaGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTc                                                      >  2:259755/1‑147 (MQ=255)
                                                                                                        ccGAAGTGGCCGCGCTGCAGGCACCGCAGAT‑GGGGATCAGCACCCAGGTGCTGGAAGAGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCAGCAGATCGTTGCCGAACAACAGAACCTGGCTGACACCTTCCAggccctg                                             >  8:223797/1‑149 (MQ=255)
                                                                                                           aaGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgggc                                          >  7:98733/1‑149 (MQ=255)
                                                                                                                         ctgGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCaa                            <  8:166424/149‑1 (MQ=255)
                                                                                                                           ggCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          >  8:32911/1‑149 (MQ=255)
                                                                                                                            gCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          >  2:227415/1‑148 (MQ=255)
                                                                                                                            gCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          <  8:98733/148‑1 (MQ=255)
                                                                                                                            gCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGc                          >  8:170856/1‑148 (MQ=255)
                                                                                                                            gCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCATGATTCCACGCAAGc                          >  7:237256/1‑148 (MQ=255)
                                                                                                                              aCCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTa                        >  3:174384/1‑148 (MQ=255)
                                                                                                                                  cAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGt                   >  1:114127/1‑149 (MQ=255)
                                                                                                                                   aGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTg                  <  7:170856/149‑1 (MQ=255)
                                                                                                                                    gCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgg                                                  >  6:188613/1‑116 (MQ=255)
                                                                                                                                    gCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgg                                                  <  5:188613/116‑1 (MQ=255)
                                                                                                                                     cT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTg                  <  8:65372/147‑1 (MQ=255)
                                                                                                                                       ccacttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTAc                       >  6:306774/8‑141 (MQ=255)
                                                                                                                                       ccacttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTAc                       <  5:306774/134‑1 (MQ=255)
                                                                                                                                         gggATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCg              <  5:199439/148‑1 (MQ=255)
                                                                                                                                         acttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGa                                                             <  6:235830/96‑1 (MQ=255)
                                                                                                                                         acttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGa                                                             >  5:235830/6‑101 (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)
                                                                                                                                                     |                                                                                                                                                   
CTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCT‑GGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTACC  >  NC_002947/3999987‑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: ‑