Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,061	0	G	A	50.9%	49.6 / 97.9	57	intergenic (‑124/+79)	ves/spy	cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (27/1);  new base A (0/29);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 2.00e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-10

ATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAT  >  NC_000913/1824912‑1825210                                                                                                                                                        |                                                                                                                                                      aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                                      <  1:268197/151‑1 (MQ=255)     tACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAgaaaagaaaa                                                                                                                                                   <  2:537412/150‑1 (MQ=255)      aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                                 >  1:136657/1‑151 (MQ=255)       cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                                >  2:438503/1‑151 (MQ=255)                  gTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                                                                                                     >  1:508820/1‑151 (MQ=255)                      tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                                 >  1:626306/1‑151 (MQ=255)                      tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                                 >  1:574447/1‑151 (MQ=255)                          tACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAt                                                                                                                             >  2:103406/1‑151 (MQ=255)                                cccccGCCCCTTTTTTTTCATTGTTATCGTCC‑AGGTAACGGTGGATGGTTATCGGCAACTGGGGTGATTTAC‑GTTAGTGGTGATCAGGTTTTTTTTACACCTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                            <  1:771953/144‑1 (MQ=255)                                cTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                       >  2:759747/1‑151 (MQ=255)                                 tCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                       >  1:175563/1‑150 (MQ=255)                                 tCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                       >  1:726920/1‑150 (MQ=255)                                   ccGGCCCTTTTTTTTTATTTTTTACCTTCGAGGTACCGGTGGTTGTTTATCTG‑AACCTGGGTTATTTTCGGTTAGGGGGGATGGGTTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                                                                                     <  2:724978/151‑1 (MQ=255)                                       cccTTTTTTTTCATTTTTATCGTTC‑AGGTACCCGTGGATGTTTACCGGCAACTGGGGTTATTTAC‑GTTGGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                                                                                <  1:88315/151‑1 (MQ=255)                                        ccGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                               >  1:195768/1‑151 (MQ=255)                                         cGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                              >  2:334015/1‑151 (MQ=255)                                           ttCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                             >  1:549478/1‑150 (MQ=255)                                              ttttCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                          >  2:469072/1‑150 (MQ=255)                                               tttCATTTTTATCGTTA‑AGGTACCCGTGGTTGTTTAGCGGCACCTCGGGTTATTTCC‑GTTAGTGGTGATCAGTTTTTTTTTACAACTGTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                                                                                         <  2:397288/150‑1 (MQ=255)                                                ttCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                          >  1:673190/1‑148 (MQ=255)                                                   aTTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                    >  2:415877/1‑151 (MQ=255)                                                    ttGTTATCGTTC‑AGGTACCCGTGGTTGTTTATCGGCAACTGGGGTGATTTCC‑GTTAGTGGTGATCAGGTTTTTTTTACAATTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                    <  1:321196/150‑1 (MQ=255)                                                    ttGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                                   >  1:632074/1‑151 (MQ=255)                                                    ttGTTATCGTTA‑AGGTAACCGTGGATGGTTACCGGCAACTCGGGTGATTTAC‑GTTGGTGGGGATCAGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAc                                                                                                      <  1:946876/148‑1 (MQ=255)                                                                     acccGGGGATGGTTATCGGCACCTGGGGTGTTTTAC‑GTTAGGGGTGATCAGTTTTTTTTTACAACTGTTAGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt                                                                                    <  2:910836/148‑1 (MQ=255)                                                                     aaCCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                   >  2:736/1‑151 (MQ=255)                                                                         gTGGATGTTTATCGGCACCTGGGGTAATTTAC‑TTTAGTGGTGATCAGTTTTTTTTTACATTTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt                                                                               <  1:911204/151‑1 (MQ=255)                                                                              tGGTTACGGGCATCTGGGGTGATTTAC‑GTTAGTGGTGATCGGGTTTTTTTTAAATGTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa                                                                           <  1:76562/150‑1 (MQ=255)                                                                                      ggCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTTGGTGAAAGattat                                                                     >  2:641529/1‑148 (MQ=255)                                                                                        cAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                                 >  1:778969/1‑150 (MQ=255)                                                                                         accgcGGGTGTTTAAC‑TTTAGTGGTGATCGGTTTTTTTTTACACCTTTTTGAGAAAAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                       <  2:717447/139‑1 (MQ=255)                                                                                         aaCTGGGGTGTTTTAC‑TTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                                 <  1:432983/149‑1 (MQ=255)                                                                                              gggTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGtt                                                            >  1:394085/1‑149 (MQ=255)                                                                                                 tGATTTAC‑GTTAGTGGTGATCAGGTTTTTTTTAAAATTGTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                        <  2:858576/150‑1 (MQ=255)                                                                                                    tttAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTTCAGGCa                                                     >  1:581312/1‑150 (MQ=255)                                                                                                      tAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCAttt                                                  >  2:841411/1‑151 (MQ=255)                                                                                                           ttAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                                <  2:763084/149‑1 (MQ=255)                                                                                                           ttAGTGGTGATCAGTTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                                <  1:594621/149‑1 (MQ=255)                                                                                                           ttAGTGGTAATCGGGTTTTTTTTACAAGTGTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAc                                               <  1:407002/150‑1 (MQ=255)                                                                                                              gtggtgATCAGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc                                              <  2:358216/148‑1 (MQ=255)                                                                                                                ggtgATCGGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                         <  1:475249/151‑1 (MQ=255)                                                                                                                ggtgATCAGGTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                         <  2:867120/151‑1 (MQ=255)                                                                                                                   gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc                                       >  1:337080/1‑150 (MQ=255)                                                                                                                     tCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTTCAGGCATTTTACCTTTTGCCGc                                     >  1:822695/1‑150 (MQ=255)                                                                                                                      cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                                   >  1:932171/1‑151 (MQ=255)                                                                                                                       gggttttttttACACTTTTTTGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                                  <  2:808123/147‑1 (MQ=255)                                                                                                                           tttttttACATCTGTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg                               <  2:678672/150‑1 (MQ=255)                                                                                                                           tttttttACAACTGTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAc                                <  1:393071/149‑1 (MQ=255)                                                                                                                            ttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                             >  1:940364/1‑151 (MQ=255)                                                                                                                            ttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                             >  1:750457/1‑151 (MQ=255)                                                                                                                              ttttACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTc                            <  2:893895/150‑1 (MQ=255)                                                                                                                                     aCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa                     <  2:136656/150‑1 (MQ=255)                                                                                                                                      cTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                    <  1:914473/150‑1 (MQ=255)                                                                                                                                         tttGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCt                 <  1:22058/150‑1 (MQ=255)                                                                                                                                          tggagaagaaaaAAGAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt                <  1:332088/148‑1 (MQ=255)                                                                                                                                             aaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtca              <  2:989334/147‑2 (MQ=255)                                                                                                                                               aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc               <  1:137393/146‑1 (MQ=255)                                                                                                                                               aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaa           <  2:243769/150‑1 (MQ=255)                                                                                                                                                agaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa         <  1:110073/151‑1 (MQ=255)                                                                                                                                                    aaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg     <  1:456551/151‑1 (MQ=255)                                                                                                                                                       agaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAt  >  2:501451/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      ATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAT  >  NC_000913/1824912‑1825210

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