Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 278,744 +A intergenic (‑239/+75) 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.A91.4% 215.1 / 13.0 70intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base A (27/37);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGG  >  NC_002947/278596‑278881
                                                                                                                                                     |                                                                                                                                            
tGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTa                                                                                                                                             <  6:116630/148‑1 (MQ=255)
 gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaa                                                                                                                                           >  3:110827/1‑148 (MQ=255)
 gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                            >  8:368321/1‑148 (MQ=255)
 gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATATATAGAAGGGGGATTAGTTACAGAGGTGACGGGCTCTTCGCGGGTaa                                                                                                                                            >  5:7849/1‑148 (MQ=255)
 gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACACTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                            >  8:190576/1‑148 (MQ=255)
  tCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                          <  5:96179/149‑3 (MQ=255)
  tCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                          <  2:81222/149‑3 (MQ=255)
    gCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacc                                                                                                                                        <  6:127731/148‑1 (MQ=255)
     cTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacccg                                                                                                                                      >  1:363248/1‑147 (MQ=255)
       aCAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCt                                                                                                                                     >  7:21441/1‑149 (MQ=255)
         aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACATAGGTGCCGGGCTCTTCGAGGGTAAAACCGCTcc                                                                                                                                   >  5:7366/1‑149 (MQ=255)
         aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                   <  6:168559/149‑1 (MQ=255)
         aGCTGATCGAGACAGACGTTTTCTGGTCAGGTTGTCAGGTATATATACCCAAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                   <  8:350927/149‑1 (MQ=255)
           cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCgctgct                                                                                                                                  >  7:189008/1‑148 (MQ=255)
              aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                                                               >  5:405361/1‑148 (MQ=255)
               tCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                              <  5:104397/148‑1 (MQ=255)
                 gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacacac                                                                                                                            <  6:219534/148‑1 (MQ=255)
                 gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTAcacaca                                                                                                                           <  2:137464/149‑1 (MQ=255)
                  cgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                          <  2:391229/149‑1 (MQ=255)
                   gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATg                                                                                                                         <  3:128386/149‑1 (MQ=255)
                     cAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAAGTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGCGCCGGGCTATTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                                        >  6:435348/1‑148 (MQ=255)
                      aGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCt                                                                                                                      <  3:175767/149‑1 (MQ=255)
                         cGTGTGCTGGTCCGGTTGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                                   <  2:343514/149‑1 (MQ=255)
                         cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCa                                                                                                                    <  7:394588/148‑1 (MQ=255)
                             tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATc                                                                                                               >  1:334311/1‑149 (MQ=255)
                                  gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                            <  8:291189/147‑1 (MQ=255)
                                  gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGc                                                                                                          >  8:234563/1‑149 (MQ=255)
                                      ggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGATACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                      <  4:99431/149‑1 (MQ=255)
                                       gtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTACTACACACATGCCTCAATATATGTGCACTg                                                                                                      >  5:33965/1‑148 (MQ=255)
                                        tggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                                    <  3:395696/149‑1 (MQ=255)
                                               gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGCGCACTGTAGGAGCgg                                                                                             <  8:325182/149‑1 (MQ=255)
                                                tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                            <  8:346175/149‑1 (MQ=255)
                                                  tatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                                           <  2:363248/148‑1 (MQ=255)
                                                       aCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAacc                                                                                     >  1:411530/1‑149 (MQ=255)
                                                         ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                   <  3:279430/149‑1 (MQ=255)
                                                             aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTACAAACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAACCGCGa                                                                                 >  5:361961/1‑147 (MQ=255)
                                                               gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCGCAAGCTCTGTGGACTGTAGGAGCGGGTTTAACCGCGAAAgg                                                                             >  6:367929/1‑149 (MQ=255)
                                                               gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAAACGCTCCTAAAAACATTCCTAAATATATGTTCACAGTAGGATCGGGTTTAAACGCGAAAgg                                                                             >  5:258706/1‑149 (MQ=255)
                                                                aTTTAAAAAATTTTTATTAACTAACTTTTTGTAATTTTTATAGATGGGGTTTTTGTTACTGCGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACTTGCCTCAAGATCTGTGCACTTTAGGAGCGGTTTTACCCGCGAAAgg                                                                             <  4:186482/148‑1 (MQ=255)
                                                                  ttATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                          >  2:294853/1‑149 (MQ=255)
                                                                  ttATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                          >  6:214333/1‑149 (MQ=255)
                                                                      aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                                  >  5:2890/1‑137 (MQ=255)
                                                                      aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                                  <  6:2890/137‑1 (MQ=255)
                                                                      aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGaa                                                                      >  5:197613/1‑149 (MQ=255)
                                                                      aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGaa                                                                      >  1:82429/1‑149 (MQ=255)
                                                                      aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGaa                                                                      >  1:75760/1‑149 (MQ=255)
                                                                                tGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGGGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACt                                                            <  5:367929/149‑1 (MQ=255)
                                                                                     taactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                             <  4:49507/111‑1 (MQ=255)
                                                                                      aactTTTGCTAATATTAAATGATGGGCGTTTTGTAACACAGGTGCGGGGCTATTCGCGGGTAACCCCGTTCCTAAACAATTCCCTCAAGATGTGTGCACTGTTGTGCGGGGTTTACCCGCTAAAGGGCCGGAACAGGTTTACTCGGAgg                                                      <  8:21441/149‑1 (MQ=255)
                                                                                       cctTTTTGGAATATTAATAGACGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATCCCTCAACCTCTGTCCACTGTAGGCGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGc                                                     <  3:72053/148‑1 (MQ=38)
                                                                                               gAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTcc                                             >  7:216287/1‑149 (MQ=255)
                                                                                               gAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTcc                                             >  8:311781/1‑149 (MQ=255)
                                                                                                    tttATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTTAAACCGCTCCAAAACACATGCCTCAAGATCTTGGCACTGTAAGAGCCGGTTTAACCGCGAAAGGGCCGGAAAAGCTTTACAACGAGGCCCCAGTCCCCAgg                                        >  6:42572/1‑149 (MQ=255)
                                                                                                    tttATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAg                                         <  3:95807/148‑1 (MQ=255)
                                                                                                          gAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCCCACCGAAACCCATGCCTCAAGATCTGTGCACAGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGAACAGGTTTAATAGGAGGCCGCACTCCCGCGGAACATg                                  >  6:280349/1‑149 (MQ=255)
                                                                                                          caaGGGGGATTTGTAACAGCGGTGCCGGGCTCTTCGCTGGAAAACCCGCTCCTACACACATGCCTCAACATCTGTGCACTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                                  <  7:368321/148‑1 (MQ=255)
                                                                                                           ccggggcATTTGTTCCAGGGTTGCGGTGTCTTTCGCGGGTAAACCGGCTCGTACACACATCCACCACGATTTTTGCTCTGTTGGAGCGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGTCCGCAGTCCCGAGGAACTTgc                                 <  6:7849/147‑1 (MQ=255)
                                                                                                           aaGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                                  <  2:82429/148‑1 (MQ=255)
                                                                                                           aaGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                                  <  2:75760/148‑1 (MQ=255)
                                                                                                             gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcgc                               <  1:294853/149‑1 (MQ=255)
                                                                                                             gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                                >  3:333810/1‑148 (MQ=255)
                                                                                                             gggggAGTGGTTACCGAGGTCCCTGTCTCTTCGCGGTTAAACCCGCTCCTACACCCATGCCTCAAGATCTGTGCCCTGTCGGTGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcgc                               <  6:361961/149‑1 (MQ=255)
                                                                                                               gggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTTAAACCGCACCAACACACATGCCACAAGAAATTTGGACAGTAGGAGCGGGGTTACCCGCCAAAGGGGCGGGACACGTTTACTGGGAGGGCGCCGGCCCCAGGAACTTGCGCaa                             >  6:24072/1‑149 (MQ=255)
                                                                                                               gggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCaa                             <  4:110827/149‑1 (MQ=255)
                                                                                                                  aattGTAACAGAGGGGCCGGGCTCTTCGCGCGTAAACCCGCTCCTACACACATCCCTCAACATCTGTGCACTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgc                                 <  8:301734/140‑1 (MQ=38)
                                                                                                                  aTTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTTCCCGCGTAAGGGGCGGATCAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgg                                 >  7:301734/1‑141 (MQ=255)
                                                                                                                       ttACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc                     >  2:44800/1‑149 (MQ=255)
                                                                                                                          cAGAGGTGCCGGGCTCTTCGCGGGTAAAACCCCTCCAAAACACCTGGCGAAAGATATGTGGACGGTAGGAGCGGGGTTACCCGCCCAAGGGGCGGAACAGGGTTTCTAGGAGGCCGCAGTCCCGAGGAACATTCGCAAGAACAGCcgcg                  >  6:22823/1‑149 (MQ=255)
                                                                                                                              ggTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGTACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTAc               <  4:333810/148‑1 (MQ=255)
                                                                                                                                  ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCAGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctc           <  8:216287/148‑1 (MQ=255)
                                                                                                                                  ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctct          <  6:197613/149‑1 (MQ=255)
                                                                                                                                    gggCTCTGCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGTTCTGTGCACTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttc        <  2:46616/149‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttccttc    >  8:322374/1‑149 (MQ=255)
                                                                                                                                         ctTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCg   >  4:205943/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttccttc    <  5:214333/147‑1 (MQ=255)
                                                                                                                                            cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCgg  <  7:322374/147‑1 (MQ=255)
                                                                                                                                            cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCg   <  2:309718/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
TGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGG  >  NC_002947/278596‑278881

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: