Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,832,060 A→T 10.4% intergenic (+35/+117) 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,0600AT10.4% 187.8 / 10.5 67intergenic (+35/+117)pvdS/pvdYalternative sigma factor/hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base T (3/4);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCT  >  NC_002947/4831912‑4832208
                                                                                                                                                    |                                                                                                                                                    
ggCGTGCCGCAGAAGGACATTGCCAAGGAACTGGCCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGa                                                                                                                                                      >  1:70567/1‑149 (MQ=255)
       cGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCGAGAGCGATCCAGGTTCGctct                                                                                                                                               >  7:44071/1‑149 (MQ=255)
        gCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCtctc                                                                                                                                              >  3:57582/1‑149 (MQ=255)
         cAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCGAGAGCGATCCAGGTTCGctctgg                                                                                                                                             <  8:44071/149‑3 (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                                                                                                               <  4:57582/149‑1 (MQ=255)
                                       tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGtttt                                                                                                               <  3:177243/149‑1 (MQ=255)
                                         gCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGTCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGc                                                                                                             <  2:70567/149‑1 (MQ=255)
                                            gACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCTCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCaaa                                                                                                          >  4:206254/1‑149 (MQ=255)
                                               cTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATAc                                                                                                       <  1:360428/149‑1 (MQ=255)
                                               cTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATAc                                                                                                       <  7:145339/149‑1 (MQ=255)
                                                   gTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGcc                                                                                                   <  2:97334/149‑1 (MQ=255)
                                                     cAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCaa                                                                                                 >  5:190515/1‑149 (MQ=255)
                                                       aCTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCg                                                                                                                                                   >  8:89974/1‑97 (MQ=255)
                                                       aCTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCg                                                                                                                                                   <  7:89974/97‑1 (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                                                                                                        <  3:296987/133‑1 (MQ=255)
                                                              gATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATa                                                                                                        >  4:296987/1‑133 (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                                                                                                                                            >  7:307715/1‑69 (MQ=255)
                                                                                          aaGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGc                                                                                                                                            <  8:307715/69‑1 (MQ=255)
                                                                                                 ccAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtg                                                     <  6:190515/149‑1 (MQ=255)
                                                                                                  cAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                                                    <  3:287875/149‑1 (MQ=255)
                                                                                                   aGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgtg                                                   <  3:170875/149‑1 (MQ=255)
                                                                                                    gTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACgg                                                                                                     >  4:325500/1‑98 (MQ=255)
                                                                                                    gTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACgg                                                                                                     <  3:325500/98‑1 (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                                   >  8:110369/1‑139 (MQ=255)
                                                                                                                             agcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                                   <  7:110369/139‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
GGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCT  >  NC_002947/4831912‑4832208

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: