BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	499,460	G→T	5.3%	L172I (CTC→ATC) ‡	aes ←	acetyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	499,460	0	G	T	5.3%	454.8 / 8.4	129	L172I (CTC→ATC) ‡	aes	acetyl esterase
Reads supporting (aligned to +/- strand): ref base G (58/64); new base T (3/4); total (61/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAA  >  NC_000913/499312‑499608

aaaTCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAg                                                                                                                                                      >  1:131001/1‑149 (MQ=255)

    ccTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGcc                                                                                                                                                  <  1:916602/149‑1 (MQ=255)

       gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAg                                                                                                                                                >  3:889625/1‑148 (MQ=255)

         tGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGCTCGCCAgtg                                                                                                                                              <  4:39526/148‑3 (MQ=255)

         tGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGCTCGCCAgtg                                                                                                                                              <  4:373863/148‑3 (MQ=255)

            gTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGc                                                                                                                                          <  1:340577/149‑1 (MQ=255)

              tAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCa                                                                                                                                         >  4:540691/1‑148 (MQ=255)

              tAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCa                                                                                                                                         >  3:81769/1‑148 (MQ=255)

              tAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCa                                                                                                                                         >  3:994456/1‑148 (MQ=255)

                aGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCt                                                                                                                                      >  2:604176/1‑149 (MQ=255)

                 gCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTg                                                                                                                                     >  1:35014/1‑149 (MQ=255)

                 gCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTg                                                                                                                                     >  2:449787/1‑149 (MQ=255)

                  ccATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTg                                                                                                                                     >  3:467233/1‑148 (MQ=255)

                  ccATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGc                                                                                                                                    >  1:989777/1‑149 (MQ=255)

                        cAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAAt                                                                                                                               <  4:977308/148‑1 (MQ=255)

                           acacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCAcc                                                                                                                           <  2:1025007/149‑1 (MQ=255)

                            cacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATg                                                                                                                                            >  4:202944/1‑131 (MQ=255)

                                gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCaa                                                                                                                      >  1:434748/1‑149 (MQ=255)

                                   cAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAACCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGc                                                                                                                   <  1:355497/149‑1 (MQ=255)

                                         acgacgAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGc                                                                                                             <  2:150498/149‑1 (MQ=255)

                                          cgacgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAAAGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGc                                                                                                             >  4:179473/1‑148 (MQ=255)

                                           gacgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCgg                                                                                                           <  1:762298/149‑1 (MQ=255)

                                            acgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCgg                                                                                                           >  4:267805/1‑148 (MQ=255)

                                             cgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCggg                                                                                                          <  4:223745/148‑1 (MQ=255)

                                             cgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCggg                                                                                                          <  3:267805/148‑1 (MQ=255)

                                                gTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACa                                                                                                       <  4:331570/148‑1 (MQ=255)

                                                gTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACa                                                                                                       <  4:842107/148‑1 (MQ=255)

                                                   aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGGACTGGCGGGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACatat                                                                                                    <  4:194246/148‑1 (MQ=255)

                                                   aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAtt                                                                                                   >  2:919356/1‑149 (MQ=255)

                                                    cGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAtt                                                                                                   >  4:707706/1‑148 (MQ=255)

                                                     ggAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttg                                                                                                  >  4:1002947/1‑148 (MQ=255)

                                                       aaTCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgat                                                                                                <  3:201009/148‑1 (MQ=255)

                                                          cccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatttg                                                                                             <  3:70377/148‑1 (MQ=255)

                                                          cccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgattag                                                                                             <  4:218140/148‑3 (MQ=255)

                                                          cccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgacttg                                                                                             <  3:70361/148‑1 (MQ=255)

                                                           ccGTAATCCGTAAAGCCCACACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatttgat                                                                                           >  2:691259/1‑149 (MQ=255)

                                                            cGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATa                                                                                          <  2:469558/149‑1 (MQ=255)

                                                             gTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATaa                                                                                         >  1:825245/1‑149 (MQ=255)

                                                              tAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAt                                                                                        >  1:903712/1‑149 (MQ=255)

                                                                  ccGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAtcctc                                                                                    <  1:663629/149‑1 (MQ=255)

                                                                  ccGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAtcct                                                                                     <  4:973237/148‑1 (MQ=255)

                                                                       aaGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCct                                                                               >  2:314251/1‑149 (MQ=255)

                                                                          cccATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGa                                                                                                         >  2:512464/1‑120 (MQ=255)

                                                                          cccATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGa                                                                                                         <  1:512464/120‑1 (MQ=255)

                                                                             aTACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCtggt                                                                         >  1:269686/1‑149 (MQ=255)

                                                                                ccACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGa                                                                      >  1:508296/1‑149 (MQ=255)

                                                                                 cacaGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGaa                                                                     >  1:231227/1‑149 (MQ=255)

                                                                                  acaGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGaa                                                                     <  4:1007067/148‑1 (MQ=255)

                                                                                   caGCAAAGCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGaaa                                                                    <  3:679921/148‑1 (MQ=255)

                                                                                        aaaCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATaacaa                                                              >  2:652761/1‑149 (MQ=255)

                                                                                         aaCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATaacaa                                                              >  3:897409/1‑148 (MQ=255)

                                                                                         aaCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAACCTCCGCCTGCTGGTGGAAATaacaa                                                              <  3:707706/148‑1 (MQ=255)

                                                                                           cGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAacaaca                                                            <  3:776198/148‑1 (MQ=255)

                                                                                             cccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagc                                                         <  1:610682/149‑1 (MQ=255)

                                                                                              ccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCAc                                                                                                                            >  2:542075/1‑81 (MQ=255)

                                                                                              ccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCAc                                                                                                                            <  1:542075/81‑1 (MQ=255)

                                                                                              ccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagca                                                        <  2:131001/149‑1 (MQ=255)

                                                                                                gCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcagc                                                      <  1:919356/149‑1 (MQ=255)

                                                                                                gCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcag                                                       >  4:629514/1‑148 (MQ=255)

                                                                                                     tttACCGCAATCGAGCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCcaca                                                  <  3:295420/148‑1 (MQ=255)

                                                                                                          cGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTcc                                            <  2:903712/149‑1 (MQ=255)

                                                                                                          cGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTcc                                            <  2:825245/149‑1 (MQ=255)

                                                                                                           gCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCct                                           <  1:160688/149‑1 (MQ=255)

                                                                                                              aTCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTa                                        >  1:23146/1‑149 (MQ=255)

                                                                                                               tCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTAt                                       >  2:152537/1‑149 (MQ=255)

                                                                                                                     tGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCtggt                                                                         >  1:430774/1‑109 (MQ=255)

                                                                                                                     tGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCtggt                                                                         <  2:430774/109‑1 (MQ=255)

                                                                                                                     cggcaaTCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCtt                                  >  4:407746/6‑148 (MQ=255)

                                                                                                                      gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTg                                 >  3:854699/1‑148 (MQ=255)

                                                                                                                      gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGc                                <  1:449787/149‑1 (MQ=255)

                                                                                                                           tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGc                                >  2:871392/1‑144 (MQ=255)

                                                                                                                           tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaa                            >  3:447356/1‑148 (MQ=255)

                                                                                                                           tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaa                            >  4:458126/1‑148 (MQ=255)

                                                                                                                            cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaa                           >  3:393097/1‑148 (MQ=255)

                                                                                                                             aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                         >  2:29030/1‑149 (MQ=255)

                                                                                                                              cGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                         >  3:650850/1‑148 (MQ=255)

                                                                                                                              cGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                         >  4:1058818/1‑148 (MQ=255)

                                                                                                                                 aaCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcag                                                       <  2:977525/115‑1 (MQ=255)

                                                                                                                                 aaCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcag                                                       >  1:977525/1‑115 (MQ=255)

                                                                                                                                 aaCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACgcgc                     <  2:396285/149‑1 (MQ=255)

                                                                                                                                  aCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACgcgc                     >  3:941001/1‑148 (MQ=255)

                                                                                                                                    cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCctgc                                                                             <  2:858436/90‑1 (MQ=255)

                                                                                                                                    cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCtt                   <  3:313986/148‑1 (MQ=255)

                                                                                                                                    cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTc                  >  2:773910/1‑149 (MQ=255)

                                                                                                                                    cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCAGCctgc                                                                             >  1:858436/1‑90 (MQ=255)

                                                                                                                                      caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCGCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAg                <  1:837594/149‑1 (MQ=255)

                                                                                                                                       aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaaa                          >  2:1057931/1‑138 (MQ=255)

                                                                                                                                       aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaaa                          <  1:1057931/138‑1 (MQ=255)

                                                                                                                                        aCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcgcg                      >  3:378351/1‑141 (MQ=255)

                                                                                                                                         cGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGt              >  3:1047935/1‑148 (MQ=255)

                                                                                                                                          atgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTACTCTATCGCTTGCGGAAAACGCGCTTCAGGTg             >  3:494839/4‑148 (MQ=255)

                                                                                                                                           tgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTc                                             <  3:276954/113‑1 (MQ=255)

                                                                                                                                           tgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcgcg                      >  4:31862/3‑138 (MQ=255)

                                                                                                                                           tgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcgcg                      <  3:31862/136‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAgg        <  2:434748/149‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg         <  3:849646/148‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg         <  3:622237/148‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg         <  3:179473/148‑1 (MQ=255)

                                                                                                                                                 cGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgtg     <  1:652761/149‑1 (MQ=255)

                                                                                                                                                    gcgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa  <  1:1011747/149‑1 (MQ=255)

AAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAA  >  NC_000913/499312‑499608

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