Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 679,882 G→T 5.6% G257C (GGT→TGT)  PP_0582 → thiolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947679,8820GT 5.6% 220.7 / 2.7 71G257C (GGT→TGT) PP_0582thiolase family protein
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base T (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CCTATGCCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTT  >  NC_002947/679736‑680030
                                                                                                                                                   |                                                                                                                                                    
ccTATGCCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGAcggcgg                                                                                                                                                      <  6:419805/147‑1 (MQ=255)
ccTATGCCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCggg                                                                                                                                                     >  3:143646/1‑148 (MQ=255)
  tATGCCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTaa                                                                                                                                                  >  4:295137/1‑149 (MQ=255)
    tGCCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCCGGGCAGGGCACGTTGACGGCGGGTaa                                                                                                                                                  <  1:218122/148‑1 (MQ=255)
      ccGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTcc                                                                                                                                              >  4:31516/1‑149 (MQ=255)
       cGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCa                                                                                                                                             >  4:365922/1‑149 (MQ=255)
        gAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCAc                                                                                                                                            <  1:190235/149‑1 (MQ=255)
              tGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGtt                                                                                                                                      <  3:272342/149‑1 (MQ=255)
                 cAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGTGGAATTACCCGCCGtcaat                                                                                                                                   <  8:321004/149‑5 (MQ=255)
                  aGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGAcc                                                                                                                                  >  6:329971/1‑149 (MQ=255)
                    gacgacCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAAACAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGa                                                                                                                                >  7:282955/1‑149 (MQ=255)
                               gACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCg                                                                                                                      <  3:365922/148‑1 (MQ=255)
                                aCGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTTACGGCGGGTAATTCCACGCCGTTGACCTTTGGTGCTTCGCt                                                                                                                    <  4:193117/149‑1 (MQ=255)
                                 cGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTg                                                                                                                   <  4:143646/149‑1 (MQ=255)
                                 cGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTg                                                                                                                   <  4:52119/149‑1 (MQ=255)
                                  gCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTgg                                                                                                                  >  6:232040/1‑149 (MQ=255)
                                     gTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACCGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGc                                                                                                               >  7:168201/1‑149 (MQ=255)
                                          tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGGTTCGCTGGTGGTGTTg                                                                                                          >  8:186847/1‑149 (MQ=255)
                                          tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCg                                                                                                                      >  2:68306/1‑137 (MQ=255)
                                          tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCg                                                                                                                      <  1:68306/137‑1 (MQ=255)
                                          tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTg                                                                                                          >  3:111176/1‑149 (MQ=255)
                                           ttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTgg                                                                                                         >  8:317951/1‑149 (MQ=255)
                                                  gCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAACCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGa                                                                                                  >  6:130096/1‑149 (MQ=255)
                                                   cACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAg                                                                                                 <  3:295137/149‑1 (MQ=255)
                                                       ggACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGGTTCGCTGGTGCTGTTGGGCAGCGAGGCGt                                                                                             >  5:363784/1‑149 (MQ=255)
                                                         acaacaACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTgg                                                                                           <  1:316835/149‑1 (MQ=255)
                                                           aacaacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgc                                                                                         >  3:143156/1‑149 (MQ=255)
                                                           aacaacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgc                                                                                         >  6:249288/1‑149 (MQ=255)
                                                             caacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgcg                                                                                        >  1:390945/1‑148 (MQ=255)
                                                              aacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGTGGAATTACCCGc                                                                                                                                          >  5:70780/1‑97 (MQ=255)
                                                              aacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGTGGAATTACCCGc                                                                                                                                          <  6:70780/97‑1 (MQ=255)
                                                              aacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgcg                                                                                        >  8:210854/1‑147 (MQ=255)
                                                                ccTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCgcag                                                                                     >  1:388220/1‑148 (MQ=255)
                                                                   gcgcCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCgcagcag                                                                                  >  8:5564/1‑148 (MQ=255)
                                                                     gcCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                               <  2:388220/149‑1 (MQ=255)
                                                                     gcCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                               <  5:399544/149‑1 (MQ=255)
                                                                      ccccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAgg                                                                              <  6:363784/149‑1 (MQ=255)
                                                                        ccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTgg                                                                                           <  8:379654/134‑1 (MQ=255)
                                                                        ccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTgg                                                                                           >  7:379654/1‑134 (MQ=255)
                                                                        ccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAgggg                                                                            <  7:317951/149‑1 (MQ=255)
                                                                                aGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAgg                                                                     <  4:143156/148‑1 (MQ=255)
                                                                                  tatttaaCACTTGGCCAGTCTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTg                                                                   <  8:168201/143‑1 (MQ=255)
                                                                                  tCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTg                                                                                                                >  4:390394/1‑103 (MQ=255)
                                                                                  tCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTg                                                                                                                <  3:390394/103‑1 (MQ=255)
                                                                                  tCTTGAACAGTTGGCCAGACACAACCCAGCCTTAGACCACAG‑CGGGCAGGCCACGTTGACGCCGGGTAATTCCACGCCGATGACCGATCGTGCTTCGCAGGTCCTGTTCGGCAGCGAGGCGTGCGCGCAGCAGCAGCGGATGCAGCTGc                                                                  >  4:322640/1‑149 (MQ=255)
                                                                                      gAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGc                                                              >  6:349560/1‑149 (MQ=255)
                                                                                            ttGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCt                                                        >  2:328754/1‑149 (MQ=255)
                                                                                               gCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGt                                                     <  1:163263/149‑1 (MQ=255)
                                                                                                 cAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTc                                                    >  1:278605/1‑148 (MQ=255)
                                                                                                  aGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGa                                                  >  4:14380/1‑149 (MQ=255)
                                                                                                   gACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGAt                                                 <  5:329971/149‑1 (MQ=255)
                                                                                                      tCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCtt                                                                >  8:27169/1‑131 (MQ=255)
                                                                                                      tCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCtt                                                                <  7:27169/131‑1 (MQ=255)
                                                                                                        aaaCCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCg                                             >  1:301170/1‑148 (MQ=255)
                                                                                                         aaCCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTCGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATgg                                               <  5:392792/145‑1 (MQ=255)
                                                                                                         aaCCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTCGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATgg                                               >  6:392792/1‑145 (MQ=255)
                                                                                                          aCCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                               <  3:179296/112‑1 (MQ=255)
                                                                                                          aCCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                               >  4:179296/1‑112 (MQ=255)
                                                                                                          aCCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGaaa                                          <  3:14380/149‑1 (MQ=255)
                                                                                                           ccAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAAc                                         >  5:101839/1‑149 (MQ=255)
                                                                                                            cAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAAcc                                        <  5:232040/149‑1 (MQ=255)
                                                                                                            cAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAAcc                                        >  6:310866/1‑149 (MQ=255)
                                                                                                              gCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGc                                      >  4:280960/1‑149 (MQ=255)
                                                                                                              gCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGc                                      <  7:5564/149‑1 (MQ=255)
                                                                                                                cTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTg                                    <  5:130096/149‑1 (MQ=255)
                                                                                                                   tGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGc                                   <  7:210854/147‑1 (MQ=255)
                                                                                                                               ggcagggcaCGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATgg                                               <  8:16607/124‑1 (MQ=255)
                                                                                                                                gcagggcaCGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCa                     >  4:173509/1‑149 (MQ=255)
                                                                                                                                  agggcaCGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACg                   <  4:111176/149‑1 (MQ=255)
                                                                                                                                    ggcaCGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACggg                 >  4:193250/1‑149 (MQ=255)
                                                                                                                                       aagTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGATGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATGTGGTCGATTGCGAAAACGCTGCGGTGGATTTTTTCACGGGgcg              >  3:250469/3‑149 (MQ=255)
                                                                                                                                         gTTGACGGCGGGTAATTCCACGGCGTTGACCGATGGTTCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGcgcg            <  3:322640/149‑1 (MQ=255)
                                                                                                                                           caaCGGCGTGGAATTAACCGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGcgcg            >  1:234872/3‑147 (MQ=255)
                                                                                                                                                gcggGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGgtt     <  3:31516/149‑1 (MQ=255)
                                                                                                                                                   ggTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGgttgtt  <  2:390945/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTATGCCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAG‑CGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTT  >  NC_002947/679736‑680030

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: ‑