Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,193,900 G→T 6.6% L324M (CTG→ATG)  PP_2802 ← amino acid permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,193,9000GT 6.6% 193.9 / 4.0 61L324M (CTG→ATG) PP_2802amino acid permease family protein
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base T (2/2);  total (32/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.90e-01

AAGGCGCACAGGTTGACGGCGGTGAAGGCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAACAGGTCGCCG  >  NC_002947/3193752‑3194042
                                                                                                                                                    |                                                                                                                                              
aaGGCGCACAGGTTGACGGCGGTGAAGGCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAg                                                                                                                                                <  3:297234/149‑1 (MQ=255)
     ggacaGGTTGACGGCGGTGGAGGCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGGTGGGATCGTCAGGCCAAACCCGGCGGTGCCGCAAAAACCGCCCCGCAAGATGTTGAAGGGCGGGGGCTTGTAGCGCGGGGGCAgcgcgc                                                                                                                                           >  8:288910/3‑149 (MQ=255)
           gTTGACGGCGGTGAAGGCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACa                                                                                                                                     >  2:182590/1‑149 (MQ=255)
                      tGAAGGCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCgggg                                                                                                                           >  4:229236/1‑148 (MQ=255)
                          ggCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCATCGCGCCGAACAGCCGCCGGGGCAgtt                                                                                                                      >  3:190946/1‑149 (MQ=255)
                          ggCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCATCGCGCCGAACAGCCGCCGGGGCAgtt                                                                                                                      >  2:28320/1‑149 (MQ=255)
                              agaAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAgttggtt                                                                                                                  <  7:288910/149‑1 (MQ=255)
                              agaAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAgttggtt                                                                                                                  >  6:352334/1‑149 (MQ=255)
                                 aaGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTc                                                                                                               <  4:28044/149‑1 (MQ=255)
                                 aaGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTc                                                                                                               <  1:182590/149‑1 (MQ=255)
                                  aGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTcc                                                                                                              >  1:35797/1‑149 (MQ=255)
                                    gcgcCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTGTTGTAGCGCAGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCg                                                                                                            <  6:176398/149‑1 (MQ=255)
                                        cGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGCCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAgttgg                                                                                                                    >  6:260998/1‑137 (MQ=255)
                                        cGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGCCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAgttgg                                                                                                                    <  5:260998/137‑1 (MQ=255)
                                            gTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAggg                                                                                                    >  2:130514/1‑149 (MQ=255)
                                               gACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCAAAGGGct                                                                                                  >  6:261584/1‑147 (MQ=255)
                                                aCCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTa                                                                                                <  3:229236/149‑1 (MQ=255)
                                                  cAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTaga                                                                                              <  3:230434/149‑1 (MQ=255)
                                                  cAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTa                                                                                                >  8:286827/1‑147 (MQ=255)
                                                     ggAAGTCGCGGTTGCGATCGTCAGGCCAAACGCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCATCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGagca                                                                                           <  4:190946/149‑1 (MQ=255)
                                                      gAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCATCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGagcag                                                                                          <  1:28320/149‑1 (MQ=255)
                                                       aaGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGagcagc                                                                                         <  1:338872/149‑1 (MQ=255)
                                                            gcgGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGcc                                                                                    >  2:283500/1‑149 (MQ=255)
                                                                    gATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTgg                                                                            >  5:249710/1‑149 (MQ=255)
                                                                     aTCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGc                                                                           <  6:208990/149‑1 (MQ=255)
                                                                           aGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGaa                                                                     <  5:352334/149‑1 (MQ=255)
                                                                           aGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGaa                                                                     <  6:249710/149‑1 (MQ=255)
                                                                             gCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTa                                                                                                >  6:4459/1‑120 (MQ=255)
                                                                             gCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAAt                                                                    >  5:42749/1‑148 (MQ=255)
                                                                                    ccGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATccc                                                            >  5:302968/1‑149 (MQ=255)
                                                                                             ccGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCg                                                     >  3:182173/1‑147 (MQ=255)
                                                                                               gAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGt                                                 >  8:61564/1‑149 (MQ=255)
                                                                                               gAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGt                                                 >  3:119742/1‑149 (MQ=255)
                                                                                                aaaaCACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAgttggtt                                                                                                                  >  1:177505/1‑83 (MQ=255)
                                                                                                aaaaCACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAgttggtt                                                                                                                  <  2:177505/83‑1 (MQ=255)
                                                                                                     acCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGcga                                           >  3:190242/1‑149 (MQ=255)
                                                                                                     acCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGCcga                                           >  5:268838/1‑149 (MQ=255)
                                                                                                      ccGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGcgac                                          >  3:151873/1‑149 (MQ=255)
                                                                                                          acaGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACgatg                                      <  1:130514/149‑1 (MQ=255)
                                                                                                          acaGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACgat                                       >  4:359088/1‑148 (MQ=255)
                                                                                                          acaGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACgat                                       <  3:359088/148‑1 (MQ=255)
                                                                                                             gCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATc                                   <  3:251431/149‑1 (MQ=255)
                                                                                                              cAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATcc                                  >  5:179226/1‑149 (MQ=255)
                                                                                                              cAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATcc                                  <  1:39552/149‑1 (MQ=255)
                                                                                                              cAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATcc                                  >  2:39552/1‑149 (MQ=255)
                                                                                                                       gAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAgc                                                                                                                                   >  7:327261/1‑43 (MQ=255)
                                                                                                                       gAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAgc                                                                                                                                   <  8:327261/43‑1 (MQ=255)
                                                                                                                       gAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAgc                                                                                                                                   <  8:322367/43‑1 (MQ=255)
                                                                                                                       gAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAgc                                                                                                                                   >  7:322367/1‑43 (MQ=255)
                                                                                                                        aaGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGa                         <  3:332440/148‑1 (MQ=255)
                                                                                                                             cGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACgg                   >  7:93933/1‑149 (MQ=255)
                                                                                                                              ggggTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGc                  <  4:190242/149‑1 (MQ=255)
                                                                                                                                 gTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCgaa               <  6:302968/149‑1 (MQ=255)
                                                                                                                                 gTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCgaa               <  2:35797/149‑1 (MQ=255)
                                                                                                                                  tCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGATTGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGCCCACGGCg                 >  3:45429/1‑146 (MQ=255)
                                                                                                                                      gTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAAc           <  6:42749/148‑1 (MQ=255)
                                                                                                                                       tAGCGCGGGTGCACGGCGCCTAACCGCCGCCGGGGATTTTGTTTGTCCGGGCCAAGGGCTTATTGCAGCCGGCCAACGCTGGCCTGTAATCACACTCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAACa          <  1:283500/148‑1 (MQ=255)
                                                                                                                                         gcgcGGGTGCAGCGCGCCGAACAGCCGCCTGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAACAGGt       >  4:155625/1‑149 (MQ=255)
                                                                                                                                         gcgcGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGTGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAACAGGt       >  1:169600/1‑149 (MQ=255)
                                                                                                                                         gcgcGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCgaagaa            <  4:100896/144‑1 (MQ=255)
                                                                                                                                         gcgcGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCgaagaa            >  3:100896/1‑144 (MQ=255)
                                                                                                                                            cgGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAACAGGTcgc    <  6:179226/149‑1 (MQ=255)
                                                                                                                                              ggTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAACAGGTcgccg  <  8:151695/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AAGGCGCACAGGTTGACGGCGGTGAAGGCGAGAAAGGCGCCGAAGTTGACCAGGGAAGTCGCGGTTGCGATCGTCAGGCCAAACCCGGCGGTGCCGAAAACACCGCACAGCAAGATGTTGAAGGCCGGGGTCTTGTAGCGCGGGTGCAGCGCGCCGAACAGCCGCCGGGGCAGTTGGTTGTCCCGGCCAAGGGCGTAGAGCAGCCGGCCCACGCTGGCCTGGAATGACATCCCCGCAGCGAAGTGGGCGACGATGATCCCGGTCAGGACCACGGCGAAGAACAGGTCGCCG  >  NC_002947/3193752‑3194042

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: ‑