Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,000,085 C→T 6.8% T242T (ACC→ACT 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,0850CT 6.8% 165.9 / 2.9 60T242T (ACC→ACTPP_3528ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

GTGACCCTTTGGGACCCGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATC  >  NC_002947/3999939‑4000229
                                                                                                                                                   |                                                                                                                                                
gTGACCCTTTGGGACCCGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACcttc                                                                                                                                               <  2:826481/150‑1 (MQ=255)
     ccTTTGGGACCCGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAGGTATCTTCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGc                                                                                                                                          >  1:29218/1‑150 (MQ=255)
          gggACCCGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGc                                                                                                                                     >  1:1117968/1‑150 (MQ=255)
           ggACCCGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCt                                                                                                                                    >  2:586065/1‑150 (MQ=255)
            gACCCGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTg                                                                                                                                   >  1:292139/1‑150 (MQ=255)
             aCCCTTTCCTCACCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGACGCTTGTCCAACAACCGCACCTTCTACCTGTCCACATACAGCTTTGCCGACCAGCACCGTGCTTTGCGTAAAACCTTCTTTGCCGTGCTgt                                                                                                                                  <  2:216515/150‑2 (MQ=255)
               ccGTTCCTCGCCGCCGCAGGAAGCCAGCCAGTGCTGAAGAACCTGCGTGACGGCAGCGGCCTGTCCAACGACCGCACCTTCTACCTGTCCACAGCCAGCTATGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGc                                                                                                                                >  2:570656/1‑150 (MQ=255)
               ccGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGc                                                                                                                                >  1:52699/1‑150 (MQ=255)
                              gcAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTgg                                                                                                                 <  2:52699/150‑1 (MQ=255)
                                    agccagccAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAAc                                                                                                           <  2:327817/150‑1 (MQ=255)
                                         gccagTGCTGAAGAACCTGCGCGCCGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGACGACCAGCACCGTGCTCTGCTGAAATTGTTCGTTTCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCaa                                                                                                      <  1:570656/150‑1 (MQ=255)
                                         gccagGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCCCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCaa                                                                                                      >  2:263523/1‑150 (MQ=255)
                                              ggCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTg                                                                                                 >  2:762301/1‑150 (MQ=255)
                                                 tGAAAAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCCTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGTGCCAACGCCAAGCCTGCCg                                                                                              >  2:299575/1‑150 (MQ=255)
                                                   aagaaCCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCTCCTTCTACCTGTCCACAGCCAGCTTTCTGCGCTGTCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGaa                                                                                            <  1:466338/150‑1 (MQ=255)
                                                   aagaaCCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGc                                                                                                                       <  1:575306/123‑1 (MQ=255)
                                                   aagaaCCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGc                                                                                                                       >  2:575306/1‑123 (MQ=255)
                                                    agaaCCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAg                                                                                           <  2:1091314/150‑1 (MQ=255)
                                                     gaaCCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGt                                                                                          >  1:264807/1‑150 (MQ=255)
                                                     gaaCCTGCGCGACGGCAGCGGCCTGTCCAACAACAACAACTTAAACCTGGCCACAGCCAGCTTTGCCGACCAGCACCGTCACATGTCCAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCCAGCCTGCCGAAGt                                                                                          >  1:523312/1‑150 (MQ=255)
                                                         ctgagcgACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCt                                                                                                  <  1:354093/134‑1 (MQ=255)
                                                         cTGAGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCt                                                                                                  >  2:354093/1‑138 (MQ=255)
                                                          tGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCcg                                                                                     >  2:538282/1‑150 (MQ=255)
                                                             gcgACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCgcgc                                                                                  >  2:686577/1‑150 (MQ=255)
                                                                actgCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCgctgc                                                                               <  1:43491/147‑1 (MQ=255)
                                                                       cGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCAcc                                                                        <  1:686577/150‑1 (MQ=255)
                                                                        ggCCTGTCTAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCGAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCg                                                                       >  2:504255/1‑150 (MQ=255)
                                                                          ccTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTgg                                                                                                                 >  1:1080307/1‑106 (MQ=255)
                                                                          ccTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTgg                                                                                                                 <  2:1080307/106‑1 (MQ=255)
                                                                          ccTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCa                                                                     <  2:323777/150‑1 (MQ=255)
                                                                                 aacaacCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACTTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTgggg                                                              >  1:272944/1‑150 (MQ=255)
                                                                                 aacaacCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACTTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTgggg                                                              >  1:272975/1‑150 (MQ=255)
                                                                                 aacaacCGCACCTTCTACCTGTCCACAGCAAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTgggg                                                              >  1:421524/1‑150 (MQ=255)
                                                                                  acaacCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGa                                                             <  2:482670/150‑1 (MQ=255)
                                                                                      ccGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGc                                                                                                    >  1:1083586/1‑107 (MQ=255)
                                                                                      ccGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGc                                                                                                    <  2:1083586/107‑1 (MQ=255)
                                                                                            cTTATACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCa                                                   >  1:714317/1‑150 (MQ=255)
                                                                                             ttCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACTTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAg                                                  <  2:272944/150‑1 (MQ=255)
                                                                                             ttCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACTTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAg                                                  <  2:272975/150‑1 (MQ=255)
                                                                                                tACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTg                                               <  1:803675/150‑1 (MQ=255)
                                                                                                   cTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTg                                            <  1:212623/150‑1 (MQ=255)
                                                                                                     gTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCTGGGTGCTgat                                          <  1:586065/150‑3 (MQ=255)
                                                                                                        cacaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGt                                       >  2:115524/1‑150 (MQ=255)
                                                                                                         acaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGc                                                                                                               >  1:268908/1‑77 (MQ=255)
                                                                                                         acaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGc                                                                                                               <  2:268908/77‑1 (MQ=255)
                                                                                                         acaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGc                                              <  2:368688/142‑1 (MQ=255)
                                                                                                         acaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGc                                              >  1:368688/1‑142 (MQ=255)
                                                                                                              cagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTg                                                                                                                  >  1:618559/1‑69 (MQ=255)
                                                                                                              cagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGCGCAAGGTCAGCCAGTg                                                                                                                  <  2:618559/69‑1 (MQ=255)
                                                                                                                   ttGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAAc                            <  1:762301/150‑1 (MQ=255)
                                                                                                                     gCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACgg                          <  1:707067/150‑1 (MQ=255)
                                                                                                                          ccAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCaa                     <  2:421524/150‑1 (MQ=255)
                                                                                                                           cAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGt                                                                                          >  2:983052/1‑80 (MQ=255)
                                                                                                                           cAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGt                                                                                          <  1:983052/80‑1 (MQ=255)
                                                                                                                                  cGTGCTCTGCTGCAACCCCTCTTTGCCGAGCTTGGTGGGGGCAGCCAGTGTGCCAACGCCAAGATACACTCAAGGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCGCCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAAc              <  1:299575/150‑1 (MQ=255)
                                                                                                                                        ctgctgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCTGGGTGCTgat                                          >  1:880847/1‑114 (MQ=255)
                                                                                                                                        ctgctgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCTGGGTGCTgat                                          <  2:880847/116‑3 (MQ=255)
                                                                                                                                            tgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTa    >  1:685751/1‑150 (MQ=255)
                                                                                                                                            tgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTa    >  2:755725/1‑150 (MQ=255)
                                                                                                                                              aaaaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATc  <  2:1117968/150‑1 (MQ=255)
                                                                                                                                              aaaaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCGGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGCAAGTGGCCAGCGCACGGCGCAACTACAACGCCGTCGCTATc  >  2:229911/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTGACCCTTTGGGACCCGTTCCTCGCCGCCGCAGAAAGCCAGCCAGGGCTGAAGAACCTGCGCGACGGCAGCGGCCTGTCCAACAACCGCACCTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATC  >  NC_002947/3999939‑4000229

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: