Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 961,710 G→C 7.0% H235Q (CAC→CAG)  selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947961,7100GC 7.0% 221.8 / 6.9 72H235Q (CAC→CAG) selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base C (3/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

TGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGG  >  NC_002947/961562‑961855
                                                                                                                                                    |                                                                                                                                                 
tGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAgg                                                                                                                                                   >  1:125343/1‑149 (MQ=255)
tGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAgg                                                                                                                                                   >  7:162274/1‑149 (MQ=255)
   tCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTgg                                                                                                                                                >  8:213158/1‑149 (MQ=255)
       gtgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTg                                                                                                                                                 <  1:238688/144‑1 (MQ=255)
        tgCCCGTAGCTGTCGTAGTTGCGGAGGGTTCCTCCGGGGATGCAGCCTTCGGCGAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGa                                                                                                                                            <  4:72263/148‑1 (MQ=255)
        tgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCGGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAg                                                                                                                                           >  5:184936/1‑149 (MQ=255)
        tgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCCTCGGCCACCTGGTCGAGCTGGCCGaa                                                                                                                                           <  3:353032/149‑2 (MQ=255)
                      ttAGTTGAGCAGGGTTACACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGACAGCTCGAACAGGTGGCCGAGCAGGACGAAGCAgg                                                                                                                             >  4:50470/2‑149 (MQ=255)
                         gTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTg                                                                                                                           <  7:213158/148‑1 (MQ=255)
                             cgcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                      >  7:204969/1‑149 (MQ=255)
                                aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGc                                                                                                                                >  5:175020/1‑136 (MQ=255)
                                aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGc                                                                                                                                <  6:175020/136‑1 (MQ=255)
                                aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCg                                                                                                                              <  7:57953/138‑1 (MQ=255)
                                aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCg                                                                                                                              >  8:57953/1‑138 (MQ=255)
                                    ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                               <  2:125343/149‑1 (MQ=255)
                                    ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                               <  4:164830/149‑1 (MQ=255)
                                      ccACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATgg                                                                                                              <  6:142461/148‑1 (MQ=255)
                                        aCCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCt                                                                                                           >  6:83141/1‑149 (MQ=255)
                                         ccGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCtt                                                                                                          <  5:83141/149‑1 (MQ=255)
                                                  cAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTc                                                                                                 <  8:204969/149‑1 (MQ=255)
                                                       ttCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGtt                                                                                            >  1:257670/1‑149 (MQ=255)
                                                         cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGttt                                                                                           <  6:67203/148‑1 (MQ=255)
                                                            ccAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCg                                                                                       >  2:270637/1‑149 (MQ=255)
                                                             cAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCg                                                                                       >  6:324289/1‑148 (MQ=255)
                                                               ggTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAg                                                                                    >  8:46993/1‑149 (MQ=255)
                                                                  aaTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                 >  3:328678/1‑149 (MQ=255)
                                                                   aTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                 <  5:324289/148‑1 (MQ=255)
                                                                     ggTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTg                                                                              >  2:252564/1‑149 (MQ=255)
                                                                     ggTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCCACCAGGTTGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCAAGGCCTTGACCCCGTCCAGTTTGCAGAAGCGGCTg                                                                              >  7:305081/1‑149 (MQ=255)
                                                                       tCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                            >  2:231197/1‑149 (MQ=255)
                                                                         gACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCgg                                                                          <  8:162274/149‑1 (MQ=255)
                                                                              tAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgtg                                                                     >  1:153814/1‑149 (MQ=255)
                                                                                ggcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCgg                                                                                                                    <  4:268508/100‑1 (MQ=255)
                                                                                ggcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCgg                                                                                                                    >  3:268508/1‑100 (MQ=255)
                                                                                ggcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGtt                                                                   >  3:248429/1‑149 (MQ=255)
                                                                                  caggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGa                                                                 <  3:50470/149‑1 (MQ=255)
                                                                                   aggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAg                                                                >  5:84330/1‑149 (MQ=255)
                                                                                     gcGTGGCACGGCGGCGTAGTTCGGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAggg                                                              <  6:184936/149‑1 (MQ=255)
                                                                                          gCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTc                                                                                                 <  8:253262/109‑1 (MQ=255)
                                                                                          gCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTc                                                                                                 >  7:253262/1‑109 (MQ=255)
                                                                                          gCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGcac                                                         <  2:257670/149‑1 (MQ=255)
                                                                                            aCGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACAt                                                       <  2:153814/149‑1 (MQ=255)
                                                                                                  gcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTc                                                 >  2:297409/1‑149 (MQ=255)
                                                                                                        ttCAGGTGCGCGGTGAGGCCGCTGCCTTCCGCCAGCACGACCAGGTTGCCGAGCAAGCCCAAACCGGTGACGTCCGTCATGGCCTTTACCCCCGCCACTTTTCCCAAACCGCTGCCGGGCGTGTTTAGGGTGCACAAACAGTACCGAGc                                           >  2:182289/1‑149 (MQ=255)
                                                                                                           aGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCCTGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                        >  8:105095/1‑149 (MQ=255)
                                                                                                           aGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                        <  1:182289/149‑1 (MQ=255)
                                                                                                           aGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                        >  8:164540/1‑149 (MQ=255)
                                                                                                             gTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGCGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGGGTTGAGGGTGCACATCCAGTCGCGAGCCAGGc                                      >  5:364164/1‑149 (MQ=255)
                                                                                                             gTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGc                                      >  3:158314/1‑149 (MQ=255)
                                                                                                               gcgcGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCt                                    >  6:220740/1‑149 (MQ=255)
                                                                                                               gcgcGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgc                                               >  3:161828/1‑137 (MQ=255)
                                                                                                               gcgcGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgc                                               <  4:161828/137‑1 (MQ=255)
                                                                                                                 gcgGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGCAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                  <  1:270637/149‑1 (MQ=255)
                                                                                                                 gcgGGGAGGCCGGTGCCTTCGGCCAGCTCGACCTGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGGCATGGCCTTGAGCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                  <  5:29950/149‑1 (MQ=255)
                                                                                                                   ggTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGGGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAgc                                                                                             >  2:307257/1‑87 (MQ=255)
                                                                                                                         gCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGc                           <  1:231197/148‑1 (MQ=255)
                                                                                                                            gCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCacac                       >  3:188566/1‑149 (MQ=255)
                                                                                                                              tGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGc                     <  5:220740/149‑1 (MQ=255)
                                                                                                                                ccTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGc                     <  6:84330/147‑1 (MQ=255)
                                                                                                                                  ttCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCt                 >  1:6573/1‑149 (MQ=255)
                                                                                                                                       cctgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGcc                                     <  7:143312/121‑1 (MQ=38)
                                                                                                                                       cctgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGcc                                     >  8:143312/4‑124 (MQ=38)
                                                                                                                                       ccAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc            >  7:105091/1‑149 (MQ=255)
                                                                                                                                       ccAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGACTTTACACCGTCCACTTAGCCGAAGACGCTGCCGGGCGTGTTGAGGGTACACATCCAGTCCCGAGCCAGGCCTAGGTCCTGCAAACGCAGcgag               >  6:283727/1‑143 (MQ=255)
                                                                                                                                       caaGCTCGACCAGGTGTCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGGTGCCGTGCGTGTTGAGGGTTCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTg               <  5:283727/144‑1 (MQ=255)
                                                                                                                                        ctgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCct           >  6:38784/3‑149 (MQ=255)
                                                                                                                                        ctgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCct           >  6:269676/3‑149 (MQ=255)
                                                                                                                                         aGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCctt          >  6:152657/1‑149 (MQ=255)
                                                                                                                                         aGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCctt          <  7:164540/149‑1 (MQ=255)
                                                                                                                                           cTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttct        <  4:188566/149‑1 (MQ=255)
                                                                                                                                              gACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCtttt     <  4:328678/149‑1 (MQ=255)
                                                                                                                                              gACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCtttt     >  3:29444/1‑149 (MQ=255)
                                                                                                                                                ccAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCg   <  8:105091/149‑1 (MQ=255)
                                                                                                                                                 cAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGTTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCgg  <  1:252564/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGG  >  NC_002947/961562‑961855

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