Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,793,968 A→T 11.3% intergenic (‑129/+27) pvdI ← / ← PP_4222 non‑ribosomal peptide synthetase/protein SyrP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,793,9680AT11.3% 166.1 / 10.3 62intergenic (‑129/+27)pvdI/PP_4222non‑ribosomal peptide synthetase/protein SyrP
Reads supporting (aligned to +/- strand):  ref base A (32/23);  new base T (4/3);  total (36/26)
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.56e-01

TGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGC  >  NC_002947/4793822‑4794115
                                                                                                                                                  |                                                                                                                                                    
tgagtCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGTGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAg                                                                                                                                                    >  5:242950/1‑149 (MQ=255)
                aTGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGTGCTTGCGGCTTGCCGCAGGCGCTCGCca                                                                                                                                   <  5:281697/149‑1 (MQ=39)
                      ggggTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTg                                                                                                                              <  3:67972/148‑1 (MQ=255)
                                 tGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAgcgc                                                                                                                  >  4:181000/1‑149 (MQ=255)
                                 tGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAgcgc                                                                                                                  >  7:111462/1‑149 (MQ=255)
                                 tGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGcgcc                                                                                                                 >  4:128790/1‑149 (MQ=255)
                                  gATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCGCAAAGcgcc                                                                                                                 >  7:252939/1‑149 (MQ=255)
                                       gCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTc                                                                                                            >  6:135110/1‑149 (MQ=255)
                                         gctttggCCGCACGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgc                                                                                                          >  4:81233/6‑149 (MQ=255)
                                         aGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgc                                                                                                          >  2:26822/1‑149 (MQ=255)
                                          gCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgcg                                                                                                         <  7:4036/149‑1 (MQ=255)
                                           cTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCgg                                                                                                        >  6:308533/1‑149 (MQ=255)
                                            tCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGt                                                                                                       >  5:294109/1‑149 (MQ=255)
                                              ggCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTg                                                                                                      <  5:100777/148‑1 (MQ=255)
                                               gCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGcc                                                                                                    >  1:177031/1‑149 (MQ=255)
                                                      ggCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgcg                                                                                                         >  5:87059/1‑137 (MQ=255)
                                                      ggCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgcg                                                                                                         <  6:87059/137‑1 (MQ=255)
                                                        ccAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCAt                                                                                          <  3:128790/149‑1 (MQ=255)
                                                                   ccAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGAc                                                                                 <  8:111462/148‑1 (MQ=255)
                                                                    cAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCa                                                                               >  4:245010/1‑149 (MQ=255)
                                                                     aGTAGAACGAGCGGAGGGAGGGCGAATTTAAAGCGGCAGCGCCATGCCCGGACGCTGCGATGGCGAGCGCCTGAGGCAAGCCGCAAGCAATAGCCACAGTGGCCTCAAAGCGACGCTTCGCGGTGCCACTGAGCAATGGCCAAGACCAc                                                                              >  6:216791/1‑149 (MQ=255)
                                                                       tAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACtt                                                                            <  3:181000/149‑1 (MQ=255)
                                                                          aaCGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGc                                                                         >  1:80982/1‑149 (MQ=255)
                                                                              aGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTgg                                                                     <  2:177031/149‑1 (MQ=255)
                                                                              aGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTgg                                                                     >  7:71714/1‑149 (MQ=255)
                                                                               gCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGc                                                                    >  7:161953/1‑149 (MQ=255)
                                                                                cGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGcc                                                                   <  6:294109/149‑1 (MQ=255)
                                                                                   agggagggCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCgg                                                                <  6:242950/149‑1 (MQ=255)
                                                                                      gagggCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTg                                                              >  1:166242/1‑148 (MQ=255)
                                                                                          gCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATg                                                                                         >  7:131685/1‑117 (MQ=255)
                                                                                          gCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATg                                                                                         <  8:131685/117‑1 (MQ=255)
                                                                                           cGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGt                                                        >  1:241047/1‑149 (MQ=255)
                                                                                            gAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACg                                                          <  5:288829/146‑1 (MQ=255)
                                                                                            gAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACg                                                          >  6:288829/1‑146 (MQ=255)
                                                                                                  aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGTCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCggg                                                 >  3:63139/1‑149 (MQ=255)
                                                                                                  aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCgg                                                  <  2:80982/148‑1 (MQ=255)
                                                                                                  aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCgg                                                  <  5:308533/148‑1 (MQ=255)
                                                                                                      cgGCAGCGCCATGCCCGGACGCTGCGCTGTCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTg                                             <  4:63139/149‑1 (MQ=255)
                                                                                                       ggCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTgc                                            <  8:71714/149‑1 (MQ=255)
                                                                                                        gCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTgc                                            >  1:67442/1‑149 (MQ=255)
                                                                                                        gCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACttt                                                                           >  4:168321/1‑117 (MQ=255)
                                                                                                               cATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAg                                    >  1:41631/1‑149 (MQ=255)
                                                                                                                   cccGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCAtg                                <  2:32045/149‑1 (MQ=255)
                                                                                                                         cgctgcgctgGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCg                          <  5:135110/149‑1 (MQ=255)
                                                                                                                          gctgcgctgGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGa                         <  8:161953/149‑1 (MQ=255)
                                                                                                                          gccactgtgGCGAGTGCTTGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTc                           >  3:254486/8‑147 (MQ=38)
                                                                                                                            cactgtgGCGAGTGCTTGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAg                        >  3:242919/6‑148 (MQ=38)
                                                                                                                                  ggCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGc                 <  3:245010/149‑1 (MQ=255)
                                                                                                                                   gCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCa                >  4:94625/1‑149 (MQ=255)
                                                                                                                                     gAGTGCTTGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGcgccgc                                       <  2:116583/124‑1 (MQ=38)
                                                                                                                                     gAGTGCTTGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAcg              >  1:2278/1‑149 (MQ=255)
                                                                                                                                     gAGTGCTTGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAcg              >  3:27512/1‑149 (MQ=255)
                                                                                                                                     gAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGc                                                      >  6:283532/1‑109 (MQ=255)
                                                                                                                                     gAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGc                                                      <  5:283532/109‑1 (MQ=255)
                                                                                                                                          cttGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGtt         <  4:254486/147‑1 (MQ=255)
                                                                                                                                           cTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGttt        >  2:134426/1‑149 (MQ=255)
                                                                                                                                              cGGCAAGCCTAAAGAACTCGCCACAGTGGCCTCAAAGCGCAGCTTCGAGGTGACCCTGAGCAATGACCACGACCACTTTGCGTGGCCCGGTTAACGGTGCCCGGGCTTGCGCCGCGAGCATGTTGTAGAGAATCAGCACGTCTTTTTcc     >  6:22808/1‑149 (MQ=255)
                                                                                                                                              cGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTcc     >  5:66551/1‑149 (MQ=255)
                                                                                                                                               ggCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCt    <  3:81233/149‑1 (MQ=255)
                                                                                                                                               ggCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCt    <  2:41631/149‑1 (MQ=255)
                                                                                                                                                gCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTg   >  5:146481/1‑149 (MQ=255)
                                                                                                                                                 cAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGc  >  3:171174/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGC  >  NC_002947/4793822‑4794115

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: