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.A92.9% 217.3 / 9.6 70intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base A (33/32);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTC  >  NC_002947/278595‑278875
                                                                                                                                                       |                                                                                                                                        
aTGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGt                                                                                                                                           >  3:281175/1‑148 (MQ=255)
aTGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGt                                                                                                                                           >  2:162271/1‑148 (MQ=255)
    cGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                      <  4:281175/148‑3 (MQ=255)
         cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTc                                                                                                                               <  1:162271/149‑1 (MQ=255)
          aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTcc                                                                                                                              <  6:242194/149‑1 (MQ=255)
          aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTcc                                                                                                                              >  1:72208/1‑149 (MQ=255)
          aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTACGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTcc                                                                                                                              >  1:65894/1‑149 (MQ=255)
            cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                            <  8:140709/149‑1 (MQ=255)
            cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                            >  6:27048/1‑149 (MQ=255)
             tGATCGCGACAGACGTGTGCTGGTCCGGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCt                                                                                                                             >  1:228744/1‑148 (MQ=255)
               aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                         <  1:326713/149‑1 (MQ=255)
               aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                                                          >  6:62959/1‑148 (MQ=255)
                tCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacaca                                                                                                                        <  2:145611/149‑1 (MQ=255)
                  gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTAcacaca                                                                                                                      >  7:199622/1‑149 (MQ=255)
                  gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAaccc                                                                                                                                  >  2:379358/1‑135 (MQ=255)
                  gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAaccc                                                                                                                                  <  1:379358/136‑2 (MQ=255)
                     aCAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                                   <  5:62959/149‑1 (MQ=255)
                      gagtcGTGTGCTTGTCCGGTGGTCCGGTATATATAGCCTTATGATTTATAAAATTTTTTTGAACTAAGTTTTTGGAATTTTTATAGAAGGGTGATTTGTTACAGAGGTGCCTGGCTGTTCGCGGGTAAACCCGCTCGTACAGACATGc                                                                                                                   <  2:956/144‑1 (MQ=255)
                         aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTc                                                                                                                >  5:384271/1‑148 (MQ=255)
                          cGTGTGCTGGTCCTGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCa                                                                                                               <  2:220505/148‑1 (MQ=255)
                           gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAg                                                                                                             <  5:27048/149‑1 (MQ=255)
                               gCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                                   <  3:425968/139‑1 (MQ=255)
                               gCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                                   >  4:425968/1‑139 (MQ=255)
                                  ggtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTACTACACACATGCCTCAAGATCtgt                                                                                                       >  5:190177/1‑148 (MQ=255)
                                   gtccggtggtcctgtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                       <  2:119829/147‑1 (MQ=255)
                                   gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgtg                                                                                                      <  6:384271/148‑1 (MQ=255)
                                     ccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTAAACACATGCCTCAAGATCTGTGCa                                                                                                    >  4:213135/1‑148 (MQ=255)
                                         gtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                >  6:29749/1‑148 (MQ=255)
                                           ggtccggtATATATACCCTAATGATTTATAAAATTTTTATGTACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                             >  2:360008/1‑149 (MQ=255)
                                            gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                             <  4:56010/148‑1 (MQ=255)
                                            gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                             <  6:61148/148‑1 (MQ=255)
                                            gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                             <  1:196737/148‑1 (MQ=255)
                                             tccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                             >  6:414826/1‑147 (MQ=255)
                                               cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTGCTACACACATGCCTCAAGATCTGTGGACTGTAGGAGc                                                                                         >  4:294496/1‑149 (MQ=255)
                                               cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGc                                                                                         <  5:29749/149‑1 (MQ=255)
                                                     atataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGcc                                                                                                                  <  1:377150/119‑1 (MQ=255)
                                                     atataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGcc                                                                                                                  >  2:377150/1‑119 (MQ=255)
                                                       ataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                       <  1:39605/143‑1 (MQ=255)
                                                       ataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                       >  2:39605/1‑143 (MQ=255)
                                                          cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAccc                                                                               <  8:199622/148‑1 (MQ=255)
                                                            cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACAGCTACTTCACAAATTCCTAAAGATCTTTGGAATGTATTAGCGGGTTTTCCcgc                                                                             >  5:406096/1‑148 (MQ=255)
                                                             tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCa                                                                                                               <  6:193156/114‑1 (MQ=255)
                                                             tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCa                                                                                                               >  5:193156/1‑114 (MQ=255)
                                                             tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTAAAGAGGTGCCGGGATCTTCGCCGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGACAGGGTTTACCcgcg                                                                            >  4:327069/1‑148 (MQ=255)
                                                                      ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                   <  5:414826/148‑1 (MQ=255)
                                                                       taAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCAGCTACTACACACATGACACAAGTTCTGTGAAAGTTAGGAGAGGGTTTAACCGCGAAAGGTCCgg                                                                  >  5:12233/1‑148 (MQ=255)
                                                                         aaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAg                                                                                          <  7:271592/122‑1 (MQ=255)
                                                                         aaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAg                                                                                          >  8:271592/1‑122 (MQ=255)
                                                                         aaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAAc                                                               <  2:72208/149‑1 (MQ=255)
                                                                         aaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCAACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAAc                                                               <  2:65894/149‑1 (MQ=255)
                                                                                  ttaactaactTTTTGTAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCCGCTCCTACACACATGCCTCAATATCTGTCCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACt                                                      <  4:231319/147‑1 (MQ=39)
                                                                                       taactTTTTGGAATATTTATAGAAGGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACa                                                              >  1:5902/1‑137 (MQ=255)
                                                                                         actTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                     >  7:283253/1‑127 (MQ=255)
                                                                                          ctTTTTGGAATATTTATACAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTCTCATGGGTAACCTCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGc                                               <  1:360008/149‑1 (MQ=38)
                                                                                            ttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTGCTACAAAAATGCCTCAAGATGTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGGCGGAACAGGTTTTCTCTGAGGcag                                             >  3:404133/1‑146 (MQ=255)
                                                                                               tGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGAGGGTAAAACCGCTCCTACACACATGACTCAAGATCTGTGCAATGTAGGAGCGGGTTTACCCGCGAAACGGCCGGAACAGGTTTACTCGGAGGCCGCAGt                                         >  1:322917/1‑149 (MQ=255)
                                                                                                ggAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATGTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTc                                        >  6:21949/1‑149 (MQ=255)
                                                                                                 gAATATTTATATAATGGGGATTTGTTTCATTGGTGACGGGATCTTCTTGTTTAAACCCCCTCATAAAAAAATGACTAAAGATCTTTTAACTGTATGAGCGGTTTTAACCGCGAAAGGGACGTAACAGGTTTACTCGGAGGCCGCAGTc                                        <  6:233451/148‑1 (MQ=255)
                                                                                                 gAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTcc                                       >  7:12176/1‑149 (MQ=255)
                                                                                                      tttATTGATGGGGTATTTTTTACAGCGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAATTTCTGTGCACTGTTGGAGGGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAg                                   <  2:228744/148‑1 (MQ=37)
                                                                                                         ataGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGACCGGGTTTACCCGCGAAAGGGCCGGACGAGGTTTACTCGGAGGCCGCAGTCCCGAGGAAc                               >  7:233524/1‑149 (MQ=255)
                                                                                                           aGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                              <  3:131413/84‑1 (MQ=255)
                                                                                                           aGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                              >  4:131413/1‑84 (MQ=255)
                                                                                                            caaGGGGTATTTTTTCCAGAGGTGCCGGGCTCTGCGCGGGAAACCCCGCTCCAACACACATGCCTCAAGATCTGTCCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                            <  3:213135/148‑1 (MQ=38)
                                                                                                                   gtaTTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACAAACATGCCTCAATATATGTGCACTGTAGGTGTGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAg                      <  4:404133/147‑1 (MQ=255)
                                                                                                                   ggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcgc                         >  1:72366/1‑146 (MQ=255)
                                                                                                                               tcaaagacGGGCTCTTCGCGCGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGt           <  4:67318/141‑1 (MQ=39)
                                                                                                                                  gTGCGGGGCTCTTCGCGGGTAAACCCGCTCCAACAAACATGCCTCAAGATCTGTGCACTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACg        <  6:406096/148‑1 (MQ=255)
                                                                                                                                  gTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACGGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACg        >  1:125424/1‑149 (MQ=39)
                                                                                                                                  gTGACGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGAGGGTTTAACCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCCAGGAACTTGCGCAACAACTGCCGCGTACGc       >  2:294805/1‑149 (MQ=255)
                                                                                                                                       gggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttc  <  8:12176/149‑1 (MQ=255)
                                                                                                                                        ggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGCGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                  <  1:375750/100‑1 (MQ=255)
                                                                                                                                        ggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                  >  2:375750/1‑100 (MQ=255)
                                                                                                                                                       |                                                                                                                                        
ATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTC  >  NC_002947/278595‑278875

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