Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,832,046 C→G 6.8% intergenic (+21/+131) pvdS → / ← pvdY alternative sigma factor/hydroxyproline acetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,832,0460CG 6.8% 182.7 / 4.1 59intergenic (+21/+131)pvdS/pvdYalternative sigma factor/hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base C (30/24);  new base G (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

ACCGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCT  >  NC_002947/4831901‑4832187
                                                                                                                                                   |                                                                                                                                             
aCCGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAG‑GACTGCAA‑GCAAACGCCa                                                                                                                                              >  2:230943/1‑147 (MQ=255)
   ggcTTAATGGCGTGCGGCAGAAGGACATTGCCAAGGAACTGGGCGTGTCGCCGGCGTTGGTCAAGGTGATGATCCGCGATGCGGGGGGGCGTTGCCGCTAGAGGGCGAGTCGCGGGGCCTGGAGGT‑GACTGCAA‑GCAAACGCCGGAGCg                                                                                                                                         >  8:167887/3‑149 (MQ=255)
    ccTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCg                                                                                                                                         >  2:268822/1‑148 (MQ=255)
     cTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGaa                                                                                                                                       >  5:137011/1‑149 (MQ=255)
         aTGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCt                                                                                                                                    <  7:34578/148‑1 (MQ=255)
            gCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGa                                                                                                                                 <  5:196518/148‑1 (MQ=255)
            gCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGAt                                                                                                                                <  5:93646/149‑1 (MQ=255)
                 ccGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGctc                                                                                                                           >  6:328225/1‑149 (MQ=255)
                   gCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCtctc                                                                                                                         >  6:112473/1‑149 (MQ=255)
                    cAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCg                                                                                                                        >  5:254252/1‑149 (MQ=255)
                     aGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCGAGAGCGATCCAGGTTCGCTCTGGc                                                                                                                       <  6:7615/149‑1 (MQ=39)
                        aGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTTCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCtt                                                                                                                     <  1:148688/148‑1 (MQ=255)
                          gACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATCTACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCtttt                                                                                                                   <  1:230943/149‑1 (MQ=255)
                             aTTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgctt                                                                                                               >  1:289395/1‑149 (MQ=255)
                              ttGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCAAGAGCGAACCTGGATCGCTCTCGCTTttgcttg                                                                                                              <  6:128509/149‑1 (MQ=255)
                                    aGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcgc                                                                                                        <  7:156054/149‑1 (MQ=255)
                                          tGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTAt                                                                                                  >  4:39668/1‑149 (MQ=255)
                                            ggCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcg                                                                                                         <  3:45680/140‑1 (MQ=255)
                                            ggCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcg                                                                                                         >  4:45680/1‑140 (MQ=255)
                                                  tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGtttt                                                                                          >  5:74027/1‑149 (MQ=255)
                                                       gACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCa                                                                                                            <  7:19990/126‑1 (MQ=255)
                                                       gACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GAAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCa                                                                                                            >  8:19990/1‑126 (MQ=255)
                                                         ccTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGc                                                                                        <  7:89785/144‑1 (MQ=255)
                                                         ccTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGc                                                                                        >  8:89785/1‑144 (MQ=255)
                                                          cTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATAc                                                                                  >  6:39725/1‑149 (MQ=255)
                                                              gTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGc                                                                               <  6:137011/148‑1 (MQ=255)
                                                                 aaCTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCaa                                                                            <  3:169541/148‑1 (MQ=255)
                                                                    ttCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAAc                                                                        >  3:74260/1‑149 (MQ=255)
                                                                     tCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACt                                                                       >  2:59398/1‑149 (MQ=255)
                                                                             cgcgATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAATGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAgcg                                                                >  1:256486/1‑149 (MQ=255)
                                                                                         gTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcgccg                                                    >  4:73615/1‑148 (MQ=255)
                                                                                         gTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcgccg                                                    >  1:283542/1‑148 (MQ=255)
                                                                                          tGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcgccgc                                                   >  6:222676/1‑148 (MQ=255)
                                                                                                gCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTAt                                            >  8:221352/1‑149 (MQ=255)
                                                                                                 ccGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATg                                           <  5:222676/149‑1 (MQ=255)
                                                                                                 ccGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATg                                           <  5:328225/149‑1 (MQ=255)
                                                                                                    cAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTc                                        >  3:212789/1‑149 (MQ=255)
                                                                                                              aGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgtg                              <  2:289395/149‑1 (MQ=255)
                                                                                                                   ccAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAaca                         <  5:39725/149‑1 (MQ=255)
                                                                                                                   ccAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAaca                         >  3:243830/1‑149 (MQ=255)
                                                                                                                    cAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAacac                        <  2:256486/149‑1 (MQ=255)
                                                                                                                      ggCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAAcacca                      <  5:112473/149‑1 (MQ=255)
                                                                                                                       gCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAAcaccac                     <  7:167887/149‑1 (MQ=255)
                                                                                                                       gCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAAcacca                      >  2:313951/1‑148 (MQ=255)
                                                                                                                            aaGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGGCGTTATTTGTGAAACACCACACTa                 <  3:313531/148‑1 (MQ=255)
                                                                                                                              gAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc              >  5:131002/1‑149 (MQ=255)
                                                                                                                              gAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc              >  4:116368/1‑149 (MQ=255)
                                                                                                                              gAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc              <  4:243830/149‑1 (MQ=255)
                                                                                                                               tagcgcTGCAA‑GCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGTGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCTTAAGTTATGGTCGTTATTTGTGa                             <  7:201918/129‑1 (MQ=21)
                                                                                                                               tagcgcTGCAA‑GCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGTGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTTTTTGTGa                             >  8:201918/6‑134 (MQ=21)
                                                                                                                                  gACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc            >  3:258920/1‑148 (MQ=255)
                                                                                                                                    cTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTAc         <  3:39668/149‑1 (MQ=255)
                                                                                                                                       caa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                               <  3:201696/124‑1 (MQ=255)
                                                                                                                                       caa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                               >  4:201696/1‑124 (MQ=255)
                                                                                                                                        aa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc            <  2:202745/142‑1 (MQ=255)
                                                                                                                                        aa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc            >  1:202745/1‑142 (MQ=255)
                                                                                                                                        aa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCa     <  4:258920/149‑1 (MQ=255)
                                                                                                                                        aa‑gcaaAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGt       >  4:255042/1‑147 (MQ=38)
                                                                                                                                            caaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCTCAGCGCAAAAAACGCAGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCACCTACGTCAGCt  >  4:180240/1‑149 (MQ=255)
                                                                                                                                                   cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc              >  6:286956/1‑130 (MQ=255)
                                                                                                                                                   cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc              <  5:286956/130‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ACCGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCT  >  NC_002947/4831901‑4832187

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