Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,610,296 C→G 7.3% I261M (ATC→ATG PP_4081 → type IV/VI secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,610,2960CG 7.3% 170.0 / 3.9 55I261M (ATC→ATGPP_4081type IV/VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base G (2/2);  total (26/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 = 9.55e-01

CACTGATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCTCG  >  NC_002947/4610148‑4610443
                                                                                                                                                     |                                                                                                                                                   
cACTGATTCACCAACAGCACGCGATCATCCGCCCACTCCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCCTGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCcgccatc                                                                                                                                                     <  7:129795/149‑1 (MQ=18)
  cTGATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCcgccatcgc                                                                                                                                                   <  2:82310/149‑1 (MQ=255)
       tCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATgg                                                                                                                                              >  2:271931/1‑149 (MQ=255)
        cACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATgg                                                                                                                                              <  7:111406/148‑1 (MQ=255)
         aCCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGCCCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCg                                                                                                                                            >  1:107775/1‑149 (MQ=255)
         aCCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCg                                                                                                                                            >  2:75457/1‑149 (MQ=255)
         aCCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCg                                                                                                                                            <  2:225532/149‑1 (MQ=255)
            aaCTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTg                                                                                                                                         >  4:107211/1‑149 (MQ=255)
            aaCTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTg                                                                                                                                         <  1:8686/149‑1 (MQ=255)
             aCTGCACGGGAGCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCACCGCCATGGCGATGGCggaga                                                                                                                                        <  5:47678/149‑4 (MQ=18)
                    gggATCATCCGCGAACTGCGCGGCCCCACCCCCGCAGAAGGCTGCGCCCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCTTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACgg                                                                                                                                 <  7:250108/149‑1 (MQ=255)
                          atcCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTa                                                                                                                           >  5:87600/1‑149 (MQ=255)
                           tcCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGTGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTa                                                                                                                           >  1:118665/1‑148 (MQ=255)
                           tcCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTAc                                                                                                                          <  6:87600/149‑1 (MQ=255)
                           tcCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTAc                                                                                                                          <  1:39062/149‑1 (MQ=255)
                             cgcgAACTGCGCGGCCCCACCCCCGAAGCAGTCTGCGACCCCTACACCAACGTCCCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCAGCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAg                                                                                                                        >  8:8021/1‑149 (MQ=255)
                                        cgGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGc                                                                                                             <  2:3949/149‑1 (MQ=255)
                                             ccACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCt                                                                                                         >  3:154265/1‑148 (MQ=255)
                                              cACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGa                                                                                                       <  6:382097/149‑1 (MQ=255)
                                              cACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCATGCACCTGa                                                                                                       >  3:248199/1‑149 (MQ=255)
                                                  cccGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATcac                                                                                                   <  3:107211/149‑1 (MQ=255)
                                                   ccGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATcac                                                                                                   >  2:366891/1‑148 (MQ=255)
                                                        gaagTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGa                                                                                             <  5:73952/149‑1 (MQ=255)
                                                         aagTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAg                                                                                            <  1:75457/149‑1 (MQ=255)
                                                               gCGACCCCTACACCAACGTCGCCACCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTAATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTg                                                                                      >  6:83175/1‑149 (MQ=255)
                                                                cGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTg                                                                                      <  1:366891/148‑1 (MQ=255)
                                                                    cccTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTc                                                                                 >  6:350782/1‑149 (MQ=255)
                                                                     ccTACACCAACGTCGCCCCAAGAAACTTCCGCATGAGCCGGATATCTACGTGGTGGAGCCCGTTGGTCATCTCCGACATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCATGCACATGATCACAACCGAGGTGCTGGAGTCg                                                                                >  3:214963/1‑149 (MQ=17)
                                                                         cacCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCt                                                                              <  5:240289/147‑1 (MQ=255)
                                                                           ccAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAAc                                                                          <  1:271931/149‑1 (MQ=255)
                                                                            cAACGTCGCCCCCAGAAACTTCCGCATGAGCGCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAAc                                                                          >  1:229209/1‑149 (MQ=255)
                                                                            cAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACg                                                                         >  6:143720/1‑149 (MQ=255)
                                                                                      cccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgca                                                               <  2:322762/149‑1 (MQ=255)
                                                                                      cccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgca                                                               <  4:259427/149‑1 (MQ=255)
                                                                                         aGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgcagca                                                            <  3:297555/149‑1 (MQ=255)
                                                                                             aCTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAg                                                                                                                        <  4:179469/85‑1 (MQ=255)
                                                                                             aCTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAg                                                                                                                        >  3:179469/1‑85 (MQ=255)
                                                                                                         gCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTca                                            <  2:118665/149‑1 (MQ=255)
                                                                                                          ccGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTcac                                           <  5:143720/149‑1 (MQ=255)
                                                                                                                    gCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACt                                  <  4:154265/149‑1 (MQ=255)
                                                                                                                          ggggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTc                            <  1:43235/147‑1 (MQ=39)
                                                                                                                           tgtaGCCCGTTGGTCATCTCCGCGATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTTGAGTCGCTGAACGGCATTTTGAAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCa                           <  4:214963/146‑1 (MQ=255)
                                                                                                                           tggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCa                           <  2:229209/149‑1 (MQ=39)
                                                                                                                               gCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACg                        >  8:199527/1‑148 (MQ=39)
                                                                                                                                 ccGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAAt                     >  5:48066/1‑149 (MQ=39)
                                                                                                                                  cGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAAtt                    >  5:284483/1‑149 (MQ=38)
                                                                                                                                   gTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTc                   >  8:202576/1‑149 (MQ=38)
                                                                                                                                       gTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGa               <  1:373426/149‑1 (MQ=37)
                                                                                                                                       gTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGa               <  2:107775/149‑1 (MQ=37)
                                                                                                                                        tCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAg              >  3:153563/1‑149 (MQ=35)
                                                                                                                                          gccaccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAAc                         >  6:375181/5‑136 (MQ=12)
                                                                                                                                          gccaccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAAc                         <  5:375181/132‑1 (MQ=12)
                                                                                                                                          aTCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTTTCACAGCACCACTGAGTTCAACGAATTCGCGa               >  1:160696/1‑146 (MQ=21)
                                                                                                                                          aTCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGa               <  2:160696/146‑1 (MQ=34)
                                                                                                                                           tctcCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCAGTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGtt           >  7:78566/1‑149 (MQ=18)
                                                                                                                                            caccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGttt          >  2:173133/3‑149 (MQ=12)
                                                                                                                                             tcCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGtttt         >  2:38717/1‑149 (MQ=32)
                                                                                                                                               ggccaccgccatGGCGATGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTg       >  5:281077/2‑149 (MQ=17)
                                                                                                                                                  catcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCgt                                                >  2:151031/1‑105 (MQ=32)
                                                                                                                                                  catcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCgt                                                <  1:151031/105‑1 (MQ=32)
                                                                                                                                                   atcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCTc   >  4:370025/1‑149 (MQ=32)
                                                                                                                                                    tcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCTCg  >  8:238420/1‑149 (MQ=32)
                                                                                                                                                     |                                                                                                                                                   
CACTGATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCTCG  >  NC_002947/4610148‑4610443

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