Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 278,742:1 +A 100% intergenic (‑237/+77) PP_0224 ← / ← sctC DszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947278,7421.A93.3% 187.8 / 6.6 60intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base A (26/30);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCG  >  NC_002947/278598‑278883
                                                                                                                                                   |                                                                                                                                                 
tCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGttaac                                                                                                                                                >  3:14694/1‑145 (MQ=255)
 cGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                               <  5:227547/148‑3 (MQ=255)
      cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTc                                                                                                                                         <  5:208857/149‑1 (MQ=255)
      cAGCTGATCGCGACAGACGTGTGCAGGTCCGGTGGTCCGGGATAAATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGc                                                                                                                                           <  7:87043/147‑1 (MQ=255)
       atcTGTTCGCGACAGAGGTTTGCTGGTCCGGGTGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGTAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                        <  4:34118/147‑1 (MQ=255)
         cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATCGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                      <  8:164208/149‑1 (MQ=255)
             tCGCGACAGACGTGTGCTGGTCCTGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa‑‑caca                                                                                                                               <  1:359612/149‑1 (MQ=255)
                cgACAGACGAGTCCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTAATGAACTAACTTTTTCGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACAt                                                                                                                            <  7:111412/149‑1 (MQ=255)
                     gACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAAACGCTCCTA‑‑CAAACATGCCt                                                                                                                        >  7:293987/1‑148 (MQ=255)
                       cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCaa                                                                                                                     <  8:78202/149‑1 (MQ=255)
                           tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATc                                                                                                                 >  3:223095/1‑149 (MQ=255)
                              tggtccggtggtccggTATATATACACTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCtg                                                                                                               <  6:21493/148‑1 (MQ=255)
                                gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCtgt                                                                                                              <  5:238568/147‑1 (MQ=255)
                                 tccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCa                                                                                                           <  6:31765/149‑1 (MQ=255)
                                  ccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTGTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACCCATGCCTCAAGATATGTTCAc                                                                                                          >  4:87213/1‑149 (MQ=255)
                                          ccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTA‑‑AAAACATGCCTAAAGATATGTGCAATGTAGGAg                                                                                                  >  3:221139/1‑149 (MQ=255)
                                             gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                               >  8:284472/1‑149 (MQ=255)
                                               ccgatctcccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                                             <  7:247014/142‑1 (MQ=255)
                                                  tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTTAAACCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTa                                                                                          >  7:334606/1‑149 (MQ=255)
                                                     aCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAccc                                                                                       >  5:66423/1‑149 (MQ=255)
                                                       ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑AACACATGCCTCAAGATCAGTGCACTGTAGGAGCCGGTTTaacagc                                                                                     >  7:103081/1‑144 (MQ=255)
                                                       ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCAA‑‑CACACATGCCTCAAGATCTGAGAACTCTAGGAGCGGGt                                                                                             >  8:247014/1‑141 (MQ=255)
                                                        cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                     >  6:329697/1‑148 (MQ=255)
                                                           aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTa                                                                                          <  1:309500/140‑1 (MQ=255)
                                                           aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTa                                                                                          >  2:309500/1‑140 (MQ=255)
                                                           aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGa                                                                                   >  2:107655/1‑147 (MQ=255)
                                                           aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCGCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGa                                                                                   >  2:103215/1‑147 (MQ=255)
                                                               tttATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                             >  3:233529/1‑149 (MQ=255)
                                                                 tataAATTTTTTTTCAACTAAGTTTTTGGAATATTTTTCCAAGGGTCTTTTGTTACTCTGTGGCGGGGCTGTTCGCGGTTAAACCCGTTCCTA‑‑CCCAATTGCTTCATCAGATGTCCACTGTAGTCGCGGGTTTACCCGCGAAAGGGcc                                                                            <  8:103081/148‑1 (MQ=255)
                                                                 tataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                            <  5:329697/148‑1 (MQ=255)
                                                                  ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa‑‑c                                                                                                                                   >  8:181996/1‑93 (MQ=255)
                                                                  ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa‑‑c                                                                                                                                   <  7:181996/93‑1 (MQ=255)
                                                                  ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                           >  6:353040/1‑148 (MQ=255)
                                                                    aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTGCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGa                                                                         >  3:355870/1‑148 (MQ=255)
                                                                    aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACACCTACTA‑‑AAAACAATCATAAAGATCTATGCACTGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGaa                                                                        >  7:226229/1‑149 (MQ=255)
                                                                        tttttATGAACTAACTTTTTGGAAAATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                               <  5:38526/122‑1 (MQ=255)
                                                                        tttttATGAACTAACTTTTTGGAAAATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                               >  6:38526/1‑122 (MQ=255)
                                                                         ttttATGAAATAACTTTTTGGAAAATTTATAGAAGGGGGATTTGTTACAGAGGGGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACCCACACATCCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAg                                                                     <  1:107655/149‑1 (MQ=38)
                                                                          tttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTATCTACACATGCCTCAAGATCTGTGCACTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAg                                                                     <  1:103215/148‑1 (MQ=255)
                                                                           ttattaACAAAATTTTTGAAATTTTTAAACATGGGTCATTTGTTACCTCGGTGCCGGGCTCTTCGCGTTAAAACCCGCCCCTA‑‑CACACATCCCTCAATATCTTTCCACTTTTGCCGCGGTTTAACCCGCGAAAGGGCCGGAACAGGtt                                                                  <  4:14694/148‑1 (MQ=255)
                                                                            tATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGttt                                                                 <  6:66423/148‑1 (MQ=255)
                                                                                 actaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCg                                                            <  5:353040/148‑1 (MQ=255)
                                                                                   taattTTTTGTAATATTTATATAAGGGGGATTTGTTACAGTGGGGCCGGGCTCTTCGCGGGTAACCCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAg                                                         <  4:233529/149‑1 (MQ=38)
                                                                                                tatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTacc                                              >  7:218160/1‑147 (MQ=255)
                                                                                                tatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGa                                            <  1:221226/149‑1 (MQ=255)
                                                                                                       aGACGGGGGATTGGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCACGATCTGTGCACTGTAGCAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACtt                                     <  3:297015/149‑1 (MQ=255)
                                                                                                       aGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                            >  6:266238/1‑110 (MQ=255)
                                                                                                       aGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                            <  5:266238/110‑1 (MQ=255)
                                                                                                         ctggggtATTTGTACCTGAGGTCCGGGTCTCTTGGGGGTAAACCCCGCTCATC‑‑CAAACATGCATCATTATGTGTGATGTTTTGTTGGGGTTTTCCCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgc                                   <  8:293987/147‑1 (MQ=255)
                                                                                                         aaGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                                             <  7:292777/91‑1 (MQ=255)
                                                                                                         aaGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                                             >  8:292777/1‑91 (MQ=255)
                                                                                                          tgggggATTTGTTACAGAGGTGCCGGGCTATTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgc                                   <  4:221139/147‑1 (MQ=255)
                                                                                                           gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcgc                                 <  7:284472/149‑1 (MQ=255)
                                                                                                           gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                                  >  1:199606/1‑148 (MQ=255)
                                                                                                                     ttACAGAGGGGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGaa                                                                        <  7:147087/100‑1 (MQ=255)
                                                                                                                         agagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGCTCTGTGGACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCgg                                                                          >  7:258977/1‑94 (MQ=38)
                                                                                                                           tggTTCCGGGCTATTCGGGGGTAAACCGGCTCCTA‑‑CACACTTTCCTCATTTTTTGTTCATTTTTGTAGCGGTTTAACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTa                  <  8:226229/147‑1 (MQ=255)
                                                                                                                                   ggCTCTTCGCGGGTAAAACCCCTGCAA‑‑AACACAATCCGCAAGATAGGTGCACGGTAGGAGCGGGGTTTCCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCCCCCTACCGAGGAACTTGCGGAAGAAAAGCCCCGTACGcgagtg          >  4:249212/1‑143 (MQ=255)
                                                                                                                                      tctTCGCGGGTAAACCCGCGCCTA‑‑CACACATGCCTCAAAATATGTGGACTGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGAACAGGTTTTCTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttccttc      >  4:334823/1‑149 (MQ=255)
                                                                                                                                        ttCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                        <  4:109644/113‑1 (MQ=255)
                                                                                                                                        ttCGCGGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                        >  3:109644/1‑113 (MQ=255)
                                                                                                                                           gcgGGTAAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGgcg  <  4:223095/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTA‑‑CACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCG  >  NC_002947/278598‑278883

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