Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,312,192	T→C	100%	M19V (ATG→GTG)	yciA ←	acyl‑CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,192	0	T	C	98.1%	333.1 / ‑3.3	100	M19V (ATG→GTG)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (46/52);  total (47/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC  >  NC_000913/1312096‑1312288                                                                                                 |                                                                                                 gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGc                                                                                               <  2:3637465/100‑1 (MQ=255) gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGc                                                                                               <  1:2825700/100‑1 (MQ=255)        ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGATATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTa                                                                                        <  1:2483381/100‑1 (MQ=255)        ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTa                                                                                        <  2:880621/100‑1 (MQ=255)            tttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTa                                                                                    >  1:3746630/1‑100 (MQ=255)              tGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaa                                                                                  >  2:1024484/1‑100 (MQ=255)              tGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaa                                                                                  >  1:817760/1‑100 (MQ=255)              tGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaa                                                                                  >  1:3464652/1‑100 (MQ=255)                ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAAc                                                                                <  1:1056427/100‑1 (MQ=255)                ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCAATGGTATCGGCGGGCACGGCTAAAGTACGTAAAAc                                                                                >  1:712424/1‑100 (MQ=255)                 cAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAAc                                                                                <  2:2174987/99‑1 (MQ=255)                 cAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACa                                                                               >  2:993248/1‑100 (MQ=255)                 cAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACa                                                                               >  2:1143736/1‑100 (MQ=255)                  aGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACa                                                                               <  2:1796417/99‑1 (MQ=255)                   gAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAg                                                                             >  1:1796881/1‑100 (MQ=255)                     aTAGCGCCGCCAATATCCATTTGTGACATTAACCAGCCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGAt                                                                           <  1:1971633/100‑1 (MQ=255)                       aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCg                                                                         >  2:2678529/1‑100 (MQ=255)                       aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCg                                                                         >  2:416799/1‑100 (MQ=255)                       aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCg                                                                         <  1:3020964/100‑1 (MQ=255)                       aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAACTACGGAACACaa                                                                              >  2:2403711/1‑95 (MQ=255)                        gcgcCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGc                                                                        <  2:1994754/100‑1 (MQ=255)                         cgccgcCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGcc                                                                       <  1:320686/100‑1 (MQ=255)                          gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGGCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGccc                                                                      <  1:3750343/100‑1 (MQ=255)                           ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGccc                                                                      >  1:775957/1‑99 (MQ=255)                                   atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                             <  2:1796871/100‑1 (MQ=255)                                   atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                             >  1:1850759/1‑100 (MQ=255)                                   atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                             >  2:1346740/1‑100 (MQ=255)                                   atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                             <  1:969132/100‑1 (MQ=255)                                      cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTa                                                          >  1:799282/1‑100 (MQ=255)                                         ttGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAtgt                                                       >  2:1081162/1‑100 (MQ=255)                                          tgtgACGTTAACCAGCCACCAGAGAGGTAAATATTCGAGTTGGTATCGGCGGGCACTGCTAAAGTACGTAAGACAAGAGGGCCCTGGGGGACGTTAttt                                                       <  2:381697/99‑3 (MQ=255)                                          tgtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAtgtg                                                      >  2:1695660/1‑100 (MQ=255)                                           gtgACATTAACCAACCGCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                     <  2:1346790/100‑1 (MQ=255)                                           gtgACATTAACCAACCACCAGAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                     >  2:2332699/1‑100 (MQ=255)                                           gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                     >  1:1445715/1‑100 (MQ=255)                                           gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                     <  2:2772903/100‑1 (MQ=255)                                           gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                     >  1:1780795/1‑100 (MQ=255)                                           gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                     <  1:2829599/100‑1 (MQ=255)                                             gACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttg                                                   <  1:3629405/100‑1 (MQ=255)                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                                  <  1:2800443/100‑1 (MQ=255)                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                                  <  2:3191407/100‑1 (MQ=255)                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                                  <  1:1388324/100‑1 (MQ=255)                                               cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                                 >  1:3698010/1‑100 (MQ=255)                                               cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                                 >  2:2734218/1‑100 (MQ=255)                                               cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                                 >  2:643037/1‑100 (MQ=255)                                                aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAg                                                <  2:2027909/100‑1 (MQ=255)                                                 ttAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGAACGGCTAAAGTACGTAAAACAAGATAGCCCCGAGGGACGTTATGTGTTGTGGa                                               >  2:1406055/1‑100 (MQ=255)                                                  tAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCCGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGGCGTTATGTGTTGTAGa                                               >  2:3691223/1‑99 (MQ=255)                                                  tAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGAc                                              <  1:1213805/100‑1 (MQ=255)                                                   aaccaaccACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACa                                             >  1:2977945/1‑100 (MQ=255)                                                     ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                           <  2:1456864/100‑1 (MQ=255)                                                      caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATgg                                          <  1:3166135/100‑1 (MQ=255)                                                      caaacaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGCCACGGCTAAACTAAGTAAAACAAAATCGCCATGAAGGACGTTATGTGTTGTAAACa                                             >  2:3025133/1‑97 (MQ=255)                                                       aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGt                                         <  1:2046209/100‑1 (MQ=255)                                                       aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGt                                         <  2:2722811/100‑1 (MQ=255)                                                       aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACCGCTAAAGTACGTAAAACAAGATCGCCCTGCGGGACGTTATGTGTTGTAAACAggg                                          >  2:3141203/1‑99 (MQ=255)                                                         ccacaaATGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAGGCTCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTa                                        <  1:1279003/99‑1 (MQ=255)                                                           accaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaaa                                     >  2:3124984/1‑100 (MQ=255)                                                            ccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAc                                    >  1:1257354/1‑100 (MQ=255)                                                            ccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAc                                    >  1:298437/1‑100 (MQ=255)                                                            ccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAc                                    >  2:1729613/1‑100 (MQ=255)                                                             caAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCGCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAcc                                   >  2:33645/1‑100 (MQ=255)                                                             caAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAcc                                   <  1:2624779/100‑1 (MQ=255)                                                              aaaGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCg                                  <  1:2829571/100‑1 (MQ=255)                                                              aaaGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCg                                  >  2:125267/1‑100 (MQ=255)                                                              aaaGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCg                                  <  1:224566/100‑1 (MQ=255)                                                                aGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGAc                                <  2:317473/100‑1 (MQ=255)                                                                 gATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACt                               <  1:3099230/100‑1 (MQ=255)                                                                  aTGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACtt                              <  2:202736/100‑1 (MQ=255)                                                                   tGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTa                             <  2:1806596/100‑1 (MQ=255)                                                                    gTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAACACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTaa                            >  2:78760/1‑100 (MQ=255)                                                                     tCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTaaa                           <  1:993257/100‑1 (MQ=255)                                                                      cACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGTCCTGAGGGACGTTATGTGTTGTAGAAATGGTAAAACCGGATTaaaa                          >  2:725471/1‑100 (MQ=255)                                                                      cACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCCGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTaaaa                          <  1:288127/100‑1 (MQ=255)                                                                      cACCATTGGCATTGGCATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTaaaa                          <  1:3760537/100‑1 (MQ=255)                                                                       aCCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAg                         <  2:1376501/100‑1 (MQ=255)                                                                       aCCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAg                         <  1:2990491/100‑1 (MQ=255)                                                                       aCCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGAAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAg                         <  1:2135777/100‑1 (MQ=255)                                                                            tggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGaaaa                    <  1:1037926/100‑1 (MQ=255)                                                                             ggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAAt                   >  1:3165570/1‑100 (MQ=255)                                                                              gcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATc                  >  1:3522488/1‑100 (MQ=255)                                                                               nattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATca                 >  2:349086/2‑100 (MQ=255)                                                                                attggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATcac                <  2:1776276/100‑1 (MQ=255)                                                                                 ttggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATcaca               >  2:764128/1‑100 (MQ=255)                                                                                 ttggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCAGGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATcaca               <  2:1446817/100‑1 (MQ=255)                                                                                  tggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAg              <  2:2977923/100‑1 (MQ=255)                                                                                  tggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAg              <  1:3425046/100‑1 (MQ=255)                                                                                  ccgatctggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAg              <  2:152343/93‑1 (MQ=255)                                                                                    gTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgc            <  2:648065/100‑1 (MQ=255)                                                                                      aTCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGGGGGGCGGTATGTGGTGTAGACAAGGTAAAAACGACTTAAAAGGAAAATTAAAGgcgc          >  2:2774474/1‑100 (MQ=255)                                                                                      aTCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGGCGTTATGTGTTGTAGACATGGTAAAACCGACTGAAAAGGAAAATCACAGgc            >  2:2257334/1‑98 (MQ=255)                                                                                      aTCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgcgc          >  1:3302273/1‑100 (MQ=255)                                                                                      aTCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgcgc          >  1:2611910/1‑100 (MQ=255)                                                                                         ggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAgagatcgg       >  1:152343/1‑93 (MQ=255)                                                                                         ggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGAt       <  1:379638/100‑1 (MQ=255)                                                                                           cggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATg      <  2:551019/99‑1 (MQ=255)                                                                                             ggCACGGCTAAAGTACGTAACACAAGATCGCCCGGACGGCCGTTATGTATGGTACACATGCTCAAACCGACTTAAAAGGAAAATCACAGGCgaga        >  2:3659669/1‑95 (MQ=255)                                                                                             ggCACGGCTAAAGTACGTAAAACAAGATCGCCCTTAGGTACGTTATGGGTTGTAGATATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGc     <  2:1870782/98‑1 (MQ=255)                                                                                             ggCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGct    <  1:1386662/99‑1 (MQ=255)                                                                                              gCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGctac  >  1:3083424/1‑100 (MQ=255)                                                                                                 |                                                                                                 GGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC  >  NC_000913/1312096‑1312288

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑