Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,737 A→T 7.4% intergenic (‑43/‑209) tktA ← / → sahR transketolase A/methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,7370AT 7.4% 195.5 / 4.1 67intergenic (‑43/‑209)tktA/sahRtransketolase A/methionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base T (3/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TCCGCCATGCCCATGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAAT  >  NC_002947/5657591‑5657885
                                                                                                                                                  |                                                                                                                                                    
tCCGCCATGCCCATGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAgg                                                                                                                                                    >  7:186553/1‑149 (MQ=255)
  cGCCATGCCCATGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAggag                                                                                                                                                  >  2:112981/1‑147 (MQ=255)
            aTGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCg                                                                                                                                        >  3:185895/1‑149 (MQ=255)
             tGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAggcggc                                                                                                                                                >  3:56902/1‑140 (MQ=255)
             tGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAggcggc                                                                                                                                                <  4:56902/140‑1 (MQ=255)
                    aCCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCCTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTgg                                                                                                                                <  5:118475/149‑1 (MQ=255)
                      cTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTgggg                                                                                                                              >  6:247026/1‑149 (MQ=255)
                         ggTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGgcg                                                                                                                           >  8:198462/1‑149 (MQ=255)
                             gCCGCTGTTGGCCTGTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGGGGTGTCACGAAAAAGGCGGCCATTTTCGCCCCCAGTGGgcgaaaa                                                                                                                       <  7:275306/149‑5 (MQ=255)
                               cGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTggg                                                                                                                     >  6:183643/1‑149 (MQ=255)
                               cGCTGTTGGCCTTTTGCACGGCATCAATGCTGGGGGCACGAATGGCGTTGGCCGGTTGACGGCGGCTGGGCATGGCTGAGCTCCTAGGGGGGTTGAATAGGTGGTGTCACGAAAATGGCGGCCATGTTCGCCCACTGGGGGCGCTTggg                                                                                                                     >  6:338867/1‑149 (MQ=255)
                                gCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCCCAGTGGgcgaaaatgg                                                                                                                    <  8:186553/149‑8 (MQ=255)
                                    ttGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCaaa                                                                                                                <  6:200477/149‑1 (MQ=255)
                                      ggCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACa                                                                                                              >  3:274318/1‑149 (MQ=255)
                                          ttttGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACAttt                                                                                                           >  6:49242/1‑148 (MQ=255)
                                            ttGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaa                                                                                                         <  5:324801/148‑1 (MQ=255)
                                               cACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgc                                                                                                     >  3:125924/1‑149 (MQ=255)
                                                  ggCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgccaa                                                                                                  <  7:198462/149‑1 (MQ=255)
                                                   gCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgccaag                                                                                                 >  3:178702/1‑149 (MQ=255)
                                                     aTCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgccaagcg                                                                                               >  5:133800/1‑149 (MQ=255)
                                                      tCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGCGCTTAAATGTttgccccaagcgc                                                                                              <  4:185895/149‑1 (MQ=255)
                                                        cATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCg                                                                                            <  4:274318/149‑1 (MQ=255)
                                                         aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                     >  4:29922/1‑104 (MQ=255)
                                                         aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                     <  3:29922/107‑4 (MQ=255)
                                                         aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                     <  3:25307/107‑4 (MQ=255)
                                                         aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                     >  4:25307/1‑104 (MQ=255)
                                                         aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCgg                                                                                           >  6:170021/1‑149 (MQ=255)
                                                            cTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCtt                                                                                        <  5:170021/149‑1 (MQ=255)
                                                               aGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGGGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCttt                                                                                       <  6:319378/147‑1 (MQ=255)
                                                               aGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCttt                                                                                       <  5:247026/147‑1 (MQ=255)
                                                               aGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTaa                                                                                     <  5:49242/149‑1 (MQ=255)
                                                                 ggCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGc                                                                                   <  4:125924/149‑1 (MQ=255)
                                                                 ggCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGc                                                                                   >  4:499/1‑149 (MQ=255)
                                                                   cACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTc                                                                                                                                         >  2:354813/1‑93 (MQ=255)
                                                                   cACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTc                                                                                                                                         <  1:354813/93‑1 (MQ=255)
                                                                           gCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTga                                                                          >  3:5670/1‑148 (MQ=255)
                                                                           gCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTga                                                                          <  4:5670/148‑1 (MQ=255)
                                                                            cGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTg                                                                           >  7:91604/1‑146 (MQ=255)
                                                                             gTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGCCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTgaaga                                                                       <  1:112981/149‑1 (MQ=255)
                                                                                   aCGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCa                                                                   >  4:198716/1‑147 (MQ=255)
                                                                                    cGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAACCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCagag                                                                >  7:98304/1‑149 (MQ=255)
                                                                                                          cTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAg                                          >  5:311841/1‑149 (MQ=255)
                                                                                                            gATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCaaa                                           <  4:260509/146‑1 (MQ=255)
                                                                                                            gATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCaaa                                           >  3:260509/1‑146 (MQ=255)
                                                                                                             aTCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGc                                         <  3:499/147‑1 (MQ=255)
                                                                                                              tctcCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCgtg                                      >  3:180307/1‑149 (MQ=255)
                                                                                                                  atAGGGTGCTTGAATAGGTGGTGTCAAGAAAAAGGCGGCCATTTTATCCCACTGGGGGCGATTGTGGCAAACATTTAATCGAAAAGCGCCGGCTTTAAGACGAAACTGAAGAAGCAGAGCGAGACTCGGGGAGACCCAAAGCgtgaga                                   >  5:74055/2‑145 (MQ=255)
                                                                                                                   tAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACt                                 >  3:20315/1‑149 (MQ=255)
                                                                                                                     gggggCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAACCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACt                                 <  8:98304/147‑1 (MQ=255)
                                                                                                                     gggggCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACtt                                >  8:201126/1‑148 (MQ=255)
                                                                                                                           ttGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTaa                                                                                     <  4:96686/89‑1 (MQ=255)
                                                                                                                           ttGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTaa                                                                                     >  3:96686/1‑89 (MQ=255)
                                                                                                                               aTAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCgg                     >  3:322603/1‑149 (MQ=255)
                                                                                                                               aTAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCg                      <  5:183643/148‑1 (MQ=255)
                                                                                                                                tAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGc                    <  4:322603/149‑1 (MQ=255)
                                                                                                                                   gtggtgTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTg                 >  4:346989/1‑149 (MQ=255)
                                                                                                                                   gtggtgTCACGAAAAAGGAGGCCATTTTCGCACACTGGGGGCGATTGGGGCAAACATTTAAGCGCCAAACGACGGATTTAAGCAGAAACTGAAGAAGCAGAGCGCGACACGGGGCGACCCAAACCGAGAAACATGAGGCATTCGGCtt                  >  2:139282/1‑148 (MQ=255)
                                                                                                                                    tggtgTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGc                <  4:20315/149‑1 (MQ=255)
                                                                                                                                      gtgtCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCATAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCa               <  1:139282/148‑1 (MQ=255)
                                                                                                                                       tgtCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGc             <  6:133800/149‑1 (MQ=255)
                                                                                                                                       tgtCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGc             <  3:346989/149‑1 (MQ=255)
                                                                                                                                         tCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCTGCTGTAAGCCGCAACTGAAGAAGCAGAGCGCGCCTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCt            >  3:80695/1‑148 (MQ=255)
                                                                                                                                           aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAg              >  5:42111/1‑144 (MQ=255)
                                                                                                                                           aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAg              <  6:42111/144‑1 (MQ=255)
                                                                                                                                           aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCAt         >  8:169436/1‑149 (MQ=255)
                                                                                                                                           aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCAt         >  7:163163/1‑149 (MQ=255)
                                                                                                                                           aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCAt         >  1:163183/1‑149 (MQ=255)
                                                                                                                                            acgaaaaGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATc        >  5:231978/4‑149 (MQ=255)
                                                                                                                                                  aGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCACATCGCCAat  <  4:180307/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCCGCCATGCCCATGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAAT  >  NC_002947/5657591‑5657885

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: