Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,455,129 A→G 6.9% V155A (GTG→GCG)  PP_4792 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,455,1290AG 6.9% 163.4 / 2.8 58V155A (GTG→GCG) PP_4792hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (2/2);  total (29/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 = 8.82e-01

CCCGCCTGTTCAAAACTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCC  >  NC_002947/5454982‑5455275
                                                                                                                                                   |                                                                                                                                                  
cccGCCTGTTCAAAACTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCAc                                                                                                                                                   <  4:143695/149‑1 (MQ=255)
cccGCCTGTTCAAAACTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCCAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCAc                                                                                                                                                   >  5:357358/1‑149 (MQ=255)
     cTGTTCAAAACTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGa                                                                                                                                              <  4:84052/149‑1 (MQ=255)
      tGTTCAAAACTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCa                                                                                                                                                    >  5:87064/1‑142 (MQ=255)
      tGTTCAAAACTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCa                                                                                                                                                    <  6:87064/142‑1 (MQ=255)
           aaaaCTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCg                                                                                                                                         >  6:254606/1‑148 (MQ=255)
             aaCTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCt                                                                                                                                                 >  2:265704/1‑138 (MQ=255)
             aaCTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGACAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCt                                                                                                                                                 <  1:265704/138‑1 (MQ=255)
                tCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCt                                                                                                                                   >  3:130667/1‑149 (MQ=255)
                  cAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCtt                                                                                                                                  >  1:340215/1‑148 (MQ=255)
                  cAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCt                                                                                                                                   <  5:254606/147‑1 (MQ=255)
                  cAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTc                                                                                                                                 <  4:130667/149‑1 (MQ=255)
                         aGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAg                                                                                                                          <  5:201262/149‑1 (MQ=255)
                                       cAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCGCCTCGACAGCGAAATCTTc                                                                                                                                 >  6:296189/1‑128 (MQ=255)
                                       cAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCGCCTCGACAGCGAAATCTTc                                                                                                                                 <  5:300932/128‑1 (MQ=255)
                                       cAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCGCCTCGACAGCGAAATCTTc                                                                                                                                 <  5:296189/128‑1 (MQ=255)
                                       cAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCGCCTCGACAGCGAAATCTTc                                                                                                                                 >  6:300932/1‑128 (MQ=255)
                                       cAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAg                                                                                                            >  3:234928/1‑149 (MQ=255)
                                              gggTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAgcctcgc                                                                                                      >  7:233395/1‑148 (MQ=255)
                                               ggTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAgcctcgc                                                                                                      <  4:278526/147‑1 (MQ=255)
                                                   aCCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCt                                                                                                >  8:250102/1‑149 (MQ=255)
                                                      aGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAAt                                                                                                                                     >  4:111910/1‑109 (MQ=255)
                                                      aGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAAt                                                                                                                                     <  3:111910/109‑1 (MQ=255)
                                                       gCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGggc                                                                                            <  7:250102/149‑1 (MQ=255)
                                                           caccCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGgcggcg                                                                                        <  4:234928/149‑1 (MQ=255)
                                                              ccGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGc                                                                                     <  8:233395/149‑1 (MQ=255)
                                                                   gAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCctcgcctcgc                                                                                                 <  6:164319/132‑1 (MQ=255)
                                                                   gAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCctcgcctcgc                                                                                                 >  5:164319/1‑132 (MQ=255)
                                                                    aaaTCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCACCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAAt                                                                               <  2:340215/149‑1 (MQ=255)
                                                                         ccAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGc                                                                             <  4:15497/146‑1 (MQ=255)
                                                                         ccAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGc                                                                             >  3:15497/1‑146 (MQ=255)
                                                                         ccAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCaa                                                                           <  2:24807/148‑1 (MQ=255)
                                                                          cAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAAcc                                                                         <  6:357358/149‑1 (MQ=255)
                                                                                 tGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAAt                                                                               <  1:173539/136‑1 (MQ=255)
                                                                                 tGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAAt                                                                               >  2:173539/1‑136 (MQ=255)
                                                                                          ttGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAgcctcgcc                                                                                                     <  6:308406/105‑1 (MQ=255)
                                                                                          ttGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAgcctcgcc                                                                                                     >  5:308406/1‑105 (MQ=255)
                                                                                                      tGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCg                                              >  5:56444/1‑148 (MQ=255)
                                                                                                       gAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCg                                              >  4:291607/1‑147 (MQ=255)
                                                                                                                cccAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGcc                                   >  5:311248/1‑149 (MQ=255)
                                                                                                                cccAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGc                                    <  3:291607/148‑1 (MQ=255)
                                                                                                                 ccAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGc                                    >  7:243638/1‑147 (MQ=255)
                                                                                                                  cAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTgg                                                      <  8:292946/128‑1 (MQ=255)
                                                                                                                  cAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTgg                                                      >  7:292946/1‑128 (MQ=255)
                                                                                                                  cAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAg                                 >  8:250868/1‑149 (MQ=255)
                                                                                                                    ggTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGcc                               >  4:285953/1‑149 (MQ=255)
                                                                                                                      tCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAAt                                                                               <  1:300003/99‑1 (MQ=255)
                                                                                                                      tCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAAt                                                                               >  2:300003/1‑99 (MQ=255)
                                                                                                                        gTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCTGGCCCGAAGCCTTGGCCAGCcgcgc                           <  4:346554/149‑1 (MQ=255)
                                                                                                                        gTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCcgcgc                           <  6:311248/149‑1 (MQ=255)
                                                                                                                        gTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCcgcgc                           >  3:43349/1‑149 (MQ=255)
                                                                                                                              cTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTACCGCCATCAGCCTCGCCTTGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCaaa                     >  8:168692/1‑149 (MQ=255)
                                                                                                                                 gCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAgcg                  <  7:250868/149‑1 (MQ=255)
                                                                                                                                  ccTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAgcg                  >  6:342388/1‑148 (MQ=255)
                                                                                                                                   ctctTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCcgcgcg                          >  6:155837/1‑139 (MQ=255)
                                                                                                                                         gcgcCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGcc                               <  8:103459/128‑1 (MQ=255)
                                                                                                                                         gcgcCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGcc                               >  7:103459/1‑128 (MQ=255)
                                                                                                                                          cgcCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTc         <  8:243638/149‑1 (MQ=255)
                                                                                                                                              aTTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGt     >  8:236666/1‑149 (MQ=255)
                                                                                                                                                  cACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgcc  <  4:43349/148‑1 (MQ=255)
                                                                                                                                                   aCCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCAt               <  7:38588/134‑1 (MQ=255)
                                                                                                                                                   aCCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCAt               >  8:38588/1‑134 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCGCCTGTTCAAAACTCCAGCGTGAGCGCTGCATGTGCCAGGCCTGGGTCACCAGCACCACCCGGCGAAATCCCAAAGGCTGCAGCACCTTGGCCGTCAGCTGAGCATTTTCCCAGGTCGTGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCC  >  NC_002947/5454982‑5455275

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