Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,596,231 C→T 7.9% A259T (GCC→ACC)  nemA ← FMN‑linked N‑ethylmaleimide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,596,2310CT 7.9% 222.5 / 5.9 76A259T (GCC→ACC) nemAFMN‑linked N‑ethylmaleimide reductase
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base T (3/3);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

TCTTCGCCTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCG  >  NC_002947/3596083‑3596378
                                                                                                                                                      |                                                                                                                                                    
tCTTCGCCTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGc                                                                                                                                                      <  2:262371/149‑1 (MQ=255)
 cTTCGCCTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGc                                                                                                                                                      >  1:23849/1‑148 (MQ=255)
    cGCCTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCa                                                                                                                                                    <  5:264768/147‑1 (MQ=255)
     gCCTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGtt                                                                                                                                                 <  8:279969/149‑1 (MQ=255)
       cTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTcg                                                                                                                                               >  3:137048/1‑149 (MQ=255)
       cTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTcg                                                                                                                                               >  3:132095/1‑149 (MQ=255)
        tttCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTcgc                                                                                                                                              >  6:12320/1‑149 (MQ=255)
          tctcGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTcgc                                                                                                                                              <  3:269859/147‑1 (MQ=255)
          tctcGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTcgc                                                                                                                                              <  1:204705/147‑1 (MQ=255)
           ctcGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGa                                                                                                                                           >  3:77961/1‑149 (MQ=255)
            tcGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGAt                                                                                                                                          >  1:53416/1‑149 (MQ=255)
              ggcggTATAGGCACCTGCAGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCATTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGTAt                                                                                                                                          <  3:279333/147‑1 (MQ=255)
               gcggTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAg                                                                                                                                       <  2:234853/149‑1 (MQ=255)
               gcggTATAGGCACCGGCCGCGATGATCAACCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAg                                                                                                                                       <  2:247244/149‑1 (MQ=255)
                 ggTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGt                                                                                                                                     <  1:24792/149‑1 (MQ=255)
                   tataGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTAc                                                                                                                                   <  4:94120/149‑1 (MQ=255)
                      aGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAgg                                                                                                                                >  7:105904/1‑149 (MQ=255)
                        gCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGcc                                                                                                                              >  8:39180/1‑149 (MQ=255)
                            cGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCcggcg                                                                                                                          <  3:20087/149‑1 (MQ=255)
                            cGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCcggcg                                                                                                                          >  4:35479/1‑149 (MQ=255)
                                cgcgATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcct                                                                                                                      >  6:107869/1‑149 (MQ=255)
                                cgcgATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcct                                                                                                                      >  5:68504/1‑149 (MQ=255)
                                cgcgATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcct                                                                                                                      >  3:158172/1‑149 (MQ=255)
                                 gcgATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCCGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcctc                                                                                                                     <  4:132095/149‑1 (MQ=255)
                                 gcgATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcctc                                                                                                                     <  4:137048/149‑1 (MQ=255)
                                  cgATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcctcc                                                                                                                    >  3:242314/1‑149 (MQ=255)
                                    atgatCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTg                                                                                                                  <  5:332432/149‑1 (MQ=255)
                                         cACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCtt                                                                                                             >  1:217497/1‑149 (MQ=255)
                                           ccccGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCg                                                                                                           <  3:117642/149‑1 (MQ=255)
                                            cccGGGGTAGGCGGCGCGAATGGCCTGGCGGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGt                                                                                                                                     >  1:207837/1‑123 (MQ=255)
                                            cccGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGt                                                                                                                                     <  2:207837/122‑1 (MQ=255)
                                            cccGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGc                                                                                                          >  4:305689/1‑149 (MQ=255)
                                            cccGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGc                                                                                                          <  4:77961/149‑1 (MQ=255)
                                            cccGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGc                                                                                                          >  2:241536/1‑149 (MQ=255)
                                             ccGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATc                                                                                                                                         <  1:2455/117‑1 (MQ=255)
                                             ccGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCCTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATc                                                                                                                                         >  2:2455/1‑117 (MQ=255)
                                              cGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCa                                                                                                        <  7:236002/149‑1 (MQ=255)
                                                gggTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCAtt                                                                                                      <  8:31147/149‑1 (MQ=255)
                                                       cggcgCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCcaca                                                                                               >  6:70733/1‑149 (MQ=255)
                                                          cgcgAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcc                                                                                                                       >  8:104208/1‑122 (MQ=255)
                                                          cgcgAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGcc                                                                                                                       <  7:104208/122‑1 (MQ=255)
                                                            cgAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCAtt                                                                                          <  5:12320/149‑1 (MQ=255)
                                                               aTGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGaa                                                                                       <  8:105904/149‑1 (MQ=255)
                                                                    cTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAaccac                                                                                  <  2:217497/149‑1 (MQ=255)
                                                                          ggAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAAc                                                                             <  6:68504/149‑1 (MQ=255)
                                                                          ggAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAAc                                                                             >  2:254964/1‑149 (MQ=255)
                                                                          ggAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAAc                                                                             <  4:158172/149‑1 (MQ=255)
                                                                                 gTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGTGAac                                                                       <  5:70733/148‑1 (MQ=255)
                                                                                 gTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAaca                                                                      >  6:310652/1‑149 (MQ=255)
                                                                                 gTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAaca                                                                      >  3:148646/1‑149 (MQ=255)
                                                                                     cGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACAcgcg                                                                  <  2:180218/149‑1 (MQ=255)
                                                                                        aGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGaa                                                                >  7:132147/1‑148 (MQ=255)
                                                                                        aGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAAt                                                              <  5:310652/149‑1 (MQ=255)
                                                                                          tGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGGAGGTTGCGGTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCGTCGTGGTGTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATa                                                             >  1:245703/1‑149 (MQ=255)
                                                                                          tGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCTGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATAc                                                            <  3:35479/149‑1 (MQ=255)
                                                                                            gCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCg                                                          >  3:323308/1‑149 (MQ=255)
                                                                                             cTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGa                                                         <  2:245703/149‑1 (MQ=255)
                                                                                              ttGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGAc                                                        >  3:28172/1‑149 (MQ=255)
                                                                                               tGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGTCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCAtt                                                                                          <  5:191423/115‑1 (MQ=255)
                                                                                               tGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGTCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCAtt                                                                                          >  6:191423/1‑115 (MQ=255)
                                                                                                gCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGaa                                                                        <  8:42852/132‑1 (MQ=255)
                                                                                                  cACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGTCAGTTCGCGGATCAgg                                                                                                                                      <  1:197503/68‑1 (MQ=255)
                                                                                                  cACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGTCAGTTCGCGGATCAgg                                                                                                                                      >  2:197503/1‑68 (MQ=255)
                                                                                                  cACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGTCAGTTCGCGGATCAgg                                                                                                                                      >  2:192265/1‑68 (MQ=255)
                                                                                                  cACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGACGTTGCGCTTGGTCAGTTCGCGGATCAgg                                                                                                                                      <  1:192265/68‑1 (MQ=255)
                                                                                                    ccGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTc                                                  >  1:332134/1‑149 (MQ=255)
                                                                                                     cGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCg                                                 <  6:85394/149‑1 (MQ=255)
                                                                                                                 gCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTg                                     <  4:242314/149‑1 (MQ=255)
                                                                                                                        cAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATgg                              <  7:39180/149‑1 (MQ=255)
                                                                                                                           gTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCg                           >  7:32079/1‑149 (MQ=255)
                                                                                                                           gTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCg                           >  2:132499/1‑149 (MQ=255)
                                                                                                                                 gTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCa                     >  6:336670/1‑149 (MQ=255)
                                                                                                                                 gTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCa                     <  4:148646/149‑1 (MQ=255)
                                                                                                                                    ggCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCg                  >  4:100665/1‑149 (MQ=255)
                                                                                                                                    ggCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCg                  <  3:305689/149‑1 (MQ=255)
                                                                                                                                              gcgcTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAACTACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAg           <  1:254964/148‑1 (MQ=255)
                                                                                                                                                    ggCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAAt     <  5:336670/147‑1 (MQ=255)
                                                                                                                                                     gCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATgcg  >  4:160602/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TCTTCGCCTTTCTCGGCGGTATAGGCACCGGCCGCGATGATCACCCCGGGGTAGGCGGCGCGAATGGCCTGGCGGAACTCGTCACGCAGTGGCTTGCCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCG  >  NC_002947/3596083‑3596378

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