BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,516	G→T	55.4%	T236K (ACG→AAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,516	0	G	T	55.4%	28.4 / inf	202	T236K (ACG→AAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (49/41); new base T (50/62); total (99/103)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

Only 100 of 205 total aligned reads displayed.

TTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGC  >  NC_000913/4116369‑4116664

     tGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTgccgc                                                                                                                                                >  4:169094/1‑149 (MQ=255)

           cACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTgc                                                                                                                                          <  3:673865/149‑1 (MQ=255)

           cACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgc                                                                                                                                          >  3:1123133/1‑149 (MQ=255)

            aCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTgcc                                                                                                                                         >  1:1121860/1‑149 (MQ=255)

            aCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTgcc                                                                                                                                         >  2:1096206/1‑149 (MQ=255)

            aCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgcc                                                                                                                                         >  4:747489/1‑149 (MQ=255)

            aCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgcc                                                                                                                                         >  3:1236784/1‑149 (MQ=255)

             ccGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccg                                                                                                                                        >  2:174247/1‑149 (MQ=255)

             ccGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccg                                                                                                                                        <  1:719453/149‑1 (MQ=255)

                cTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCa                                                                                                                                     >  2:1339543/1‑149 (MQ=255)

                 tttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCaa                                                                                                                                    <  1:174247/149‑1 (MQ=255)

                      cGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCTGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTa                                                                                                                               <  4:1074393/149‑1 (MQ=255)

                      cGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTa                                                                                                                               <  4:761672/149‑1 (MQ=255)

                       gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAg                                                                                                                              >  4:1386092/1‑149 (MQ=255)

                       gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAg                                                                                                                              >  2:1248053/1‑149 (MQ=255)

                            tgtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTCCACTCACAACTGACCAAACAGCTCTGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGGCTTTGCCGCCAATTTTAGTCTGa                                                                                                                         >  1:1295069/1‑149 (MQ=255)

                              ttCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAcc                                                                                                                       >  2:747434/1‑149 (MQ=255)

                              gtcatcaGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGAcc                                                                                                                       >  4:1525680/2‑149 (MQ=255)

                                    aGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCaa                                                                                                                                    >  2:1409375/1‑130 (MQ=255)

                                       aTAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTc                                                                                                              >  3:57219/1‑149 (MQ=255)

                                        tAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCg                                                                                                             >  1:1482354/1‑149 (MQ=255)

                                        tAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCg                                                                                                             >  4:1494157/1‑149 (MQ=255)

                                            gcagcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaag                                                                                                         >  3:140560/1‑149 (MQ=255)

                                              agcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGgaagaa                                                                                                       <  3:194358/149‑1 (MQ=255)

                                               gcCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac                                                                                                      <  1:851543/149‑1 (MQ=255)

                                                  aGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                                                                                                   <  1:1145984/149‑1 (MQ=255)

                                                   gTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAacgacg                                                                                                  >  1:625347/1‑149 (MQ=255)

                                                   gTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAacgacg                                                                                                  >  3:101679/1‑149 (MQ=255)

                                                       cATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                                                                                              >  2:445824/1‑149 (MQ=255)

                                                         tAGATGTTCTTCGCCATCCCTTCTTTCTCGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc                                                                                            <  4:1013612/149‑1 (MQ=255)

                                                          aGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGGGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATTCTTCGGAAGAACGACTCCCTTCt                                                                                           >  3:404231/1‑149 (MQ=255)

                                                            gtgtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTgg                                                                                         <  1:1579952/149‑1 (MQ=255)

                                                            gtgtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGgaa                                                                                                          >  1:1076071/1‑132 (MQ=255)

                                                            gtgtTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGgaa                                                                                                          <  2:1076071/132‑1 (MQ=255)

                                                                  ttcGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAt                                                                                   <  2:1281335/149‑1 (MQ=255)

                                                                  ttcGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAt                                                                                   >  4:954791/1‑149 (MQ=255)

                                                                    cGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtct                                                                                 >  3:1530273/1‑149 (MQ=255)

                                                                     gCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtctc                                                                                <  3:587329/149‑1 (MQ=255)

                                                                         tCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtc                                                                                  <  1:394291/143‑1 (MQ=255)

                                                                         tCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtc                                                                                  >  2:394291/1‑143 (MQ=255)

                                                                         tCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgcgc                                                                            >  3:869641/1‑149 (MQ=255)

                                                                          cccTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgcgcg                                                                           >  2:1446685/1‑149 (MQ=255)

                                                                          cccTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgcgcg                                                                           >  4:400733/1‑149 (MQ=255)

                                                                           ccTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg                                                                          >  3:95351/1‑149 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCTCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTAGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                         >  1:305318/1‑149 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAGTTCTGGCAGCATCTCGCGCGGa                                                                         <  1:967664/149‑1 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                         >  2:1512300/1‑149 (MQ=255)

                                                                            cttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa                                                                         >  3:1158897/1‑149 (MQ=255)

                                                                                 ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATc                                                                    <  1:357111/149‑1 (MQ=255)

                                                                                  tCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATcc                                                                   <  3:1093340/149‑1 (MQ=255)

                                                                                  tCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATcc                                                                   <  1:903022/149‑1 (MQ=255)

                                                                                   cACGCACAACTGACCAAACAGCGCGGCCTGATGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCAGTTCTGGCAGCATCTCGCGCGGAATATCCa                                                                  <  3:1277452/149‑1 (MQ=255)

                                                                                       cacaACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCAGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGAACGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGGGCGGAATATcgaacag                                                              <  4:604859/149‑7 (MQ=255)

                                                                                        acaACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                                 <  1:974426/145‑1 (MQ=255)

                                                                                        acaACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACt                                                             >  3:264237/1‑149 (MQ=255)

                                                                                        acaACTGAACAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                                 >  2:974426/1‑145 (MQ=255)

                                                                                          aaCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTc                                                           <  2:614373/149‑1 (MQ=255)

                                                                                          aaCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTc                                                           <  1:565873/149‑1 (MQ=255)

                                                                                          aaCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTc                                                           >  1:754784/1‑149 (MQ=255)

                                                                                           aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATATCGCGCGGAATATCCAGCACTTcc                                                          <  1:1274390/149‑1 (MQ=255)

                                                                                            cggACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGAATCTCGCGCGGAATATCCAGAACTTCCa                                                         <  4:372965/147‑1 (MQ=255)

                                                                                            cTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCa                                                         <  4:264237/149‑1 (MQ=255)

                                                                                                ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCTCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCa                                                                  >  2:772567/1‑136 (MQ=255)

                                                                                                ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAt                                                     >  2:556071/1‑149 (MQ=255)

                                                                                                 cAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtt                                                    >  1:1410760/1‑149 (MQ=255)

                                                                                                  aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAttt                                                   >  3:648067/1‑149 (MQ=255)

                                                                                                  aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTCCCAGCAttt                                                   >  3:1192416/1‑149 (MQ=255)

                                                                                                   aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtttt                                                  <  2:248880/149‑1 (MQ=255)

                                                                                                   aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAtttt                                                  <  1:1180169/149‑1 (MQ=255)

                                                                                                    aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaa                                                                                                          >  2:693379/1‑92 (MQ=255)

                                                                                                    aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTg                                                 <  1:1339543/149‑1 (MQ=255)

                                                                                                      aGCGCGGCCTGCTGGTCACTGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtc                                               <  4:1134568/149‑1 (MQ=255)

                                                                                                      aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtc                                               <  2:669490/149‑1 (MQ=255)

                                                                                                           ggCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTccc                                          <  3:1494157/149‑1 (MQ=255)

                                                                                                            gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCa                                         <  1:1446685/149‑1 (MQ=255)

                                                                                                            gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCa                                         <  2:1430337/149‑1 (MQ=255)

                                                                                                               tgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGGATATCCAGCACTTCCAGCATTTTg                                                 <  1:1418210/138‑1 (MQ=255)

                                                                                                                  tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAg                                   <  3:784469/149‑1 (MQ=255)

                                                                                                                   ggTCACCGGCGATCCCGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAgg                                  <  2:467370/149‑1 (MQ=255)

                                                                                                                   ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAgg                                  <  2:625347/149‑1 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGCGATTGGAATACGCTTGACGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGCAAGTACGACGCACTTTTGGCACCATCTCGCGCGGCATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                              <  3:970754/149‑1 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                              >  3:1570224/1‑149 (MQ=255)

                                                                                                                       aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAt                              >  2:312947/1‑149 (MQ=255)

                                                                                                                                 ccGGAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaa                    <  1:1232068/149‑1 (MQ=255)

                                                                                                                                 ccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaa                    >  4:1202329/1‑149 (MQ=255)

                                                                                                                                 ccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaa                    >  3:1322603/1‑149 (MQ=255)

                                                                                                                                 ccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGACGTTGaa                    <  2:1129093/149‑1 (MQ=255)

                                                                                                                                 ccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaa                    >  2:927196/1‑149 (MQ=255)

                                                                                                                                  cGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTgtcgtc                                            >  1:208879/1‑124 (MQ=255)

                                                                                                                                   ggAGATTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGaaca                  >  1:783533/1‑149 (MQ=255)

                                                                                                                                   ggAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAtt                                                                                                                                  >  4:128028/1‑36 (MQ=255)

                                                                                                                                    gagaTTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                                 <  3:476415/101‑1 (MQ=255)

                                                                                                                                    gagaTTGGAATACGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAg                                                                 >  4:476415/1‑101 (MQ=255)

                                                                                                                                      gaTTGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAacaaca               <  2:1103465/149‑1 (MQ=255)

                                                                                                                                        ttGGAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATg             <  1:1096206/149‑1 (MQ=255)

                                                                                                                                           gAATACGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTa          <  1:1587932/149‑1 (MQ=255)

                                                                                                                                               aCGCTTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACgag      <  4:1260779/149‑1 (MQ=255)

                                                                                                                                                   gTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGc  >  1:654065/1‑149 (MQ=255)

                                                                                                                                                   gTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGc  <  3:536244/149‑1 (MQ=255)

                                                                                                                                                   gTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGc  <  4:869641/149‑1 (MQ=255)

TTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGC  >  NC_000913/4116369‑4116664

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