Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	39.2%	109.8 / 55.4	79	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/8);  new base A (0/31);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 7.02e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAT  >  NC_000913/1435579‑1435876                                                                                                                                                        |                                                                                                                                                       gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGaa                                                                                                                                                      <  1:587738/151‑1 (MQ=255) gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGaa                                                                                                                                                      <  1:1179063/151‑1 (MQ=255) gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGaa                                                                                                                                                      <  2:1054787/151‑1 (MQ=255) gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGa                                                                                                                                                       >  1:454776/1‑150 (MQ=255)    tCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGaaaa                                                                                                                                                    >  1:399854/1‑150 (MQ=255)       gAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑Gcc                                                                                                                                               <  2:127205/151‑1 (MQ=255)       gAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑Gcc                                                                                                                                               <  1:723864/151‑1 (MQ=255)       gAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑Gcc                                                                                                                                               <  2:460684/151‑1 (MQ=255)          atatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGc                                                                                                                                            <  1:687152/151‑1 (MQ=255)           tatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGcc                                                                                                                                           >  2:560382/1‑151 (MQ=255)                      aCAAGATTCGTTTTTTTCATAAACCCCCCTTTTTTGGTATCCGAATAAATTGTACCCCCGTATTCCTTTTTTTGTA‑TAATGAGGGCGTTCACTTTTTTATTTG‑AAATAAAAGGAGA‑AAAAATAAGCAAAAAAA‑GCCCGCCATAACAGCgg                                                                                                                                 <  2:810505/150‑1 (MQ=255)                       cAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGc                                                                                                                               >  2:800996/1‑151 (MQ=255)                         aGAATCGTTTTGTTCAAAAACCCTCCTGTTTTAGGACCCGAAAAAAGTAAACCGCCAGATGACATTTTTTGTA‑AAAAGAGGGAGTTCAAATTTTTAATTA‑AGAAAAAAGGAGA‑AAAAAAAAGAAAAAAAA‑CCCCGCCATAACAGCGGGCag                                                                                                                             <  1:496842/151‑1 (MQ=255)                           aaTCGTTTTGTACAAAAACCCTCTTGTTTTAGTATCCACATAAAGTGAAACGCCAGATGACTTTTTTTGTA‑TAATGACGGATTTCAATTTTTAAATTA‑AGAAAAAGGGAGA‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCagga                                                                                                                           <  1:1113796/151‑1 (MQ=255)                              cGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCaggagg                                                                                                                         >  2:52839/1‑150 (MQ=255)                               gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                       >  2:1231668/1‑151 (MQ=255)                                  cTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                                     >  2:1037037/1‑150 (MQ=255)                                   tGTACACAAACCCTCTTGTTTAAGTTCCCAAATAAAGTGAAACGCCAGATGTCTTTTTTTGTG‑GAAGCGAGGATTTCACATTTTTATTTT‑AGA‑CGAAGGAGGAAAAAAAAAGAAGAAAAA‑GCCCGCCATAACAGCGGGCAggagga                                                                                                                        <  1:412781/146‑1 (MQ=255)                                     ttaaaaaacCCCCCTTTTTTAGTATCCGCAAAAAGTGTACCGCCAGATGACATTTTTTGTG‑TAAGGAGGGAGTTCAAATTTTTATTTT‑AAATAAAAGGAGG‑AAAAATAAGCAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                                     <  2:1061431/142‑1 (MQ=255)                                        aaaaacCCTCTTGTTTTAGTACCCGAAAAAAGTTAAACCCAAGATAACTTTTTTGGTG‑TAATGGCGGATTTCATTTTTTAATTTT‑AAAAAAAGGGAGA‑AAAAATAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACt                                                                                                               <  1:1252878/148‑1 (MQ=255)                                        aTAAACCCTCCTGTTTTAGTACCCGCAAAAATTGTAACGCAAGATAACATTTTTTTTT‑TAATGCCGAATTTCAATTTTTTAATTA‑AAACAAAAGGAGA‑AAAAAAAAGCAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAAc                                                                                                                <  2:600629/149‑1 (MQ=255)                                        aTAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                              >  1:455852/1‑151 (MQ=255)                                        aTAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                              >  1:870475/1‑151 (MQ=255)                                        aTAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                              >  2:454544/1‑151 (MQ=255)                                         tAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                              >  2:153624/1‑150 (MQ=255)                                          aaacCCTCTTTTTTTAGTACCCGAAAAAGTTGTACCCCCAGATAACTTTTTTTTTG‑AAATGAGGGATTTCATTTTTTTATTTA‑AAAAAAAAGGAGA‑AAAAATAAGCAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGaa                                                                                                                 <  2:431976/146‑1 (MQ=255)                                          aaacCCTCCTTTTTTAGTACCCGAAAAAAGTGTAACCCAAGATGACTTTTTTTGTG‑AAAGGGCGGAGTTCACTTTTTTAATTT‑AGATAAAAGGAGA‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                             <  2:979348/150‑1 (MQ=255)                                           aaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                            >  2:287774/1‑150 (MQ=255)                                              cccctttttttAGTACCCGAAAAAAGTGAAACCCCAGATAACTTTTTTTGTT‑AAAGGACGAATTTCATATTTTTAATTT‑AAAAAAAAGGAGA‑AAAAATAAGCAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaa                                                                                                          <  1:1186867/142‑1 (MQ=255)                                                tcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                                                                       >  2:852312/1‑150 (MQ=255)                                                    gTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                                                  >  1:1175394/1‑151 (MQ=255)                                                      tttAGTCTCCAAAAAAAGTGTAACCCCAAAAGACATTTTTGGTG‑AAAGGGGGGAGTTCATTTTTTTAATTT‑AGATAAAAGGGGG‑AAAAAGAAGCAAAAAAC‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                                 <  2:74201/150‑1 (MQ=255)                                                      tttAGTACCCGAAAAAAGTGTACCGCCAGATGACTTTTTTTGTG‑AAATGACGGAGTCCACATTTTTAATTA‑AGAAAAAAGGAGA‑AAGAAAAAGAAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                                 <  2:1156260/150‑1 (MQ=255)                                                                gCAAAAAGTGAAACGCGAAATAACTCTTTTTGTG‑AAAGTACGGCATTAACATTTTTAATTT‑GGACAAAAGGGGA‑AAAAATAACAAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTa                                                                                                <  1:654488/141‑1 (MQ=255)                                                                  aaaaaGTGTACCGCAAAATAACATTTTTTGTG‑AAAGGAGGGATTTCA‑ATTTTTAATTTAAAAAAAAAGGAGA‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                                    <  2:109746/149‑1 (MQ=255)                                                                  aTAAAGTGAAACGAAAAATGAC‑TTTTTTGTGAAAATAAAGGAGTTCACTTTTTTATTTA‑AAATAAAAGGAAA‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAata                                                                                     <  1:254359/150‑1 (MQ=255)                                                                        tgtAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCgg                                                                              >  1:2710/1‑151 (MQ=255)                                                                         gtAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                             >  2:468182/1‑151 (MQ=255)                                                                            aCGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                             >  1:862912/1‑148 (MQ=255)                                                                               ccAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                           >  2:794088/1‑147 (MQ=255)                                                                                cAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AATAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                                      >  2:46495/1‑151 (MQ=255)                                                                                   aTGACTTTTTTTGTG‑TAAGGGCGGAGTTCATATTTTTAATTA‑AGAAAAAAGGAGA‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGcc                                                                    <  1:449899/150‑1 (MQ=255)                                                                                     gacttttttttTA‑TAATGGCGCGTTTCCCTTTTTAAATTA‑AAATAAAGGGAGG‑AAAAATAAGAAAAAAAG‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAgaga                                                                        <  2:1109549/139‑1 (MQ=255)                                                                                     gACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAtt                                                                 >  2:1084163/1‑151 (MQ=255)                                                                                         ttttttGTG‑AAATGCGGGATTTCAATTTTTTAATTA‑AAAAAAAGGGAGA‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                               <  1:1029073/148‑1 (MQ=255)                                                                                          tttttGTG‑TAAATCAGGGGTTACCATTTTTAATTA‑AGAACAAAGGAGA‑AGGAATAAGAAAAAAAG‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                              <  2:1048992/149‑1 (MQ=255)                                                                                            ttgttg‑tAATGACGGATTTCACTTTTTTATTTT‑AGAAAAAGGGAGG‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                                           <  2:759471/146‑1 (MQ=255)                                                                                              ttgttaaaTGAGGGAGTTCATTTTTTAAATTA‑AAAAAAAGGGAGA‑AAAAAAAAAAAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCg                                                          <  1:112884/144‑1 (MQ=255)                                                                                                 g‑tAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAg                                                       >  1:681844/1‑149 (MQ=255)                                                                                                   tAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                                      >  2:513319/1‑149 (MQ=255)                                                                                                   tAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                                      >  2:442489/1‑149 (MQ=255)                                                                                                    aaTGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAG‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAg                                                    >  1:578277/1‑150 (MQ=255)                                                                                                         cgttcttgatttttttAATTT‑GACTAAAGGGAGA‑AAAAAAAATAAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGaa                                                     <  1:150119/133‑1 (MQ=255)                                                                                                            aGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                           >  2:456460/1‑151 (MQ=255)                                                                                                            aGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACATAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                           >  1:700724/1‑151 (MQ=255)                                                                                                              ttaatttttttAATTA‑AAATAAAAGGAGG‑AAGAAAAAACAAAAAAG‑ACCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATc                                            <  1:993449/142‑1 (MQ=255)                                                                                                              ttCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc                                         >  2:389797/1‑151 (MQ=255)                                                                                                                cacaTTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                       >  2:942747/1‑151 (MQ=255)                                                                                                                  aaTTTTTATTTA‑ACAAAAAAGGGGA‑AAAAAAAAGAAAAAAAG‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                       <  1:415756/148‑1 (MQ=255)                                                                                                                   aTTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATcc                                     >  2:940445/1‑150 (MQ=255)                                                                                                                    tttttAATTA‑AGATAAAAGGAGA‑AAAAATAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa                                    <  2:957569/150‑1 (MQ=255)                                                                                                                        tAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                               >  2:1182805/1‑151 (MQ=255)                                                                                                                         aaTTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                              >  1:640668/1‑151 (MQ=255)                                                                                                                          aTTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGtt                             >  1:842963/1‑151 (MQ=255)                                                                                                                               agagaaaaGGAGG‑AAGAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtt                         <  2:744837/146‑1 (MQ=255)                                                                                                                               aaaaaaaaGGAGA‑AAAAATAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                            <  2:754093/143‑1 (MQ=255)                                                                                                                               aaaaaaaaGGAGA‑AAAAAAAAGAAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGa                           <  1:947077/144‑1 (MQ=255)                                                                                                                                gATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAttt                        >  2:872746/1‑151 (MQ=255)                                                                                                                                aaaaaaaGGAGA‑AAAAAAAAGAAAAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAttt                        <  1:150133/147‑1 (MQ=255)                                                                                                                                   aaaaGGGGA‑AAAAAAAAGAAAAAAAC‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                     <  2:1207314/150‑1 (MQ=255)                                                                                                                                   aaaaGGAGA‑AAAAAAAAGAAAAAAAA‑CCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGa                           <  1:444700/144‑1 (MQ=255)                                                                                                                                    aaaGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                     >  1:896402/1‑150 (MQ=255)                                                                                                                                     aaGGAAG‑AAAAAAAAACAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAt                          <  1:82450/145‑1 (MQ=255)                                                                                                                                         agg‑aAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCa                >  1:379481/1‑150 (MQ=255)                                                                                                                                             aagaaaaagaaaAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAAc              <  1:765800/149‑1 (MQ=255)                                                                                                                                                aaTAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAgg         >  2:1192594/1‑151 (MQ=255)                                                                                                                                                 aTAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAgg         >  1:869615/1‑150 (MQ=255)                                                                                                                                                 aTAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGt        >  2:1021418/1‑151 (MQ=255)                                                                                                                                                       aGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAt  >  2:1030077/1‑151 (MQ=255)                                                                                                                                                       aGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAt  >  2:975228/1‑151 (MQ=255)                                                                                                                                                        gAAAAA‑GCCCGTCATAACCGCGGGCAGACGGATTTAGAACTGATAAACCAGACCTAAACCGCCAATATCATCGGTAGGGATGCCATTGCAAGCGTAGACGCTGTCATCTTCATCCACCAGGGTGATTTTATAGTCAACGTAGGAGGatat  >  1:129070/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                       GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTG‑TAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGG‑AAGAATAAGCAGAAAAA‑GCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAT  >  NC_000913/1435579‑1435876

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