Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,925,396 A→G 6.6% L250P (CTG→CCG) ‡ gcvT‑II ← aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,925,3960AG 6.6% 270.2 / 6.4 90L250P (CTG→CCG) ‡gcvT‑IIaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base G (2/4);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

ACGACCTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGC  >  NC_002947/5925252‑5925538
                                                                                                                                                  |                                                                                                                                               
acgACCTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCCATGCTCCAGCACAGGtt                                                                                                                                             >  3:101322/1‑149 (MQ=255)
    ccTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAg                                                                                                                                        <  1:332995/149‑1 (MQ=255)
      tGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGt                                                                                                                                      >  7:4793/1‑149 (MQ=255)
       gATGGGCGCGCAGTACAACGCGCTCTTCCAGCACCAGACACACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCCACACTGCCGAAGAAGTTGCGCTCGGCCGGAT‑‑CCCAGGCGATGCTGCTGCCCAGGTTTGAGGt                                                                                                                                      >  4:238401/1‑148 (MQ=255)
        aTGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCTGGAgcatc                                                                                                                                    <  3:221020/149‑5 (MQ=255)
        aTGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCTGGAgcatc                                                                                                                                    <  2:54198/149‑5 (MQ=255)
           ggCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCTGGAgcatcgc                                                                                                                                  <  1:2463/148‑7 (MQ=255)
            gcgcgcAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCgg                                                                                                                                <  1:181446/149‑1 (MQ=255)
                 aaGTACACCGAGCTCTTCCAGCACGAGACCCACTAGCTTTTCCTGCACGCCGTGATCGATTTCAGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTGGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGgtg                                                                                                                            <  4:101322/147‑1 (MQ=255)
                  aGTACACCGCGCTCTTCCAGTACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGtgtg                                                                                                                          <  2:159241/149‑1 (MQ=255)
                  aGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGtgtg                                                                                                                          <  5:246946/149‑1 (MQ=255)
                    tACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGt                                                                                                                        >  6:326643/1‑149 (MQ=255)
                       acCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGtttt                                                                                                                     >  4:333845/1‑149 (MQ=255)
                       acCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGtttt                                                                                                                     <  5:326643/149‑1 (MQ=255)
                        ccGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTtc                                                                                                                    >  7:263292/1‑149 (MQ=255)
                        ccGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTtc                                                                                                                    >  1:96895/1‑149 (MQ=255)
                         cgcgCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTtcg                                                                                                                   >  3:172391/1‑149 (MQ=255)
                            gcTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTtcgtcg                                                                                                                <  3:18082/149‑1 (MQ=255)
                              tctTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              >  6:299348/1‑149 (MQ=255)
                                ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              <  2:273367/147‑1 (MQ=255)
                                ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              <  3:313805/147‑1 (MQ=255)
                                ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              >  1:273367/1‑147 (MQ=255)
                                ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGt                                                                                                            >  5:263229/1‑149 (MQ=255)
                                    aGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTg                                                                                                        <  4:172391/149‑1 (MQ=255)
                                    aGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCTGGAgcatcgc                                                                                                                                  >  3:318915/1‑117 (MQ=255)
                                    aGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCTGGAgcatcgc                                                                                                                                  <  4:318915/123‑7 (MQ=255)
                                     gCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTgg                                                                                                       <  2:96895/149‑1 (MQ=255)
                                       aCCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCACGCAGGCCGACACGGCCGATGAAGTTGCGCTCGCCCGGTT‑‑CCCAGGCGATGCTCCAGCCGAGGTTTGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGcc                                                                                                     >  2:316302/1‑149 (MQ=255)
                                         cAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGATGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGATCGGGGTGTGGTTTTCGTCGATGTCCTGGcc                                                                                                     <  6:200357/147‑1 (MQ=255)
                                         cAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGt                                                                                                   <  4:238932/149‑1 (MQ=255)
                                            aCCCACCAGCTTTTCCTGCACGCCGTGTTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGGGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTa                                                                                                  >  8:151705/1‑147 (MQ=255)
                                              ccaccaGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAgg                                                                                              >  3:115607/1‑149 (MQ=255)
                                               caccaGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGt                                                                                             >  6:104704/1‑149 (MQ=255)
                                                    gCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGTTGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATg                                                                                        <  3:11683/149‑1 (MQ=255)
                                                          ccTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCt                                                                                  >  2:118604/1‑149 (MQ=255)
                                                           cTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCtt                                                                                 <  8:263292/149‑1 (MQ=255)
                                                           cTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCt                                                                                  <  1:146470/148‑1 (MQ=255)
                                                             gCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTc                                                                                >  6:357720/1‑148 (MQ=255)
                                                               aCGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAg                                                                             <  8:78329/149‑1 (MQ=255)
                                                                cGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAgg                                                                            <  5:357720/149‑1 (MQ=255)
                                                                 gCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGgc                                                                           >  6:157085/1‑149 (MQ=255)
                                                                     tGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAg                                                                       <  7:350851/149‑1 (MQ=255)
                                                                     tGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAg                                                                       <  3:273915/149‑1 (MQ=255)
                                                                     tGCTCGATTTCCGAAAGCAGGCAGAAACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CACAGGCGATGCTCCAGCCCAGGTTGGAGGTGCGCGGGGTGTGGTTTTCGACGATGTCATGGCAGTACAGGTTCATGCCGGCTTCCAGGCGCAg                                                                       >  6:95624/1‑149 (MQ=255)
                                                                        tCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCTCGCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGt                                                                    <  3:200184/149‑1 (MQ=255)
                                                                         cGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACGGGTTCATGCCGGCTTCCAGGCGCAGGGTa                                                                   >  5:333701/1‑149 (MQ=255)
                                                                         cGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGGTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTa                                                                   >  7:343938/1‑149 (MQ=255)
                                                                          gATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTAt                                                                  >  3:32785/1‑149 (MQ=255)
                                                                                cgaaatcAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCTTGCCGGCTTCCAGGCGCAGGGTATCGCggg                                                            <  3:238401/143‑1 (MQ=255)
                                                                                cGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCgg                                                             <  1:118604/148‑1 (MQ=255)
                                                                                 gCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgc                                                           >  7:150395/1‑149 (MQ=255)
                                                                                 gCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgc                                                           >  4:144924/1‑149 (MQ=255)
                                                                                 gCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgc                                                           >  4:304204/1‑149 (MQ=255)
                                                                                   aaGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgcgc                                                         >  5:76734/1‑149 (MQ=255)
                                                                                         gCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTAGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAgg                                                    >  1:293465/1‑149 (MQ=255)
                                                                                         gCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGc                                                   <  8:4793/149‑1 (MQ=255)
                                                                                              cacGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGgcgc                                                                         >  3:312001/1‑122 (MQ=255)
                                                                                              cacGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGgcgc                                                                         <  4:312001/122‑1 (MQ=255)
                                                                                               acGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCTGCTTg                                             <  5:299348/149‑1 (MQ=255)
                                                                                               acGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTg                                             >  6:200107/1‑149 (MQ=255)
                                                                                                cGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTgg                                            >  4:140353/1‑149 (MQ=255)
                                                                                                 ggCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTggg                                           >  8:133050/1‑149 (MQ=255)
                                                                                                 ggCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGGTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTggg                                           <  3:52753/149‑1 (MQ=255)
                                                                                                  gCCGATGAAGGAGCGCTCGGCCGGTT‑‑CCCAGGCGCTGCTCCAGCCCAGGGTGGAGGTGAGCGGGGTGGGGGTTTCGTCGATGTCCTGGCCGTACAGGTTCCTGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCGTgggg                                          <  4:297446/149‑1 (MQ=255)
                                                                                                     gATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTgcggcg                                        >  8:50653/1‑148 (MQ=255)
                                                                                                         atgTTGCGCTGGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGTGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATAc                                    <  5:157085/146‑1 (MQ=255)
                                                                                                          aGTTGCGCTCGGCCGGTA‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCg                                  >  7:171731/1‑149 (MQ=255)
                                                                                                            ttGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCTGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGgc                                <  1:230581/149‑1 (MQ=255)
                                                                                                               cgcTCGGCCGGCTCCCCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTATCGTCGATCTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGgc                                <  1:316302/148‑1 (MQ=255)
                                                                                                                  tCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGCTACCGGCGCCGAc                          >  5:287870/1‑149 (MQ=255)
                                                                                                                      ccGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAgg                      <  5:200107/149‑1 (MQ=255)
                                                                                                                      ccGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAgg                      >  2:126982/1‑149 (MQ=255)
                                                                                                                       cGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTgg                                            <  3:233764/126‑1 (MQ=255)
                                                                                                                       cGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGt                     >  4:326145/1‑149 (MQ=255)
                                                                                                                         gTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCg                   >  5:163987/1‑149 (MQ=255)
                                                                                                                              cccAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGtt                 <  6:163987/148‑1 (MQ=255)
                                                                                                                              cccAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGtt                 <  6:287870/148‑1 (MQ=255)
                                                                                                                                cAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTg                >  7:189660/1‑147 (MQ=255)
                                                                                                                                 aGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTgaag             <  3:144924/149‑1 (MQ=255)
                                                                                                                                   gctcaccTCCAACCTGGGCTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTgaag             >  3:102617/7‑147 (MQ=255)
                                                                                                                                   gCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGTTGTCCTGGGCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTTTCGCGGGCGCCGAGGCCGGTTGGGGCGATACCGGCGCCGACCAGGTCGTTgaagaa           <  7:50653/149‑1 (MQ=255)
                                                                                                                                    cGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGaagaag          <  2:293465/149‑1 (MQ=255)
                                                                                                                                      aTGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGc        <  7:133050/149‑1 (MQ=255)
                                                                                                                                        gCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCa      >  5:213992/1‑149 (MQ=255)
                                                                                                                                          tCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGCCCAGGTCGTTGAAGAAGGCCAcc    >  4:238018/1‑149 (MQ=255)
                                                                                                                                          tCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCAcc    >  6:341889/1‑149 (MQ=255)
                                                                                                                                            cAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  >  3:283833/1‑149 (MQ=255)
                                                                                                                                            cAGCCCAGGTGGGCGGGGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  <  6:263229/149‑1 (MQ=255)
                                                                                                                                                 cAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  >  7:50411/1‑144 (MQ=255)
                                                                                                                                                 cAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  <  8:50411/144‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACGACCTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGC  >  NC_002947/5925252‑5925538

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