Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,832,068 C→G 7.8% intergenic (+43/+109) 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,0680CG 7.8% 198.6 / 7.1 64intergenic (+43/+109)pvdS/pvdYalternative sigma factor/hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base G (2/3);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCG  >  NC_002947/4831920‑4832212
                                                                                                                                                    |                                                                                                                                                
gCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCtctc                                                                                                                                                  >  3:57582/1‑149 (MQ=255)
   gAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCt                                                                                                                                               >  4:177243/1‑149 (MQ=255)
   gAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGcg                                                                                                                                               >  5:59726/1‑148 (MQ=38)
    aaGGAGATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCtt                                                                                                                                              >  7:322161/1‑149 (MQ=255)
     aGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCttt                                                                                                                                             <  7:185580/149‑1 (MQ=255)
      ggACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCgttt                                                                                                                                            <  7:42942/149‑1 (MQ=38)
        aCATTGTCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgc                                                                                                                                          >  4:170875/1‑149 (MQ=255)
            tGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgcttg                                                                                                                                       >  2:360428/1‑148 (MQ=255)
             gCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTttgcttgc                                                                                                                                      <  6:59726/148‑1 (MQ=37)
              ccAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGATGATGACTGCAAGCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCa                                                                                                                                     <  7:66085/148‑1 (MQ=37)
              ccAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAg                                                                                                                                    <  4:163406/149‑1 (MQ=255)
              ccAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAACCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCa                                                                                                                                     <  4:191208/148‑1 (MQ=255)
               cAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgc                                                                                                                                   >  2:311980/1‑149 (MQ=255)
               cAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgc                                                                                                                                   >  7:173902/1‑149 (MQ=255)
                 aGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcgc                                                                                                                                 >  1:97334/1‑149 (MQ=255)
                         ggCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATgg                                                                                                                         <  1:311980/149‑1 (MQ=255)
                               tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGtttt                                                                                                                   <  3:177243/149‑1 (MQ=255)
                               tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGtttt                                                                                                                   <  4:57582/149‑1 (MQ=255)
                                 gCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGTCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGc                                                                                                                 <  2:70567/149‑1 (MQ=255)
                                    gACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCTCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCaaa                                                                                                              >  4:206254/1‑149 (MQ=255)
                                       cTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATAc                                                                                                           <  7:145339/149‑1 (MQ=255)
                                       cTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATAc                                                                                                           <  1:360428/149‑1 (MQ=255)
                                           gTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGcc                                                                                                       <  2:97334/149‑1 (MQ=255)
                                             cAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCaa                                                                                                     >  5:190515/1‑149 (MQ=255)
                                                 ttCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGc                                                                                                        >  7:85444/1‑142 (MQ=255)
                                                 ttCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGc                                                                                                        <  8:85444/142‑1 (MQ=255)
                                                 ttCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAAc                                                                                                 >  3:157920/1‑149 (MQ=255)
                                                  tCATGATCCGCGATGCGCTGTTTGATTGCCGCATGACTGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGGCAGAGCGAAGCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTATTATGGTTTTTTGCAAATACGGCCAATAAc                                                                                                 <  8:110217/148‑1 (MQ=255)
                                                    aTGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCg                                                                                              >  5:309757/1‑149 (MQ=255)
                                                      gATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATa                                                                                                            >  4:296987/1‑133 (MQ=255)
                                                      gATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATa                                                                                                            <  3:296987/133‑1 (MQ=255)
                                                                  gcTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcg                                                                                >  3:33400/1‑149 (MQ=255)
                                                                          aTTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGt                                                                        >  4:287875/1‑149 (MQ=255)
                                                                          aTTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGt                                                                        >  5:345474/1‑149 (MQ=255)
                                                                           ttGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGt                                                                        >  1:52382/1‑148 (MQ=255)
                                                                               cGCCAGACGGCCAGTCGCCAGGCCTGACGATGACGGCAAGCAACCGCCAGACCGAACCTGGAGCGCTCTCGCTTGTGCTTGCAGCGCTCGTAGGGTGGTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATgg                                                                   <  8:322161/149‑1 (MQ=255)
                                                                                gCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGt                                                                  >  1:82927/1‑149 (MQ=255)
                                                                                  aaGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGc                                                                                                                                                <  8:307715/69‑1 (MQ=255)
                                                                                  aaGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGc                                                                                                                                                >  7:307715/1‑69 (MQ=255)
                                                                                         ccAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtg                                                         <  6:190515/149‑1 (MQ=255)
                                                                                          cAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                                                        <  3:287875/149‑1 (MQ=255)
                                                                                           aGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgtg                                                       <  3:170875/149‑1 (MQ=255)
                                                                                            gTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACgg                                                                                                         <  3:325500/98‑1 (MQ=255)
                                                                                            gTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACgg                                                                                                         >  4:325500/1‑98 (MQ=255)
                                                                                                 cAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAacac                                                 <  8:173902/149‑1 (MQ=255)
                                                                                                     ccTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACcaca                                             <  2:52382/149‑1 (MQ=255)
                                                                                                       tGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAAcacca                                               <  2:309460/145‑1 (MQ=255)
                                                                                                       tGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACcacag                                            >  1:309460/1‑147 (MQ=255)
                                                                                                            aTGACTGCAAGCAAACGCCAGAGCTAACCCGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGt                                                               <  7:168822/124‑1 (MQ=255)
                                                                                                            aTGACTGCAAGCAAACGCCAGAGCGAACCCGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCCCAGCGCAAAAATCGCCGCAACTTATGGtaga                                                               >  8:168822/1‑121 (MQ=255)
                                                                                                                 tGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACg                                 >  3:121698/1‑149 (MQ=255)
                                                                                                                 tGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACg                                 >  1:361077/1‑149 (MQ=255)
                                                                                                                  gcaagcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGt                                >  3:85971/1‑149 (MQ=255)
                                                                                                                    aagcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGt                                >  1:357260/1‑147 (MQ=255)
                                                                                                                     agcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                                       <  7:110369/139‑1 (MQ=255)
                                                                                                                     agcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                                       >  8:110369/1‑139 (MQ=255)
                                                                                                                      gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCATGTTATGGTCGTGATTAGTGAAACACCACACTAGCCCCTACGTCCGc                            >  6:196199/1‑149 (MQ=255)
                                                                                                                      gcaaAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGc                            >  4:96175/1‑149 (MQ=38)
                                                                                                                           cGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGc                        >  5:173306/1‑148 (MQ=255)
                                                                                                                              cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGc                                        >  1:299542/1‑129 (MQ=255)
                                                                                                                              cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGc                                        <  2:299542/129‑1 (MQ=255)
                                                                                                                                        cTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAAc          >  5:188217/1‑149 (MQ=255)
                                                                                                                                          ggATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCg        >  4:82923/1‑149 (MQ=255)
                                                                                                                                            aTCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCt      <  4:33400/149‑1 (MQ=255)
                                                                                                                                                ctctCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCg  >  4:48908/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCG  >  NC_002947/4831920‑4832212

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: