Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 285,098 G→T 5.4% P118T (CCT→ACT)  atsK ← alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947285,0980GT 5.4% 239.7 / 2.7 74P118T (CCT→ACT) atsKalpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base T (2/2);  total (39/35)
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.57e-01

GCGGCGTGGTGCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGG  >  NC_002947/284951‑285246
                                                                                                                                                    |                                                                                                                                                    
gcggcgTGGTGCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAgg                                                                                                                                                      >  3:170772/1‑148 (MQ=255)
gcggcgTGGTGCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAgg                                                                                                                                                      <  6:228643/148‑1 (MQ=255)
   gcgTGGTGCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAAt                                                                                                                                                  >  8:36625/1‑149 (MQ=255)
        gtgCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCtt                                                                                                                                             <  4:115930/149‑1 (MQ=255)
         tgCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTg                                                                                                                                            <  8:83883/149‑1 (MQ=255)
          gCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCCGCGGCGCCGACAAGGCCTCGCAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTgg                                                                                                                                           >  4:240822/1‑149 (MQ=255)
             aaaGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGgcgc                                                                                                                                        >  5:323061/1‑149 (MQ=255)
             aaaGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGgcgc                                                                                                                                        >  1:146861/1‑149 (MQ=255)
                   cTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTa                                                                                                                                  >  7:288114/1‑149 (MQ=255)
                                     gttGAAGTCGTGGGTAGCGGTCAGCCCGTAGAGCATTTCAGGCAGCGGCGCCGACTAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGGTGAGTACCGCGCCCAGTGCCggg                                                                                                                 <  6:130478/146‑1 (MQ=255)
                                      gTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGt                                                                                                               >  7:313564/1‑149 (MQ=255)
                                         aaGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCgg                                                                                                            >  6:196912/1‑149 (MQ=255)
                                          aGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTTgg                                                                                                            >  2:190388/1‑148 (MQ=255)
                                            tCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCa                                                                                                          >  2:144401/1‑148 (MQ=255)
                                                 ggTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAgg                                                                                                    >  7:213387/1‑149 (MQ=255)
                                                 ggTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAgg                                                                                                    <  4:170772/149‑1 (MQ=255)
                                                 ggTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAg                                                                                                     <  5:240209/148‑1 (MQ=255)
                                                 ggTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAg                                                                                                     <  1:144401/148‑1 (MQ=255)
                                                  gTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGt                                                                                                               >  4:8515/1‑137 (MQ=255)
                                                  gTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGt                                                                                                               <  3:8515/137‑1 (MQ=255)
                                                    aGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCa                                                                                                 >  2:92503/1‑149 (MQ=255)
                                                       ggTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGt                                                                                              >  5:74113/1‑149 (MQ=255)
                                                        gTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTc                                                                                             >  6:233532/1‑149 (MQ=255)
                                                           aGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCgg                                                                                           <  8:313564/148‑1 (MQ=255)
                                                             cccGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACcg                                                                                                                               <  2:304003/110‑1 (MQ=255)
                                                             cccGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACcg                                                                                                                               >  1:304003/1‑110 (MQ=255)
                                                               cGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGc                                                                                      >  2:321950/1‑149 (MQ=255)
                                                                   cAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCcaca                                                                                  >  8:192339/1‑149 (MQ=255)
                                                                    aGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCacac                                                                                 >  1:136837/1‑149 (MQ=255)
                                                                    aGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCacac                                                                                 <  2:339518/149‑1 (MQ=255)
                                                                     gCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACg                                                                                >  2:329664/1‑149 (MQ=255)
                                                                      cATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACgg                                                                               <  7:192339/149‑1 (MQ=255)
                                                                        tttCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCacac                                                                                 <  4:305849/145‑1 (MQ=255)
                                                                        tttCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCacac                                                                                 >  3:305849/1‑145 (MQ=255)
                                                                        tttCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCg                                                                             <  7:334091/149‑1 (MQ=255)
                                                                            aCGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAg                                                                                                     <  4:169716/121‑1 (MQ=255)
                                                                            aCGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAg                                                                                                     >  3:169716/1‑121 (MQ=255)
                                                                            aCGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCCCGTCGGTGTGCCACACGGCGTTGt                                                                         >  3:154788/1‑149 (MQ=255)
                                                                            aCGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGt                                                                         >  2:342502/1‑149 (MQ=255)
                                                                              gCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAgg                                                                                                    <  1:172768/120‑1 (MQ=255)
                                                                              gCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAgg                                                                                                    >  2:172768/1‑120 (MQ=255)
                                                                               cAGCGGCGCCGACAAGGCCTCGAAGGCAGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTcgc                                                                      <  8:288114/149‑1 (MQ=255)
                                                                                             aGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAAcc                                                        <  1:92503/149‑1 (MQ=255)
                                                                                               gCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGTGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAAccg                                                       >  1:49238/1‑149 (MQ=255)
                                                                                               gCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAACGTCACGTCGGTGTTCCACACCGCCTTGTCGCGCACAACAGTAAccgc                                                      >  8:277782/1‑149 (MQ=255)
                                                                                                ccTCGAAGGCCGCAATGCCTCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGTGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAAccgcc                                                     <  4:236583/149‑1 (MQ=255)
                                                                                                 cTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAAccgccg                                                    <  6:323061/149‑1 (MQ=255)
                                                                                                  tCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCgt                                                   >  6:2130/1‑149 (MQ=255)
                                                                                                  tCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCgt                                                   >  2:228955/1‑149 (MQ=255)
                                                                                                   cGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCgtg                                                  >  7:244173/1‑149 (MQ=255)
                                                                                                      aGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTc                                                >  6:338952/1‑148 (MQ=255)
                                                                                                        gCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAg                                             <  1:228955/149‑1 (MQ=255)
                                                                                                          cGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGca                                           <  1:34979/149‑1 (MQ=255)
                                                                                                           gCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTTGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGAACATCGGTAACCGCCGTGTCGAGcac                                          <  1:190388/149‑1 (MQ=255)
                                                                                                            cAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGcac                                          >  5:271133/1‑148 (MQ=255)
                                                                                                             aaTGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACaa                                        <  2:136837/149‑1 (MQ=255)
                                                                                                              aTGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCTGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAg                                       <  5:196912/149‑1 (MQ=255)
                                                                                                              aTGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCGTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAg                                       <  8:244173/149‑1 (MQ=255)
                                                                                                               tGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAg                                       >  4:20557/1‑148 (MQ=255)
                                                                                                                  cGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCAcc                                   >  7:302234/1‑149 (MQ=255)
                                                                                                                       gCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCtgtg                               <  1:342502/148‑1 (MQ=255)
                                                                                                                       gCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCtgt                                <  3:20557/147‑1 (MQ=255)
                                                                                                                        cccAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTggg                             <  6:276774/149‑1 (MQ=255)
                                                                                                                        cccAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTggg                             <  1:321950/149‑1 (MQ=255)
                                                                                                                         ccAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTgggg                            >  8:255855/1‑149 (MQ=255)
                                                                                                                           aaaGGTTATGGCCAGCGTAAGCAGGCTATTGCTTGGAGCTGATTACCGCGCCCAGTGCCGGGGTCGGCAAAAATGTCACGTGGGTGTGCCACACGGCGTTGTCGCGCAGATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGgtg                          <  7:277782/149‑1 (MQ=255)
                                                                                                                             aGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTc                        <  1:329664/149‑1 (MQ=255)
                                                                                                                             aGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTc                        >  6:146560/1‑149 (MQ=255)
                                                                                                                               ggTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTg                      <  5:146560/149‑1 (MQ=255)
                                                                                                                                   aTCCCCACCGTAAGCAGGCAATTGCTTGGCGCAGACTACCGCGCCCAGAGCCGGGGTCGGCACAACGCTCACGTCCGAGAGCCACACGGCGTTGTCGCGCACAACGCTAACCGCCGTGTCGAGCACAAGCACCTGGGCGCTGTCTGGCa                   >  5:114117/1‑149 (MQ=255)
                                                                                                                                      gCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGt                >  8:119439/1‑149 (MQ=255)
                                                                                                                                        cACCGTAAGCCGGCAATTGCTTGGCGCTGAGTACGGCGCCCAGTGCCGGGGTCGGCAAAAAGGGCACGTCGGGGTGCCACACGGCGTTGTCGCGCACTTCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGGCTGGCACGTTg              <  6:114117/149‑1 (MQ=255)
                                                                                                                                           cGTAAGCAGTCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTg                      <  4:83229/138‑1 (MQ=255)
                                                                                                                                           cGTAAGCAGTCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTg                      <  4:78432/138‑1 (MQ=255)
                                                                                                                                           cGTAAGCAGTCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTg                      >  3:78432/1‑138 (MQ=255)
                                                                                                                                           cGTAAGCAGTCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTg                      >  3:83229/1‑138 (MQ=255)
                                                                                                                                                    gCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCggg  <  7:119439/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGGCGTGGTGCCAAAGCGCTCCAGCGGGAACGACTTGGTGAAGTCGTGGGTAGCGGTCAGCCCGTCCAGCATTTCACGCAGCGGCGCCGACAAGGCCTCGAAGGCCGCAATGCCGCTGGCCCAAAGGGTATCGCCACCGTAAGCAGGCAATTGCTTGGCGCTGAGTACCGCGCCCAGTGCCGGGGTCGGCAAAAAGGTCACGTCGGTGTGCCACACGGCGTTGTCGCGCACATCGGTAACCGCCGTGTCGAGCACAAGCACCTGTGGGGTGTCTGGCACGTTGGGGTAAATCGGG  >  NC_002947/284951‑285246

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: