Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 679,889 C→A 5.4% S259Y (TCC→TAC)  PP_0582 → thiolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947679,8890CA 5.4% 239.6 / 2.8 74S259Y (TCC→TAC) PP_0582thiolase family protein
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base A (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

CCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGATGG  >  NC_002947/679742‑680035
                                                                                                                                                   |                                                                                                                                                  
ccGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTcc                                                                                                                                                   >  4:31516/1‑149 (MQ=255)
 cGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCa                                                                                                                                                  >  4:365922/1‑149 (MQ=255)
  gAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCAc                                                                                                                                                 <  1:190235/149‑1 (MQ=255)
        tGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGtt                                                                                                                                           <  3:272342/149‑1 (MQ=255)
           cAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGTGGAATTACCCGCCGtcaat                                                                                                                                        <  8:321004/149‑5 (MQ=255)
            aGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGAcc                                                                                                                                       >  6:329971/1‑149 (MQ=255)
              gacgacCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAAACAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGa                                                                                                                                     >  7:282955/1‑149 (MQ=255)
                         gACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCg                                                                                                                           <  3:365922/148‑1 (MQ=255)
                          aCGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTTACGGCGGGTAATTCCACGCCGTTGACCTTTGGTGCTTCGCt                                                                                                                         <  4:193117/149‑1 (MQ=255)
                           cGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTg                                                                                                                        <  4:143646/149‑1 (MQ=255)
                           cGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTg                                                                                                                        <  4:52119/149‑1 (MQ=255)
                            gCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTgg                                                                                                                       >  6:232040/1‑149 (MQ=255)
                               gTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACCGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGc                                                                                                                    >  7:168201/1‑149 (MQ=255)
                                    tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGGTTCGCTGGTGGTGTTg                                                                                                               >  8:186847/1‑149 (MQ=255)
                                    tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCg                                                                                                                           <  1:68306/137‑1 (MQ=255)
                                    tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCg                                                                                                                           >  2:68306/1‑137 (MQ=255)
                                    tttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTg                                                                                                               >  3:111176/1‑149 (MQ=255)
                                     ttcGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTgg                                                                                                              >  8:317951/1‑149 (MQ=255)
                                            gCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAACCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGa                                                                                                       >  6:130096/1‑149 (MQ=255)
                                             cACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAg                                                                                                      <  3:295137/149‑1 (MQ=255)
                                                 ggACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGGTTCGCTGGTGCTGTTGGGCAGCGAGGCGt                                                                                                  >  5:363784/1‑149 (MQ=255)
                                                   acaacaACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTgg                                                                                                <  1:316835/149‑1 (MQ=255)
                                                     aacaacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgc                                                                                              >  3:143156/1‑149 (MQ=255)
                                                     aacaacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgc                                                                                              >  6:249288/1‑149 (MQ=255)
                                                       caacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgcg                                                                                             >  1:390945/1‑148 (MQ=255)
                                                        aacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGTGGAATTACCCGc                                                                                                                                               >  5:70780/1‑97 (MQ=255)
                                                        aacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGTGGAATTACCCGc                                                                                                                                               <  6:70780/97‑1 (MQ=255)
                                                        aacCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGgcg                                                                                             >  8:210854/1‑147 (MQ=255)
                                                          ccTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCgcag                                                                                          >  1:388220/1‑148 (MQ=255)
                                                             gcgcCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCgcagcag                                                                                       >  8:5564/1‑148 (MQ=255)
                                                               gcCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                                    <  2:388220/149‑1 (MQ=255)
                                                               gcCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                                    <  5:399544/149‑1 (MQ=255)
                                                                ccccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAgg                                                                                   <  6:363784/149‑1 (MQ=255)
                                                                  ccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTgg                                                                                                <  8:379654/134‑1 (MQ=255)
                                                                  ccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTgg                                                                                                >  7:379654/1‑134 (MQ=255)
                                                                  ccGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAgggg                                                                                 <  7:317951/149‑1 (MQ=255)
                                                                          aGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAgg                                                                          <  4:143156/148‑1 (MQ=255)
                                                                            tatttaaCACTTGGCCAGTCTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTg                                                                        <  8:168201/143‑1 (MQ=255)
                                                                            tCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTg                                                                                                                     <  3:390394/103‑1 (MQ=255)
                                                                            tCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTg                                                                                                                     >  4:390394/1‑103 (MQ=255)
                                                                            tCTTGAACAGTTGGCCAGACACAACCCAGCCTTAGACCACAGCGGGCAGGCCACGTTGACGCCGGGTAATTCCACGCCGATGACCGATCGTGCTTCGCAGGTCCTGTTCGGCAGCGAGGCGTGCGCGCAGCAGCAGCGGATGCAGCTGc                                                                       >  4:322640/1‑149 (MQ=255)
                                                                                gAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGc                                                                   >  6:349560/1‑149 (MQ=255)
                                                                                      ttGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCt                                                             >  2:328754/1‑149 (MQ=255)
                                                                                         gCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGt                                                          <  1:163263/149‑1 (MQ=255)
                                                                                           cAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTc                                                         >  1:278605/1‑148 (MQ=255)
                                                                                            aGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGa                                                       >  4:14380/1‑149 (MQ=255)
                                                                                             gACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGAt                                                      <  5:329971/149‑1 (MQ=255)
                                                                                                tCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCtt                                                                     >  8:27169/1‑131 (MQ=255)
                                                                                                tCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCtt                                                                     <  7:27169/131‑1 (MQ=255)
                                                                                                  aaaCCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCg                                                  >  1:301170/1‑148 (MQ=255)
                                                                                                   aaCCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTCGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATgg                                                    >  6:392792/1‑145 (MQ=255)
                                                                                                   aaCCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTCGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATgg                                                    <  5:392792/145‑1 (MQ=255)
                                                                                                    aCCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                                    >  4:179296/1‑112 (MQ=255)
                                                                                                    aCCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGcagcagcag                                                                                    <  3:179296/112‑1 (MQ=255)
                                                                                                    aCCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGaaa                                               <  3:14380/149‑1 (MQ=255)
                                                                                                     ccAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAAc                                              >  5:101839/1‑149 (MQ=255)
                                                                                                      cAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAAcc                                             <  5:232040/149‑1 (MQ=255)
                                                                                                      cAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAAcc                                             >  6:310866/1‑149 (MQ=255)
                                                                                                        gCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGc                                           >  4:280960/1‑149 (MQ=255)
                                                                                                        gCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGc                                           <  7:5564/149‑1 (MQ=255)
                                                                                                          cTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTg                                         <  5:130096/149‑1 (MQ=255)
                                                                                                             tGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGc                                        <  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)
                                                                                                                                              tAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGa     <  1:328754/149‑1 (MQ=255)
                                                                                                                                                aTTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGAt    <  2:278605/148‑1 (MQ=255)
                                                                                                                                                 ttCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGATgg  <  2:301170/149‑1 (MQ=255)
                                                                                                                                                 ttCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGATgg  >  1:378700/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGAGGGTTGGCAGGACGACCTGTTGACGCCGTTTCTTTCGCTGGCACGGGACAACAACCTGCGCCCCGACCTGAGTCTTGAACAGTTGGCCAGACTCAAACCAGCCTTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGATGG  >  NC_002947/679742‑680035

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