Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	61.0%	12.9 / 63.5	59	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (5/18);  new base A (36/0);  total (41/18)
Fisher's exact test for biased strand distribution p-value = 5.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-10

ATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTA  >  NC_000913/2386728‑2386997                                                                                                                                        |                                                                                                                                                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                                                  <  2:424895/150‑1 (MQ=255)                                                                         gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCt                                                  <  1:673827/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:905706/1‑148 (MQ=255)                                                                                                                           ggTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa  >  1:414426/1‑148 (MQ=255)                                                                                                                           ggTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTCCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa  >  2:785704/1‑148 (MQ=255)                                                                                                                                        |                                                                                                                                                ATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCC‑AGAAAAAGCGAC‑ATGC‑AC‑ATTTTATTTTCCCTCATC‑T‑G‑TTTTAAGATGC‑GAAACATGAAATTCAACGGAGGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTA  >  NC_000913/2386728‑2386997

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