Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 994,951 T→A 7.8% intergenic (+27/+98) PP_0859 → / ← PP_0860 ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947994,9510TA 7.8% 331.2 / 18.2 115intergenic (+27/+98)PP_0859/PP_0860ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (55/51);  new base A (4/5);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

Only 100 of 117 total aligned reads displayed.

CCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGCA  >  NC_002947/994804‑995100
                                                                                                                                                   |                                                                                                                                                     
ccTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCCTCCTTGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGCCAGCAGGCCCGGCTTCTTCGCGGGCTTGc                                                                                                                                                     <  1:1091105/150‑1 (MQ=255)
    cTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCcgc                                                                                                                                                 >  2:214658/1‑150 (MQ=255)
    cTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCcgc                                                                                                                                                 <  1:1118730/150‑1 (MQ=255)
       aGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:633684/150‑1 (MQ=255)
       aGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  2:480232/150‑1 (MQ=255)
        gTGCGGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                             <  2:1155158/150‑1 (MQ=255)
            cGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCGTTCTCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggc                                                                                                                                         <  1:1423825/150‑1 (MQ=255)
              gggAGTCGGACGGGGGGTTCACTGCAGTACTGGACGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAACCTTGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccg                                                                                                                                       <  2:786868/150‑1 (MQ=255)
              gggAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTCCCCGGCATATCTGGATGCCGATACCTTCTAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccg                                                                                                                                       <  1:235648/150‑1 (MQ=255)
               ggAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccgg                                                                                                                                      <  2:1695437/150‑1 (MQ=255)
                gAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggc                                                                                                                                     >  1:738586/1‑150 (MQ=255)
                         ggggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGcc                                                                                                                                                    <  2:796832/126‑1 (MQ=255)
                         ggggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGcc                                                                                                                                                    >  1:796832/1‑126 (MQ=255)
                         ggggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGACCGCGAAAGGGCCGGCCCTgccagc                                                                                                                            >  1:885664/1‑150 (MQ=255)
                            gtgtTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTGCCGGGCCAAGTTCCCGGCCAATCTGGAAGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGACCTGCCAGCCaa                                                                                                                         >  1:1848491/1‑150 (MQ=255)
                             tgtTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCACGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAAt                                                                                                                        >  1:1095556/1‑150 (MQ=255)
                                  aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCa                                                                                                                   >  2:441201/1‑150 (MQ=255)
                                  aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCa                                                                                                                   >  2:649216/1‑150 (MQ=255)
                                  aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              >  2:1341907/1‑123 (MQ=255)
                                  aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:1341907/123‑1 (MQ=255)
                                    tGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCAcc                                                                                                                 >  1:949045/1‑150 (MQ=255)
                                      cAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCgg                                                                                                               >  2:456934/1‑150 (MQ=255)
                                         tACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCg                                                                                                            >  2:606377/1‑150 (MQ=255)
                                            tGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTc                                                                                                         >  1:1188858/1‑150 (MQ=255)
                                            tGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTc                                                                                                         <  1:1714882/150‑1 (MQ=255)
                                                 gcggcgGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGt                                                                                                    <  1:606377/150‑1 (MQ=255)
                                                   ggcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGTCCGGCCCTCCCAGCCAACTTCCACCGGCCGTTACCCGttt                                                                                                  >  2:148526/1‑149 (MQ=255)
                                                    gcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgc                                                                                                 <  1:1543385/150‑1 (MQ=255)
                                                      ggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTGCACCGTTgccg                                                                                               >  1:213272/1‑150 (MQ=255)
                                                        aaCTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgcc                                                                                             <  1:333060/150‑1 (MQ=255)
                                                             gcggcgTACCGGGCAAAGTTCACGGCCAATCTGGATGTCGATACCTTCGAGCTGCACTTACAGCAGGCCCTGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCAc                                                                                        <  2:213272/150‑1 (MQ=38)
                                                               ggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCAc                                                                                                       >  1:1443288/1‑133 (MQ=255)
                                                                gcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGt                                                                                     >  2:690416/1‑150 (MQ=255)
                                                                  gTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTAt                                                                                   <  2:1027998/150‑1 (MQ=255)
                                                                    aCCGGGCCAAGTTCCCGGCCAATCTGGAAGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTc                                                                                 <  2:1848491/150‑1 (MQ=255)
                                                                        ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgcc                                                                                                <  2:928600/131‑1 (MQ=255)
                                                                        ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgcc                                                                                                >  1:928600/1‑131 (MQ=255)
                                                                        ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                             <  1:1513366/150‑1 (MQ=255)
                                                                        ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                             >  2:1469392/1‑150 (MQ=255)
                                                                         gCCAAGTTCCCGGCCGATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCggg                                                                            <  2:189796/150‑1 (MQ=255)
                                                                           cAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCGGGCCCGGCTTCTTCGCGGGATTGCCCGCGAAAGGGCCGGCCCTGCCAGCCCATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCTTAGTCGCGGTTa                                                                          >  1:702969/1‑150 (MQ=38)
                                                                           cAAGTTCCCGGCCAATCTGGATGCCGATACCTTAGAGCTGCACTGACAGCAGACCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTa                                                                          >  1:1494925/1‑150 (MQ=255)
                                                                              gTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAc                                                                       <  1:1228980/150‑1 (MQ=255)
                                                                                 cccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgc                                                                                                                                <  1:1347895/90‑1 (MQ=255)
                                                                                 cccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgc                                                                                                                                >  2:1347895/1‑90 (MQ=255)
                                                                                  ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGc                                                                   >  2:713459/1‑150 (MQ=255)
                                                                                      ccAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCCCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCt                                                               >  2:1462007/1‑150 (MQ=255)
                                                                                       cAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                              >  2:1578477/1‑150 (MQ=255)
                                                                                          tCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacac                                                           >  1:381522/1‑150 (MQ=255)
                                                                                                 gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                                    <  1:1578477/150‑1 (MQ=255)
                                                                                                 gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                                    <  2:1502547/150‑1 (MQ=255)
                                                                                                   cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGa                                                      >  2:290992/1‑146 (MQ=255)
                                                                                                   cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGa                                                      <  1:290992/146‑1 (MQ=255)
                                                                                                   cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                  >  2:651298/1‑150 (MQ=255)
                                                                                                     aTACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGACAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGCATTCGCGGGTAAACCCGCTCCTCCACGCCGATCAAtc                                                >  2:922617/1‑149 (MQ=255)
                                                                                                       aCCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATaca                                              <  2:381522/150‑1 (MQ=255)
                                                                                                         cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                 >  1:682590/1‑129 (MQ=255)
                                                                                                          ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCAcc                                                                                                      <  2:906248/91‑1 (MQ=255)
                                                                                                              aGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCg                                       >  1:713674/1‑150 (MQ=255)
                                                                                                               gCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGAAAACCCGCTCCTACACACCGATCAATACACATCTCGc                                      >  1:1638086/1‑150 (MQ=255)
                                                                                                                cTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                                     >  1:1586783/1‑150 (MQ=255)
                                                                                                                 tGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCg                                    >  2:869382/1‑150 (MQ=255)
                                                                                                                   cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGat                                  >  1:554527/1‑150 (MQ=255)
                                                                                                                    aCTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                  <  2:1211612/133‑1 (MQ=255)
                                                                                                                    aCTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                  >  1:1211612/1‑133 (MQ=255)
                                                                                                                    aCTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGACGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCGACACACCGATCAATACACATCTCGCCGata                                 >  2:1446921/1‑150 (MQ=255)
                                                                                                                      tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTc                                                                                 <  2:392367/100‑1 (MQ=255)
                                                                                                                      tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTc                                                                                 >  1:392367/1‑100 (MQ=255)
                                                                                                                      tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacacac                                                         <  1:1541952/124‑1 (MQ=255)
                                                                                                                      tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                             >  1:1528320/1‑136 (MQ=255)
                                                                                                                      tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                             <  2:1528320/136‑1 (MQ=255)
                                                                                                                           gcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCa                          >  2:1201489/1‑150 (MQ=255)
                                                                                                                              ggCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATcc                       <  1:412022/150‑1 (MQ=255)
                                                                                                                               gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccg                                                                                                                                       >  2:1217151/1‑37 (MQ=255)
                                                                                                                               gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCt                      >  1:1727733/1‑150 (MQ=255)
                                                                                                                                cccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCTTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc                     <  1:516639/150‑1 (MQ=255)
                                                                                                                                cccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc                     <  1:1201489/150‑1 (MQ=255)
                                                                                                                                 ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                 >  2:1571334/1‑105 (MQ=255)
                                                                                                                                 ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGa                                                      >  2:1936474/1‑116 (MQ=255)
                                                                                                                                 ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGa                                                      <  1:1936474/116‑1 (MQ=255)
                                                                                                                                  cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGACGATATCGCCCATCCTCGa                   <  2:949045/150‑1 (MQ=255)
                                                                                                                                   ggCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                  >  1:455106/1‑150 (MQ=255)
                                                                                                                                   ggCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                  <  2:888487/150‑1 (MQ=255)
                                                                                                                                   ggCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                  >  2:1229876/1‑150 (MQ=255)
                                                                                                                                    gCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                 <  2:546753/102‑1 (MQ=255)
                                                                                                                                    gCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                 >  1:546753/1‑102 (MQ=255)
                                                                                                                                    gCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATc                 <  2:738586/150‑1 (MQ=255)
                                                                                                                                      ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                 <  1:1595706/100‑1 (MQ=255)
                                                                                                                                      ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                 >  2:1595706/1‑100 (MQ=255)
                                                                                                                                       ctttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAc              >  1:1557830/3‑150 (MQ=255)
                                                                                                                                        tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAct             >  1:1431434/2‑150 (MQ=255)
                                                                                                                                        tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAct             >  2:59056/2‑150 (MQ=255)
                                                                                                                                        cttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATGCACATCTCGCCGATATCGCCCATCCTCGATCAAct             >  1:1571457/1‑150 (MQ=255)
                                                                                                                                         ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActt            >  2:381760/1‑149 (MQ=255)
                                                                                                                                         ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg            >  1:1667323/1‑150 (MQ=255)
                                                                                                                                         ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg            <  1:1737148/150‑1 (MQ=255)
                                                                                                                                         ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg            <  2:554527/150‑1 (MQ=255)
                                                                                                                                              gggCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTa       <  2:1095556/150‑1 (MQ=255)
                                                                                                                                               ggCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTAc      <  2:1586783/150‑1 (MQ=255)
                                                                                                                                                   tGCCCGCGAAAGGGCCGGCCATGACACCCCAATCCCACCGACCTTTCACCCTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTAcagca  <  2:1667323/150‑1 (MQ=37)
                                                                                                                                                   |                                                                                                                                                     
CCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGCA  >  NC_002947/994804‑995100

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: