Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,868	0	T	A	53.4%	31.8 / 86.2	60	intergenic (+37/‑66)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base T (10/17);  new base A (31/0);  total (42/18)
Fisher's exact test for biased strand distribution p-value = 4.27e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-07

ATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATA  >  NC_000913/2386728‑2387018                                                                                                                                               |                                                                                                                                                              ataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑atttta                                                                                                                                                   <  1:716743/151‑1 (MQ=255)   aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAAAGAC‑AAAC‑AC‑aatttttt                                                                                                                                                 >  1:475058/1‑151 (MQ=255)     aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑attttatttt                                                                                                                                               <  2:654771/151‑1 (MQ=255)      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATg                                                                                                                                                              <  1:415682/137‑1 (MQ=255)      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑attttatt                                                                                                                                                 <  1:888673/148‑1 (MQ=255)      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAAAGAC‑AATC‑AC‑atttttttttt                                                                                                                                              >  2:953640/1‑150 (MQ=255)      cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCCCGGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTAC‑AAAAAAAAGGAA‑ATGC‑AC‑attttattttt                                                                                                                                              >  2:888980/1‑150 (MQ=255)        aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCC‑AGAAAAAAAGAA‑AAGA‑AC‑aatttatttttt                                                                                                                                             >  2:154439/1‑148 (MQ=255)              cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAACGAC‑AAGG‑AC‑CATTTTTTTTTCCTCAt                                                                                                                                        >  1:246654/1‑149 (MQ=255)                    cGGTGATACCGCAGACAGATTGAATTATCAGAAAATTGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑AGGC‑AA‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTa                                                                                                                             <  2:7836/151‑1 (MQ=255)                      gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTC‑TGAAAAAATGAC‑AAGC‑AC‑CCTTTTTTTTTCCTCTTC‑T‑C‑TTTTaaa                                                                                                                           >  2:266929/1‑150 (MQ=255)                      gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAAAGAA‑AAGC‑AC‑CATTTTTTTTTCCTCCTC‑T‑G‑TTTTTAg                                                                                                                           >  2:300641/1‑151 (MQ=255)                      gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCC‑CGAAAAAAAGAC‑AAGC‑AA‑AATTTTTTTTTCCTCATC‑T‑T‑TTTTaaa                                                                                                                           >  2:391661/1‑150 (MQ=255)                      gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATAAAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGC‑AGAAAAAACGAA‑AAGC‑AC‑ATTTTATTTTTTCTCATC‑T‑G‑TTTTata                                                                                                                           >  1:761550/1‑149 (MQ=255)                      gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATAAAACAGTAAATAATAAAAATTTCTCGAGGCGGGATC‑CTAAAAAAAGAG‑ATGT‑AC‑AATTTTTTTTTCCttttttttttttaa                                                                                                                              >  2:415240/1‑138 (MQ=255)                                agacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGC‑AAAAAAAAAGAA‑AAGA‑AC‑ATTTTTTTTTTCCTCCTC‑T‑T‑TTTTTAGAAGC‑Ggaaaa                                                                                                                >  2:80135/1‑149 (MQ=255)                                  acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAAc                                                                                                                 <  2:918640/148‑1 (MQ=255)                                      aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAAt                                                                                                          <  1:241833/151‑1 (MQ=255)                                             tATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACgg                                                                                                   <  2:8600/151‑1 (MQ=255)                                                       tGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAaaat‑a                                                                                        <  1:501637/151‑1 (MQ=255)                                                        ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCC‑CGAAAAAAAGAC‑AAGC‑AA‑ATTTTTTTTTTTCTCCTC‑T‑G‑GTTTTAGATGA‑GAAAAAAGAAAATTAAAAGAGGGAaaaa‑aa                                                                                       >  1:313685/1‑151 (MQ=255)                                                         gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAaaat‑aa                                                                                       <  2:278462/150‑1 (MQ=255)                                                            cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTC‑CGAAAAAAAGAC‑AAAC‑AC‑ATTTTTTTTTTTCTCCTC‑T‑G‑TTTTTTTATTC‑GGAAAATGAAAATTAACGGAGGAAaaat‑aaaatg                                                                                   >  1:821706/1‑148 (MQ=255)                                                             aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTC‑AGAAAAAACGAA‑AAGC‑AA‑ATATTTTTTTTCCTCACC‑T‑G‑TTTTTAGAAGG‑TTAACATGAAAATCTAACGAGGAaaaaa‑aaaaagg                                                                                  >  2:918759/1‑147 (MQ=255)                                                              ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTC‑CGAAAAAAAAAA‑AAGC‑AC‑ACTTTTTTTTTCCTCTTC‑T‑GTTTTTAAAAAGG‑GAAA‑AAGAAAATTAACAGAGGGAaaat‑aaaatggg                                                                                 >  1:855056/1‑146 (MQ=255)                                                              ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTC‑CGAAAAAAAGAC‑AAGCAAC‑ATTTT‑TTTTTCCTTATC‑C‑G‑GTTTTAGATGCGGAAAAAAAAATTTAAAGGGGGGAAAAAA‑AAATggg                                                                                  >  1:28333/1‑151 (MQ=255)                                                               gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAAAGAA‑AAGC‑AC‑AATTTATTTTTTCTCATC‑T‑G‑TTTTAAGATGC‑GAAACAAGAAAATCAACAGAGGGAAAAT‑AAATGGGaa                                                                                >  1:691110/1‑151 (MQ=255)                                                                  gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑CGAAAAAAAGAA‑AAGA‑AA‑ATTTTATTTTTCCTCATC‑T‑G‑TTTTTAAATGA‑GAAAAAAAAAATTCAACAGGGGGAAAAT‑AAAAGGGGATGa                                                                             >  1:652822/1‑151 (MQ=255)                                                                    ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGC‑AGAAAAAAAGAC‑AAAC‑AA‑ATTTTATTTTTTCTCCTC‑T‑G‑TTTTTAAATGA‑GAAAAAAAAAAATTAAAGGAGGGAAAAT‑AAATAGGGAAGAtt                                                                           >  1:305600/1‑151 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTC‑CGAAAAAACGAC‑AAGC‑CCAATTTTTTTTTCCCCCTCT‑T‑T‑TTTTAAAATGG‑AAAAAAGAAAATTAAAAGAGGGAAAAAA‑AAATGGGAATGTTTggtt                                                                       >  2:166650/1‑149 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTC‑CGAAAAAAAGAC‑AAGA‑AA‑ATTTTTTTTTTCCTCATC‑T‑G‑TTTTTAGAAGC‑GAAAAAAGAAAATTAAACGAGGGAAAAT‑AAATTGGGATGAATGGCTa                                                                      >  1:403992/1‑151 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTC‑CGAAAAAAAGAA‑AAAC‑AC‑AATTTTTTTTTCCCCTTCTT‑T‑TTTTAAAATCG‑GAAAATAAAAATTAAACGAGGGAAAAAA‑AATTGGGAATgttgggta                                                                       >  2:1042354/1‑144 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑CGAAAAAAAGAC‑AAGC‑AC‑ATTTTTTTTTTCCTCATC‑T‑T‑TTTTTAAAAAC‑GAAAAAAGAAAATTAAAGGGGGAAAAAA‑AAAAAGGAATGATTGGGTa                                                                      >  1:344950/1‑151 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCt                                                                       <  1:673827/150‑1 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCt                                                                       <  2:424895/150‑1 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAACGAC‑AAGA‑AC‑ATTTTTTTTTTCCTCCTC‑T‑G‑TTTTTAGATGC‑GTAACAAAAAAATTAAAGGAGGAAAAAT‑AAATAGGGATGATTTGCtt                                                                      >  2:329827/1‑150 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTAC‑AGAAAAAAAGAC‑ATGC‑AC‑AATTTATTTTTCCCCATC‑T‑T‑TTTTTAGAAGC‑GAAAAAAGAAAATTAACGGAGGAAAAAT‑AAATGGGAATTATTGGCTa                                                                      >  2:363618/1‑151 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTC‑CGAAAAAAAAAC‑AAGC‑AC‑ATTTTTTTTTTCCCCCTC‑T‑G‑TTTTTAGATGC‑GGAAAAAAAAATTAAAAGGCGGGAAAAT‑AAATTGGAAAGaattggta                                                                      >  1:546963/1‑144 (MQ=255)                                                                          tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTa                                                                      <  1:920445/150‑1 (MQ=255)                                                                           gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAAACAAA‑AAGA‑AA‑ATTTTTTTTTTCC‑CATC‑T‑G‑TTTTTAAAAGCGGAAAAATAAAATTAAACGGGGGAAAAAT‑AAAGGGCAAAAATTGGcatt                                                                     >  1:591579/1‑148 (MQ=255)                                                                                    ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑CCAAAAAAGAAC‑AACC‑CA‑ATTTTTTTTTTCCTCCTC‑T‑T‑GTTTTAGATGG‑GGAACATGAAAATTAACGGGGGGAAAATAAAATGGGATGGATTGG‑TTTTTTTCAGaaa                                                           >  2:511788/1‑151 (MQ=255)                                                                                      gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCC‑CGAAAAAAAGAC‑AAGG‑AC‑AATTTTTTTTTCCTCCTT‑T‑T‑TTTTTAAAAGG‑GGAAAAAGAAAATAAACGGAGGAAAAAA‑AAAATGGGATTATTGGGTATTTTTTAGaaaat                                                         >  2:544933/1‑149 (MQ=255)                                                                                       ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCC‑AGAAAAAACGAC‑AAGC‑AC‑ATTTTTTTTTTCCTCATC‑T‑G‑TTTTTAGAAGC‑GAAAAATGAAAATCAAAGGAGGGAAAAA‑AAAAAGGAATGATTTGCTATTTTTTAGAAAATg                                                        >  2:116850/1‑151 (MQ=255)                                                                                         ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGC‑AGAAAAAAAGGG‑AAGC‑AA‑AATTTATTTTTCCTCATC‑TCG‑TTTTAAATGGG‑GAAAAAGAAAATTAAAAGGGGGAAAAA‑‑AAATAGGAATAATTGGTTCTTTTTCAGaaaataaa                                                      >  2:129691/1‑146 (MQ=255)                                                                                         ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGAc                                                      <  2:30262/151‑1 (MQ=255)                                                                                              acaGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAAAGAC‑AAGC‑AA‑ATTTTTTTTTTCCTCATC‑T‑G‑TTTTTAGATGC‑GAAAAAAGAAAATCAAAGGAGGGAAAAA‑AAATAGGGATGATTTGCTATTTTTGAGAAAATGAATCAGa                                                 >  1:857572/1‑151 (MQ=255)                                                                                              acaGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAAAGAA‑AAGG‑AA‑ATTTTATTTTTCCTCATC‑T‑G‑TTTTTAGAAGC‑GAAAAAAGAAATTTAACCGGGGAAAAAT‑AAAAGGGGATGATTTGGTATTTTTTAGAAAATTAATCAGa                                                 >  1:822664/1‑151 (MQ=255)                                                                                              acaGTAAATAATAAAAATTTCTCGAGGCGGGGTC‑AGAAAAAAAGAA‑AAGC‑AC‑ATTTTTTTTTTCCTCATC‑T‑G‑TTTTTAGATTG‑GAAACATGAAAATTAAAAGAGGGAAAAT‑AAATAGGGATGATTTGCTATTTTTGAGAAAATGACTCAGa                                                 >  1:930860/1‑151 (MQ=255)                                                                                                 gTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AAAAAAAACGGCAAGGC‑CC‑TTTTTTTTTTCCCCCTTT‑T‑T‑TTTTAA‑ATGT‑GGAACAAAAAAATTAACGGAGGAAAAAA‑AAATTGGAATGATTTTTTTTTTTTCAAAAAATGAATCaaaaaa                                              >  2:349249/1‑151 (MQ=255)                                                                                                      taataaAAATTTCTCGAGGCGGGTCC‑CGGAAAAAAGAC‑ATTG‑AC‑AATTTTTTTTTCCTCCTC‑C‑G‑TTTTTAAAAGC‑GAAAAAAGAAAATTAAAGGGGGGAAAAA‑AAATGGGGAAGATTTGGTTTTTTTTCGAAAATTACTTAGaaaaaaata                                         >  2:645978/1‑146 (MQ=255)                                                                                                         taaAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                                      <  1:154403/151‑1 (MQ=255)                                                                                                          aaaaaTTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                                     <  1:895699/151‑1 (MQ=255)                                                                                                                      aGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGt                          <  2:857863/150‑1 (MQ=255)                                                                                                                        gCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc                         <  2:657541/149‑1 (MQ=255)                                                                                                                           ggTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                       >  2:837266/1‑148 (MQ=255)                                                                                                                           ggTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                       >  2:88652/1‑148 (MQ=255)                                                                                                                           ggTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                       >  1:414426/1‑148 (MQ=255)                                                                                                                           ggTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                       >  1:905706/1‑148 (MQ=255)                                                                                                                           ggTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTCCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                       >  2:785704/1‑148 (MQ=255)                                                                                                                                            c‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGa    >  1:52105/1‑151 (MQ=255)                                                                                                                                               tGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGata  >  1:923235/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                              ATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATA  >  NC_000913/2386728‑2387018

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