Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,116,536 G→T 12.6% N229K (AAC→AAA glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,116,5360GT12.6% 451.4 / 44.1 135N229K (AAC→AAAglpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (61/57);  new base T (9/8);  total (70/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 138 total aligned reads displayed.

TCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACAT  >  NC_000913/4116390‑4116684
                                                                                                                                                  |                                                                                                                                                    
 cGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTa                                                                                                                                                   <  3:330913/149‑1 (MQ=255)
 cGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTa                                                                                                                                                   <  4:591046/149‑1 (MQ=255)
   ccAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAGCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGt                                                                                                                                                 <  4:655958/149‑1 (MQ=255)
     aGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCt                                                                                                                                               <  2:400985/149‑1 (MQ=255)
        gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGt                                                                                                                                                 >  2:128452/1‑144 (MQ=255)
        gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGt                                                                                                                                                 <  1:128452/144‑1 (MQ=255)
        gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                                                                                                            >  3:781418/1‑149 (MQ=255)
             tcaGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                                                                                                                       >  3:68263/1‑149 (MQ=255)
               aGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATAc                                                                                                                                     >  2:63940/1‑149 (MQ=255)
                     aaGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCt                                                                                                                                               >  3:600846/1‑133 (MQ=255)
                     aaGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGga                                                                                                                               <  1:539176/149‑1 (MQ=255)
                      agcagcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGgaa                                                                                                                              <  3:396454/149‑1 (MQ=255)
                      agcagcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGgaa                                                                                                                              <  1:554811/149‑1 (MQ=255)
                          gcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac                                                                                                                          >  3:728498/1‑149 (MQ=255)
                            cAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacga                                                                                                                        <  2:499367/149‑1 (MQ=255)
                             aGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                                                                                                                       <  4:68263/149‑1 (MQ=255)
                             aGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                                                                                                                       <  2:39291/149‑1 (MQ=255)
                                  cATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                                                                                                                  >  1:516497/1‑149 (MQ=255)
                                    tAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc                                                                                                                >  1:41865/1‑149 (MQ=255)
                                       gtgtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTgg                                                                                                             >  4:425787/1‑149 (MQ=255)
                                        tgtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                                                                                                            <  2:218764/117‑1 (MQ=255)
                                        tgtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                                                                                                            >  1:218764/1‑117 (MQ=255)
                                         gtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                                                                                                           <  4:750309/149‑1 (MQ=255)
                                             ttcGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAt                                                                                                       >  3:706931/1‑149 (MQ=255)
                                                gCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGtt                                                                                                                                                    >  1:783787/1‑101 (MQ=255)
                                                gCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGtt                                                                                                                                                    <  2:783787/101‑1 (MQ=255)
                                                 ccATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcg                                                                                                   >  1:322259/1‑149 (MQ=255)
                                                 ccATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcg                                                                                                   >  1:511225/1‑149 (MQ=255)
                                                      ccTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtctc                                                                                                    >  4:342989/1‑143 (MQ=255)
                                                       cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                                             >  1:112604/1‑149 (MQ=255)
                                                       cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                                             >  3:454628/1‑149 (MQ=255)
                                                           tttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat                                                                                         >  3:37606/1‑149 (MQ=255)
                                                            ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATc                                                                                        >  1:422823/1‑149 (MQ=255)
                                                               aCGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                                                     <  3:489354/149‑1 (MQ=255)
                                                               aCGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                                                     >  1:258904/1‑149 (MQ=255)
                                                               aCGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                                                     <  2:511225/149‑1 (MQ=255)
                                                                 gCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGAAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCa                                                                                   >  4:146719/1‑149 (MQ=255)
                                                                   acaACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACt                                                                                 >  1:219250/1‑149 (MQ=255)
                                                                     aaCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTc                                                                               >  3:325438/1‑149 (MQ=255)
                                                                          aCCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCa                                                                          >  1:347280/1‑149 (MQ=255)
                                                                          aCCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCa                                                                          >  1:408986/1‑149 (MQ=255)
                                                                            cAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtt                                                                        >  1:31051/1‑149 (MQ=255)
                                                                              aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtttt                                                                      >  1:271941/1‑149 (MQ=255)
                                                                              aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtttt                                                                      >  3:575085/1‑149 (MQ=255)
                                                                                 aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtc                                                                   <  3:339691/149‑1 (MQ=255)
                                                                                 aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtc                                                                   <  2:148586/149‑1 (MQ=255)
                                                                                  gcgcGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtcg                                                                  >  4:717539/1‑149 (MQ=255)
                                                                                         ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGa                                                                                                                                             <  4:28523/67‑1 (MQ=255)
                                                                                         ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGa                                                                                                                                             >  3:28523/1‑67 (MQ=255)
                                                                                              ggTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAgg                                                      >  1:122983/1‑149 (MQ=255)
                                                                                              ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacga                                                                                                                        >  2:643462/1‑83 (MQ=255)
                                                                                              ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacga                                                                                                                        <  1:643462/83‑1 (MQ=255)
                                                                                              ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAgg                                                      <  4:427772/149‑1 (MQ=255)
                                                                                               gTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAggg                                                     >  3:449936/1‑149 (MQ=255)
                                                                                                 cACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTa                                                   <  2:422823/149‑1 (MQ=255)
                                                                                                  aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                                                  >  3:686099/1‑149 (MQ=255)
                                                                                                  aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                                                  >  4:319958/1‑149 (MQ=255)
                                                                                                  aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                                                  <  2:516497/149‑1 (MQ=255)
                                                                                                  aCCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                                                  <  2:288549/149‑1 (MQ=255)
                                                                                                    cGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                                                                                                                                <  2:60163/69‑1 (MQ=255)
                                                                                                    cGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                                                                                                                                >  1:60163/1‑69 (MQ=255)
                                                                                                    cGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTcc                                                               >  1:528953/1‑134 (MQ=255)
                                                                                                      gCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                                                                                                           >  4:194025/1‑88 (MQ=255)
                                                                                                      gCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                                                                                                           <  3:194025/88‑1 (MQ=255)
                                                                                                      gCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGAt                                              <  2:806325/149‑1 (MQ=255)
                                                                                                      gCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGAt                                              <  1:244069/149‑1 (MQ=255)
                                                                                                        gATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGt                                            >  3:742140/1‑149 (MQ=255)
                                                                                                        gATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGt                                            >  2:149610/1‑149 (MQ=255)
                                                                                                         aTCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGtt                                           >  3:257226/1‑149 (MQ=255)
                                                                                                         aTCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGtt                                           >  3:54963/1‑149 (MQ=255)
                                                                                                         aTCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGtt                                           <  1:715183/149‑1 (MQ=255)
                                                                                                            ccGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaa                                        >  4:552879/1‑149 (MQ=255)
                                                                                                                 gaTTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAacaaca                                   <  4:325438/149‑1 (MQ=255)
                                                                                                                  aTTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACAt                                  >  4:214796/1‑149 (MQ=255)
                                                                                                                  aTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACAt                                  >  4:376699/1‑149 (MQ=255)
                                                                                                                   ttGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATg                                 <  2:347280/149‑1 (MQ=255)
                                                                                                                   ttGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATg                                 <  4:70907/149‑1 (MQ=255)
                                                                                                                   ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATg                                 <  4:54963/149‑1 (MQ=255)
                                                                                                                     tgAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGt                               <  2:271941/148‑1 (MQ=255)
                                                                                                                     ggAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGt                               <  4:575085/149‑1 (MQ=255)
                                                                                                                     ggAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGt                               >  3:301623/1‑149 (MQ=255)
                                                                                                                      gAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTa                              >  2:431351/1‑149 (MQ=255)
                                                                                                                      gAATACGCTGGCCGCATTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGCCGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTAAAGCATTTTGACGTCCCAGTCCAGGGTATGGATGTTCAACAACATGGtt                              >  2:231010/1‑148 (MQ=255)
                                                                                                                      gAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTa                              <  4:707347/149‑1 (MQ=255)
                                                                                                                      gAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTACTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTa                              <  4:454628/149‑1 (MQ=255)
                                                                                                                           cgctTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACgaga                         <  3:425787/149‑1 (MQ=255)
                                                                                                                            gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGagag                        <  2:408986/149‑1 (MQ=255)
                                                                                                                              ttGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGc                      <  4:359508/148‑1 (MQ=255)
                                                                                                                              ttGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGc                      >  3:767036/2‑149 (MQ=255)
                                                                                                                                   gccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTgg                 >  2:546504/1‑149 (MQ=255)
                                                                                                                                    ccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                <  3:146719/149‑1 (MQ=255)
                                                                                                                                           cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATc                                                                                        <  3:511637/70‑1 (MQ=255)
                                                                                                                                           cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATc                                                                                        >  4:511637/1‑70 (MQ=255)
                                                                                                                                            gccAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCgg        >  2:531503/1‑149 (MQ=255)
                                                                                                                                            gccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAgg                                                      <  1:610769/103‑1 (MQ=255)
                                                                                                                                             ccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGt       <  4:231857/149‑1 (MQ=255)
                                                                                                                                               aaTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTca     >  2:784864/1‑149 (MQ=255)
                                                                                                                                                  gTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTc      >  3:797903/1‑145 (MQ=255)
                                                                                                                                                  gTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTc      <  4:797903/145‑1 (MQ=255)
                                                                                                                                                  gTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACAt  <  4:33003/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACAT  >  NC_000913/4116390‑4116684

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: