Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,213	0	G	A	44.1%	72.5 / 42.1	69	L12F (CTT→TTT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (7/31);  new base A (30/0);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 7.44e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCT  >  NC_000913/1312069‑1312363                                                                                                                                                 |                                                                                                                                                          cAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGc                                                                                                                                                      <  2:917475/150‑1 (MQ=255)  aaCCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATACGTAAAAAAAAATCGccc                                                                                                                                                    >  2:405225/1‑151 (MQ=255)  aaCCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAAAACGGAAAAAAAAATCGccc                                                                                                                                                    >  2:174034/1‑151 (MQ=255)     ccGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGa                                                                                                                                                 <  2:489781/151‑1 (MQ=255)      cGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGCTAAAATAAGGAAAAAAAAAACGccc                                                                                                                                                    >  1:591925/1‑147 (MQ=255)       gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAgg                                                                                                                                               <  2:771151/151‑1 (MQ=255)         acaGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGCTAAAAAAAGTAAAAAAAAAACGCCCCGAggg                                                                                                                                              >  2:1087842/1‑150 (MQ=255)         acaGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCCAAAGAAGTAAAAAAAAAACACCCCGAggg                                                                                                                                              >  2:833229/1‑150 (MQ=255)           aGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                                                                                                           <  2:1260522/151‑1 (MQ=255)           aGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGTAAAACAAAATCGCCCTGAggg                                                                                                                                              >  1:492150/1‑148 (MQ=255)               aCTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCTGAGGGagggtgt                                                                                                                                       >  1:407389/1‑145 (MQ=255)                   cgcgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCAGGCTAAAAAAAGTAAAAAAAAAAAGCCCTCAggg                                                                                                                                              >  1:60315/1‑140 (MQ=255)                   cgcgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                                                                                                   <  1:639227/151‑1 (MQ=255)                        ccGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCCAAAATAAGTAAAAAAAAAACGCCCTGCGGGGGGGTATGtgttgt                                                                                                                                >  2:389270/1‑149 (MQ=255)                         cGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGACTGGCTAAAGGAAGGAAAAAAAAATAGCCCCGAGGGGCGGTATTtgttgt                                                                                                                                >  2:789315/1‑148 (MQ=255)                          gTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGCTTAAAAAAGTTAAACAAAAACGCCCCGAGGGGCGtttttt                                                                                                                                     >  1:455223/1‑138 (MQ=255)                            gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAt                                                                                                                         <  2:1132570/151‑1 (MQ=255)                            gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATCGCCGAAGTACGGAAAACAAAATTCCCCCGTGGGGCGGTAtgt                                                                                                                                     >  2:780956/1‑140 (MQ=255)                              gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGGTAAAAAAAGGAAAAAAAAAACGCCCCCAGGGGCGGTTTGTGTTGTAG‑GAAT‑gg                                                                                                                      >  2:716596/1‑151 (MQ=255)                               cAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGt                                                                                                                     <  1:869750/151‑1 (MQ=255)                                    tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGGTAAAGAACGTAAAAAAAAAACGCCCCCAGGGGGGGTATGTGTTGTTG‑AAAA‑GGGaaaa                                                                                                                 >  2:1046333/1‑150 (MQ=255)                                       tttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCg                                                                                                              <  2:1257707/150‑1 (MQ=255)                                        tttctagatgaatgccgccAATATCCATTTGTGACGTTAACCAACCACCAAAGAAGTCAACATTGGCATTGGTAGCGGCGGGCATGGCTAAAGTACGTAAAACAAGAGCGTCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAAcc                                                                                                               <  1:233964/135‑1 (MQ=255)                                          gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACtt                                                                                                          <  2:685687/151‑1 (MQ=255)                                              gAATAGCGCCGCTAATATCCATTTGTGACATTAACCAACCACAAAAGCTGTCATCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTaaa                                                                                                       <  1:660590/150‑1 (MQ=255)                                              gAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTaaaa                                                                                                      <  1:1000673/151‑1 (MQ=255)                                              gAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTaaaa                                                                                                      <  1:202060/151‑1 (MQ=255)                                               aaTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGTAAAACAAAATCGCCCCGAGGGGCGTTATGTGTTGTAG‑ACAA‑GGTAAAAAACACTTaaaa                                                                                                      >  2:244506/1‑150 (MQ=255)                                                  aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCGAGGAGAAAGTACGTAAAAAAAAAACTCACGCGGGGAAGTTTTGTTTTGTAG‑AGAA‑GGGAAAACCAAATTtaaa                                                                                                      >  2:1159205/1‑147 (MQ=255)                                                  aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGCAAAAAAAGTAAAAAAAAAACGCCCTGAGGGGCGGTATGTGGTGTTG‑AAAT‑GGGAAAAACGACTTAAAAAGa                                                                                                   >  1:568045/1‑150 (MQ=255)                                                    cgccgcCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGaaaa                                                                                                <  2:673046/151‑1 (MQ=255)                                                      ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCAAAATAAGTAAAAAAAAAACGCCCTGAGGGGCGGTTTGTGTTTTAG‑AAAA‑GGGAAAAACCACCTAAAAAGaaaa                                                                                                >  2:144067/1‑149 (MQ=255)                                                         ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGGTAAAAAAAATAAAAAAAAAAAGCCCTGAGGGGGGGGATGTGTTGTAG‑GAAA‑GGGAAAAAAAACTTAAAAAGaaaaaaa                                                                                             >  1:38686/1‑146 (MQ=255)                                                          cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAAGAAAAAAAAAAAGCCCCTAGGGGGGGTTTGTGTTGTTG‑ACATGGGAAAAAACCCCTTAAAAGAaaaaaaaa                                                                                            >  1:240357/1‑146 (MQ=255)                                                            atatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACAGGTAAAGAAAGGAAAAAAAAAAAGCCCTGAGGGGGGGTATGTGGTTTGG‑AAAT‑GGTAAAAACGACTTAAAAGGGAAAAca                                                                                             >  2:657077/1‑146 (MQ=255)                                                              atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGgcg                                                                                       <  1:863624/150‑1 (MQ=255)                                                                  aTTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAGACGTAAAAAAAAAACGCCCCGTGGGGCGGTATGTGTTGTAG‑ACAA‑GGGAAAAACGACTTAAAAAGAAAATAACAGGCGCGGTg                                                                                  >  2:648030/1‑151 (MQ=255)                                                                   tttGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCAAAAAACAGAAAAAAAAAAAGCCCCGGGGGGCGGTTTTGGTTGTAGAAAAG‑GGAAAAACCACCTTaaaaaaaaa                                                                                                 >  1:529843/1‑131 (MQ=255)                                                                     tgtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGct                                                                                <  2:1099853/150‑1 (MQ=255)                                                                      gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGcta                                                                               <  2:43716/150‑1 (MQ=255)                                                                        gACATTAATCAACCCCCAAAGATGTAACCATTGGCATTGGTATCGGCGTGCATGGCTAAAGTACGTAAAACCAGATCGCCCTGACGGACGTTATGAGTTGTAG‑ACAT‑GGTAAAACCGCCTTAAAATGAAAATCACAGGCGCGCTGcta                                                                               <  1:277394/148‑1 (MQ=255)                                                                          cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATACGGAAAAAAAAAAAGCCCCGAGGGGCGGTTTGTGTTTTTG‑ACAA‑GGGAAAAACCACTTAAAAAGAAAAAAACAAGCGCGAGGGTTC‑TAt                                                                          >  2:561027/1‑150 (MQ=255)                                                                           aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCCGAGGGGAGGTTTTTGTTTTTG‑ACAA‑GGGAAAAAAGCCTTTAAAAGAAAATAACAAGGGCGATg                                                                                  >  2:366182/1‑142 (MQ=255)                                                                           aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGCTAAAATAAGTAAAAAAAAAACGCCCCGAGGGGGGGTATGTGTTTTAGAAAAG‑GGAAAAAACACCTTAAAAGAaaaaaacaagggcga                                                                                     >  2:785467/1‑129 (MQ=255)                                                                              aaccaaccACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGAttt                                                                     <  1:348241/151‑1 (MQ=255)                                                                                 caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGCAGGCTTAAAAAAGTAAAACAAAAAAGCCCCGCGGGGGGGTTTGTGGTGTAT‑AAAT‑GGGAAAATCCACGTAAAAGGAAAAACAAAAGGGCGCTGCCAG‑TAATAtttt                                                                    >  1:887212/1‑149 (MQ=255)                                                                                  aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAAGACGTAAAAAAAAATCGCCCCGCGGGGGGGTATGTGTTGTTG‑ACAA‑GGGAAAAACGACTTAAAAAGAAAAACAAAGGCGCGATGCTAC‑TATGAttttt                                                                   >  2:869315/1‑148 (MQ=255)                                                                                       ccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCGTGGGTAAAATACGTAAAACAAAATCGGCCTGAGGGGGGTTGTGTGTTGTTG‑ACAT‑GGGAAAAACAACCTAAAAAGAAAATCACAGGGGCGCTGCTAC‑TATGATTTTTCttggg                                                             >  1:766664/1‑146 (MQ=255)                                                                                        caAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAAAACGGAAAAAAACAAAGGCCCGAGGGGGGGTAATGGGTGTGG‑AGAT‑GGGAAAAAAAACTTAAAAAGAAAAAAACAAGGGCGATGATACATATAAttt                                                                     >  2:328253/1‑142 (MQ=255)                                                                                        caAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTgaggag                                                           <  2:350942/151‑1 (MQ=255)                                                                                               gTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACagag                                                    <  1:101395/151‑1 (MQ=255)                                                                                               gTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACagag                                                    <  2:426684/151‑1 (MQ=255)                                                                                                  aCCATTGGCATTGGTATCGGCGGGCAAGGCTAAAGAAAATAAAACAAAAACGACCCGGGGGGCGGTTTGTGTTGTAG‑ACAA‑GGGAAAAACCACCTAAAAAAAAAAACACAGGCGGGATGCTAC‑TATTATTTTTCTGGGGGGaaaaaa                                                     >  1:1062435/1‑143 (MQ=255)                                                                                                        ggcattggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACa                                            <  2:1144554/150‑1 (MQ=255)                                                                                                         gcattggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAgc                                          <  2:1149116/151‑1 (MQ=255)                                                                                                           attggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAgcgc                                        <  1:386448/151‑1 (MQ=255)                                                                                                              ggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCa                                      <  1:403160/150‑1 (MQ=255)                                                                                                                     gcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGa                               >  2:1252943/1‑150 (MQ=255)                                                                                                                       gggCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAgg                            <  1:1283608/151‑1 (MQ=255)                                                                                                                           aTGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTg                        >  2:252214/1‑151 (MQ=255)                                                                                                                           aTGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGAGTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAg                             <  1:858946/146‑1 (MQ=255)                                                                                                                                tAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTTTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGt                   >  2:800049/1‑151 (MQ=255)                                                                                                                                 aaaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTc                  <  2:1248160/151‑1 (MQ=255)                                                                                                                                 aaaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTc                  >  1:1275353/1‑151 (MQ=255)                                                                                                                                       cGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGtt             >  1:1111593/1‑150 (MQ=255)                                                                                                                                       cGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGtt             >  1:1016753/1‑150 (MQ=255)                                                                                                                                        gTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGtt             >  1:441937/1‑149 (MQ=255)                                                                                                                                                aGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATc   <  2:259804/151‑1 (MQ=255)                                                                                                                                                aGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATc   <  2:915428/151‑1 (MQ=255)                                                                                                                                                 gATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCt  <  2:1003170/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                          CAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAG‑ACAT‑GGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTAC‑TATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCT  >  NC_000913/1312069‑1312363

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