Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,056	0	G	A	42.3%	73.7 / 24.4	52	intergenic (‑119/+84)	ves/spy	cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (27/3);  new base A (0/22);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 8.50e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAA  >  NC_000913/1824907‑1825203                                                                                                                                                        |                                                                                                                                                    ccGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagaga                                                                                                                                                    <  2:122299/151‑1 (MQ=255) ccGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagag                                                                                                                                                     <  1:578547/150‑1 (MQ=255)  cGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagagaa                                                                                                                                                   <  2:449564/151‑1 (MQ=255)      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:726920/1‑150 (MQ=255)                                      tCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                >  1:175563/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:750457/1‑151 (MQ=255)                                                                                                                                 ttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                      >  1:940364/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)                                                                                                                                                        |                                                                                                                                                    CCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTC‑AGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAA  >  NC_000913/1824907‑1825203

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: ‑