Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,832,046 C→G 11.2% 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,0460CG11.2% 199.4 / 13.7 70intergenic (+21/+131)pvdS/pvdYalternative sigma factor/hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base G (4/4);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGC  >  NC_002947/4831902‑4832190
                                                                                                                                                |                                                                                                                                                
ccGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCaga                                                                                                                                               <  6:256273/148‑1 (MQ=255)
 cGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGTCCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGc                                                                                                                                             <  8:137664/149‑1 (MQ=255)
 cGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAAAGCGAGAGc                                                                                                                                             >  8:42942/1‑149 (MQ=255)
  gCCTTAATGGCGTGCCGCAGAAGGACATTCCCAACGAACTCGGAGTATCGCCGACCTAGGTCAACTTCATCATCCGCCATGCGCTGGTGCATTCCCGCAACACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCACAGc                                                                                                                                             >  7:110217/1‑148 (MQ=255)
       aaTGGCGTGCCGCAGAAGGACATTGCCAAGTAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAAcc                                                                                                                                        <  5:292460/148‑1 (MQ=255)
       aaTGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCt                                                                                                                                       >  8:185580/1‑149 (MQ=255)
          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                                                                                                             >  7:173902/1‑149 (MQ=255)
                                 cAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgc                                                                                                             >  2:311980/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                                                                                     <  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)
                                                                 aCTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCg                                                                                                                                 <  7:89974/97‑1 (MQ=255)
                                                                 aCTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCg                                                                                                                                 >  8:89974/1‑97 (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                                                  >  5:345474/1‑149 (MQ=255)
                                                                                            aTTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGt                                                  >  4:287875/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                                                                                   <  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                  <  2:299542/129‑1 (MQ=255)
                                                                                                                                                cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGc                  >  1:299542/1‑129 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCGCCTTAATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGC  >  NC_002947/4831902‑4832190

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: