Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	42.3%	102.9 / 51.5	78	intergenic (‑108/+33)	uspF/ompN	stress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (40/5);  new base A (0/33);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 4.60e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATG  >  NC_000913/1435583‑1435877                                                                                                                                                    |                                                                                                                                                        cGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGaaaaa                                                                                                                                                    <  2:396065/150‑1 (MQ=255)    aaaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑ccc                                                                                                                                               <  2:335418/151‑1 (MQ=255)        atCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGcaa                                                                                                                                           <  2:902771/151‑3 (MQ=255)        atCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCa                                                                                                                                           <  1:534262/151‑1 (MQ=255)           gATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATa                                                                                                                                         >  2:1017900/1‑150 (MQ=255)           gATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATa                                                                                                                                         >  2:962735/1‑150 (MQ=255)                  aaaaGAATCTTTCTTTTAAAAAACCCTCTTTTTTTAGTACCCGCATAAAGTAAAACGCC‑AGATGACTTTTTTTGTGAAATGGCGGATTTCACTTTTTTAATTA‑AGAAAAAA‑GGAGAAAAAAAAAGCAAAAAAAG‑CCCGCCATAACAGCggg                                                                                                                                 <  1:622135/149‑1 (MQ=255)                   cAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGc                                                                                                                                >  1:359099/1‑151 (MQ=255)                       attcttttttttCATAAACCCTCTTTTTTTGGTATCCGAATAAAGTGAACCCCC‑AGATGACTTTTTTTGTGTAAGGAGGGATTTCACTTTTTTAATTT‑AGACAAAA‑GGAGAAAGAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCagg                                                                                                                             <  2:673978/140‑1 (MQ=255)                       attctttCTTTTCAAAAACCCTCTTTTTTTAGTTCCCGAAAAAAGTGAACCGCA‑AAATGACACTTTTTGTATAATGAGGGAGTTCAATTTTTTAATTA‑AGATAAAG‑GAAGAAAAAATAAAAAAAAAAGG‑CCCGCCATAACAGCGGGCagga                                                                                                                            <  1:19351/146‑1 (MQ=255)                       aaTCTTTCTGTTCAAAAACCCTCTTTTTTTATTACCCAAAAAAAGTGTAACCAC‑AGATGACATTTTTTGTAAAATGAGGGATTTCACATTTTTAATTA‑AGAAAAAA‑GGAGAAAAAAAAAAAAAAAAAAG‑CCCGCCATAACAGCGGGCagga                                                                                                                            <  1:357481/151‑1 (MQ=255)                           gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                        >  1:475818/1‑151 (MQ=255)                           gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                        >  2:380541/1‑151 (MQ=255)                            ttCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGAtt                                                                                                                       >  1:153389/1‑151 (MQ=255)                                gTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                                      >  1:1023785/1‑148 (MQ=255)                                 ttCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                                   >  1:382375/1‑150 (MQ=255)                                 ttCAAAAACCCTCCTTTTTTATTACCCGAAAAAAGTGAACCGCC‑AGATGACTTTTTTTGTGTAATGAGGGATTTCACTTTTTTAATTT‑AGATAAAA‑GGAGAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                                   <  2:782268/150‑1 (MQ=255)                                 ttCAAAAACCCTCCTTTTTTAGTCCCCGAATAAAGTGTACCGCC‑AGATGACTTTTTTTGTGTAATGACGGATTTCATTTTTTTAATTA‑AGAAAAAA‑GGGGAAAGAATAAGCAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAg                                                                                                                    <  2:212252/149‑1 (MQ=255)                                    aaaaacCCTCTTTTTTTATTACCCGAATAAAGT‑TAACCCCAAGATGACTTTTTTTGTGAAAGGAGGGATTTCA‑ATTTTTAATTTAAGAGAAAG‑GGGGAAAAAAAAAAAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                                   <  1:850862/145‑1 (MQ=255)                                    aTAAACCCTCTTTTTTTAGTATCCGCATAAAGTGAAACCCC‑AGATGACTTTTTTTGTGTAATGACGGATTTCACTTTTTTAATTT‑AGAAAAAA‑GGAGAAAGAATAAGAAAAAAAGG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                               <  2:573927/151‑1 (MQ=255)                                    aTAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                               >  1:319388/1‑151 (MQ=255)                                      aaacCCTCTTTTTTTAGTACCCGAATAAAGTGTAACGCC‑AGAGGACATTTTTTGTGTAATGACGGATTCCACATTTTTAATTA‑AGAAAAAA‑GGAGAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                              <  1:795521/150‑1 (MQ=255)                                      aaacCCTCTTGTTTTTGTACCCACAAAATGGGAAACGCC‑AGATGACTTTTTTTGTATAATGACGGAGTTCACATTTTAATTTT‑AGTAAAAA‑GGGGAAAAAATAAAAAAAAAAGC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                              <  1:873909/150‑1 (MQ=255)                                      aaaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                              >  1:241409/1‑150 (MQ=255)                                          cctcctTTTTTAGTACCCGCAAAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGAGGGAGTTCACATTTTTAATTT‑AGAAAAAA‑GGAGAAAAAATAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa                                                                                                          <  2:277390/150‑1 (MQ=255)                                          cctcctGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAc                                                                                                         >  2:472829/1‑151 (MQ=255)                                          cctcctGTTTTAGTATCCGCATAAAGAGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAc                                                                                                         >  1:1005148/1‑151 (MQ=255)                                             cctGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                                                                        >  1:858048/1‑149 (MQ=255)                                               ttttttAGTACCCGCAAAAATTGAAACCCC‑AGATACCATTTTTTGTGAAATGAGGGAGTTCAAATTTTTAATTA‑AGAAAAAG‑GGAGAAAAAAAAAAAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCa                                                                                                       <  1:926959/146‑1 (MQ=255)                                               ttttttAGTACCCGCAAAAAGTGAAACGCA‑AGATAACATTTTTTGTGTAAGGACGGATTTCACATTTTTAATTA‑AGATAAAA‑GGAGAAAAAATAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCa                                                                                                       <  1:478986/146‑1 (MQ=255)                                                 ttttAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                                  >  2:91749/1‑151 (MQ=255)                                                  tttAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                                  >  2:420807/1‑150 (MQ=255)                                                   ttAGTATCCGAATAAAGTGAAACGCC‑AGATGACACTTTTTGTGTAATGAGGGAGTTCATATTTTTAATTT‑AGAAAAAA‑GGAGGAAAAATAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTa                                                                                                 <  2:554634/150‑1 (MQ=255)                                                   ttAGTACCCGCAAAAAGTGAACCCCC‑ATATTACTTTTTTTGTTAAATGACGGATTTCACATTTTTAATTT‑AAATCAAA‑GGAGAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTa                                                                                                 <  2:23824/150‑1 (MQ=255)                                                     aGTACCCGCATAAAGTGTAACGCC‑AGATTACTTTTTTTGTGTAATGAGGGATTTCACTTTTTTAATTT‑AGAAAAAG‑GGAGGAAAAATAAGAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa                                                                                               <  2:430329/150‑1 (MQ=255)                                                          ccgaaaaaaGTGTAACGCC‑AGATGACTTTTTTTGTGTAATGAGGGATTTCACATTTTTAATTA‑AGAAAAAG‑GGAGAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGa                                                                                           <  2:973939/143‑1 (MQ=255)                                                          ccGCATAAAGTGTAACGCC‑AGATGACTTTTTTTGTGAAATGACGGAGTTCATATTTTAAATTT‑AGAAAAAA‑GGAGAAAAAAAAAGAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGa                                                                                           <  2:437598/149‑1 (MQ=255)                                                              aTAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                                     >  2:557648/1‑151 (MQ=255)                                                              aTAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                                     >  1:481287/1‑151 (MQ=255)                                                                aaaTTGTAACGCC‑AAATGACACTTTTTTTATAATGCCGGTTTCCCTTTTTTTAATTA‑AGATAAAA‑GGAGGAAGAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatc                                                                                    <  2:447793/150‑1 (MQ=255)                                                                  aGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcatc                                                                                 >  1:200928/1‑151 (MQ=255)                                                                    tgtAACCCC‑AGATGACACTTTTTGGATAATGGCGGATTCCACTTTTTTAATTT‑AGAAAAAA‑GGAGAAAAAAAAAGCAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                                                <  2:951209/150‑1 (MQ=255)                                                                       aaCGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                              >  2:383374/1‑149 (MQ=255)                                                                           cc‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                            >  1:581595/1‑147 (MQ=255)                                                                                  gactttttttGTGTAAGGACGGAGTTCACATTTTTAATTT‑AGATAAAA‑GGAGGAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCa                                                                    <  2:678587/144‑1 (MQ=255)                                                                                    cattttttGTGTAAGGAGGGATTTCAATTTTTTATTTT‑AGAAAAAG‑GGAGAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAtt                                                                  <  2:576135/146‑1 (MQ=255)                                                                                       tttttGTGTAAGGCCGGATTTCA‑TTTTTTAATTT‑AAAAAAAAGGGAGAAAGAAAAAAAAAAAAACG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                                <  1:391587/148‑1 (MQ=255)                                                                                          ttGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCg                                                           >  2:669756/1‑150 (MQ=255)                                                                                          ttGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGA‑TTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAg                                                             >  1:428103/1‑147 (MQ=255)                                                                                             tgaaaTGACGGAGTTAACATTTTTAATTA‑AGAAAAAA‑GGAGAAAAAAAAAAAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                                          <  1:1021924/145‑1 (MQ=255)                                                                                              gtAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                                       >  1:715540/1‑150 (MQ=255)                                                                                               tAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                                       >  2:143082/1‑149 (MQ=255)                                                                                                aaggACGGTGTTAACTTTTTTAATTT‑AGAAAAAG‑GGAGGAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAg                                                        <  1:1017754/144‑1 (MQ=255)                                                                                                aaggACGGAGTAATCATTTTTAATTT‑AAAAAAAA‑GGAGAAAAAAAAAAAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa                                                         <  1:230540/143‑1 (MQ=255)                                                                                                aaTGACGTAGTTCACTTTTTTATTTA‑AAAAAAAG‑GGAGAAAAAAAAAGAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAg                                                     <  1:111629/150‑1 (MQ=255)                                                                                                aaTGACGGATTCAACTTTTTAAATTA‑AGAAAAAA‑GGAGAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGaa                                                      <  1:1050952/149‑1 (MQ=255)                                                                                                aaTGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAg                                                     >  2:783956/1‑150 (MQ=255)                                                                                                      ggggtccaaaTTTTTAATTA‑AAAGAAAA‑GGAGAAAAAAAAAAAAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTg                                                  <  2:166526/138‑1 (MQ=255)                                                                                                      ggATTTCACTTTTTTAATTA‑AGAAAAAA‑GGAGGAAAAAAAAAAAAAAAACG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCAt                                              <  1:983378/151‑1 (MQ=255)                                                                                                       gAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATc                                             >  2:862133/1‑151 (MQ=255)                                                                                                        aGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                            >  2:944949/1‑151 (MQ=255)                                                                                                           tCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc                                          >  2:966157/1‑150 (MQ=255)                                                                                                                 ttttAATTTAAGATAAAA‑GGAGGAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa                                     <  2:409369/150‑1 (MQ=255)                                                                                                                  tttAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                   >  1:982802/1‑150 (MQ=255)                                                                                                                    tAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTTATCCAACAg                                 >  1:383219/1‑150 (MQ=255)                                                                                                                    tAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                                >  1:92845/1‑151 (MQ=255)                                                                                                                       ttt‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAggat                              >  1:829532/1‑148 (MQ=255)                                                                                                                       ttt‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                             >  1:499935/1‑151 (MQ=255)                                                                                                                       atttaaaaaaaa‑GGAAGAAAAAAAAAAAAAAAAACCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                               <  1:643265/142‑1 (MQ=255)                                                                                                                           agaaaaaa‑GGAGAAAAAAAAAAAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAt                           <  1:962596/145‑1 (MQ=255)                                                                                                                           aGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtt                          >  1:48228/1‑151 (MQ=255)                                                                                                                             aTCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                        >  1:830707/1‑151 (MQ=255)                                                                                                                              tCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                        >  2:395671/1‑150 (MQ=255)                                                                                                                                 aa‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                      >  1:956649/1‑149 (MQ=255)                                                                                                                                      agaaaaaaTAAGAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCa                 <  1:557547/144‑1 (MQ=255)                                                                                                                                       ggaAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAAc               >  1:82240/1‑151 (MQ=255)                                                                                                                                         aagaaaaagaaaAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAAc               <  2:1003247/149‑1 (MQ=255)                                                                                                                                            aaTAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAgg          >  1:413181/1‑151 (MQ=255)                                                                                                                                                aGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTgg       >  2:854098/1‑150 (MQ=255)                                                                                                                                                    gAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATg  >  2:782377/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCC‑AGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAA‑GGAGGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATG  >  NC_000913/1435583‑1435877

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