Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,642,689	0	T	G	43.1%	83.7 / 12.0	58	A52A (GCA→GCC)	rodZ	MreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base T (9/24);  new base G (25/0);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 3.10e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGT  >  NC_000913/2642544‑2642838                                                                                                                                                   |                                                                                                                                                          aaaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCtt                                                                                                                                                      <  1:1107771/151‑1 (MQ=255)  aaaaCTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTa                                                                                                                                                     <  1:103161/151‑1 (MQ=255)        ctGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTAtctt                                                                                                                                                 <  2:508108/149‑1 (MQ=255)          gCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTAtcttc                                                                                                                                                <  2:1049397/148‑1 (MQ=255)          gCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATcttctt                                                                                                                                              <  1:954194/150‑1 (MQ=255)            aTCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCtttttttttt                                                                                                                                              >  1:933691/1‑140 (MQ=255)             tCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGCCTTATTTTTTTCTAt                                                                                                                                          >  2:70212/1‑151 (MQ=255)                  gcAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTTTCTTCTTCTAt                                                                                                                                          >  2:1071228/1‑146 (MQ=255)                  gcAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTcgc                                                                                                                                     <  2:933662/151‑1 (MQ=255)                     aCTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑a                                                                                                                                <  2:833756/151‑1 (MQ=255)                          tGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGGCTTTTTTTTTTTTATGTCTC‑GT‑A‑c                                                                                                                              >  2:862446/1‑147 (MQ=255)                           gCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTgg                                                                                                                        <  2:390662/151‑1 (MQ=255)                           gCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTg                                                                                                                         <  2:610098/150‑1 (MQ=255)                           gCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTg                                                                                                                         <  2:946841/150‑1 (MQ=255)                           gCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCCTTTTTTCTTTTATTTTGC‑GT‑A‑CC‑CTgg                                                                                                                        >  1:909968/1‑151 (MQ=255)                             aGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCCTTTCTTTTTTTATGTCGC‑GG‑A‑CC‑GCgg                                                                                                                        >  1:190866/1‑149 (MQ=255)                              gCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTATTTTTTTTAATTTCGC‑GG‑A‑CC‑CTgg                                                                                                                        >  1:710268/1‑148 (MQ=255)                                   aaGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTTTCTTCTTTTATTTcgc‑ggtccgggggaa                                                                                                                        >  2:409970/1‑135 (MQ=255)                                      tGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTTTTTTTTTTTATGTCTC‑GT‑A‑CC‑CGGGGAAACTTCAgg                                                                                                             >  2:773459/1‑151 (MQ=255)                                      tGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGCCCTATTTTTTTCTATGtgccctagccgggg                                                                                                                           >  2:403303/1‑128 (MQ=255)                                       ggAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAgg                                                                                                             <  1:823063/150‑1 (MQ=255)                                        gAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCGTATCTTTTTCTATATCGC‑GT‑ACCC‑GGGGAAACCtcatg                                                                                                              >  2:626655/1‑145 (MQ=255)                                        gAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCCTTTTTTCTTTCATGTCGCGGA‑A‑CC‑TGGGAAACCtccgg                                                                                                              >  1:224976/1‑145 (MQ=255)                                              gCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTTTCTTTTTCAATTTCGC‑GT‑A‑CC‑GTGGGAACCTTCAGGCaaaaac                                                                                                      >  1:1033622/1‑147 (MQ=255)                                                ttttCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGCCTTTTTTTTTTCTATTTCGC‑CT‑A‑CC‑GGGGGAAACCTCAGGCaaaaaccc                                                                                                    >  1:541843/1‑145 (MQ=255)                                                ttttCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGCCTTATCTTCTTCAAAGTTGC‑GTGA‑CG‑GGGGAACCTTTCAGGCaa                                                                                                          >  1:310067/1‑145 (MQ=255)                                                 tttCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCt                                                                                                   <  2:633505/150‑1 (MQ=255)                                                 tttCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTTTTTTTTTTTATGTCGC‑GT‑G‑TC‑GTGGAGAACTTCTCGCaaaagcg                                                                                                     >  2:352170/1‑144 (MQ=255)                                                 tttCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTTTTTTTTTTAATGTTGC‑GT‑A‑AC‑CTGGGAAACTTCTGGCAAAAACGCTc                                                                                                  >  2:756279/1‑151 (MQ=255)                                                   tCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCg                                                                                                 <  1:740533/150‑1 (MQ=255)                                                       gCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAAc                                                                                            <  1:907010/151‑1 (MQ=255)                                                              cagcagTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCctgc                                                                                     <  2:967275/151‑1 (MQ=255)                                                                gcagTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTATCTTTTTCTATTTCGC‑GT‑A‑CC‑GCGGGAACCTTCTAGCAAAAACGCTCGGccaacgc                                                                                         >  1:597677/1‑139 (MQ=255)                                                                 cagTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGa                                                                                  <  2:779167/151‑1 (MQ=255)                                                                  agTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTTTTTTTTTTTAAGTCGC‑GC‑A‑CC‑GGGGGAAACCTCACGCAAAAACGCTCGGCAAccggccgccgac                                                                                 >  1:247688/1‑140 (MQ=255)                                                                              ggAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTTTTCAATGTCGC‑GG‑A‑CC‑GGGGGAAACCTCCGGCAAAAACGCTCGGCAACGGCCTGCCGGCTAAGTcccaatt                                                                     >  1:979483/1‑146 (MQ=255)                                                                                aaTATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTTTTTTTTTTCATTTCTC‑GG‑A‑AC‑CTGGGAAAACTCAGGGAAAGAAGCTCGGCAACGGCCCGCTGGCTAAGTCCTAggt                                                                     >  2:786353/1‑147 (MQ=255)                                                                                   atGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTc                                                                 <  2:1088103/150‑1 (MQ=255)                                                                                     gTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTATCTTTTTCTATTTCGC‑GG‑A‑CCGGGGGAAACCTT‑AGGCAAAAAAGCTCGGCAAAGGCCCGCTGACTCAATccccgttttt                                                                  >  2:638624/1‑139 (MQ=255)                                                                                        ccAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGCCCTTTTTTCTTCAATGTCGC‑GG‑A‑CC‑GGGGAAAACTTCCTGCAAAAACGCCCGGCAACGGCCTGCTGACTAAGACCTAGTTGTTtgcggg                                                            >  1:355621/1‑148 (MQ=255)                                                                                         cAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTTTTTTCTTTTATGTCGC‑GG‑A‑CC‑GGGGGAAACCTCAGGGAAAATCGCCCCGCCAAGGCCCGCTGACTAAAACCTTGTTTTTtgcgcg                                                            >  2:346523/1‑147 (MQ=255)                                                                                           gACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCAtt                                                        <  2:503380/151‑1 (MQ=255)                                                                                           gACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGGCTTATTTTTTTTCATGTTGC‑GT‑A‑AC‑GCGGGAAACCTCAGGCAAAAGCGCTCCGCAACGGGCCGCCGACCAAGTCCCTGGTGTTTGCGGGCAtt                                                        >  1:484211/1‑151 (MQ=255)                                                                                             cgcgCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTAc                                                      <  1:401564/151‑1 (MQ=255)                                                                                                gcATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGGGCCTTATTTTCTTCTATGTTGC‑GT‑A‑CC‑GTGGGAACCCTCATGCAAAAACGCCCGGCAACGGCCTGCTGACTAAGTCCTAGTTTTTTGCGAGCAtt                                                        >  1:732083/1‑146 (MQ=255)                                                                                                             tatatCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGGCTTTTCTTTTTCAATGTCGC‑GG‑A‑CC‑GGGGGAAACCTCAGGGAAAAACGCTCGGCAACGGGCCGCCGACTAAATCCCAATTTTTTGCGAGCATTAAGCAGGGGAGCGCCggg                                      >  1:24979/1‑150 (MQ=255)                                                                                                                atCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagt                                   <  2:695887/151‑1 (MQ=255)                                                                                                                 tCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCCTTTTTTCTTTAATGTTGC‑GG‑A‑CC‑GTGGGAAACTTCCGGGCAAAAGGCGCGGGCACGGCCCGCCGTCTCAGTCCTCGTTGTTTGCGGGCATTACTCAGGGGGGCGCCGggag                                    >  2:336571/1‑149 (MQ=255)                                                                                                                  ccGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTaa                                 <  1:102848/151‑1 (MQ=255)                                                                                                                     cgcAGGAATGTTGAAGCAAGATCGGC‑GGGTGGCTTATCTTCTTTTATGTCGC‑GT‑A‑CC‑GTGGGAAACTTCAGGCAAAAGCGCTCGGCCACCGCCCGCTGGCTAAAGCCTAGTTGTTTGCGGGGAATAAGCAGGCGAGCGCCGGGAGGaa                                 >  2:1072088/1‑148 (MQ=255)                                                                                                                      gcAGGAATGTTGAAGCAAGATCGGCGGGGGGCCTTTTTTTTTTC‑ATGTAGC‑GG‑ACCC‑GGGGAAACC‑TCAGGCCAAAACGCGCGGCAACAGGCCGCTGACTGAGGGCTCGGTGTTTGCGGGCCTTAAGCAAGGGGGcggcggggggaaatgg                             >  2:1094685/1‑141 (MQ=255)                                                                                                                      gcAGGAATGTTGAAGCAAGATCGGC‑GGGTGGCCTTTCTTTTTCCATGTTGC‑GG‑G‑AC‑GCGGGAAACTTTAGGCAAAAACGCCCCGCCACCGCCCGCCGGCTAAATCCTAGTTGTTTGCGGGGATTACGCCGGCGGGGGCCCGGAGGAAgggg                             >  2:935276/1‑148 (MQ=255)                                                                                                                        aGGAATGTTGAAGCAAGATCGGC‑GGGGGCCTTTTCTTTTTTCATTTCGC‑GT‑A‑CC‑GTGGAAAACTTCTGGCAAAATCGCTCGGCAACGGCCCGCTGACTAAATCCTAGTTTTTTGCGCGCATTTCGCAGGCGCGGGCCGGGAGTAAAGGCtt                           >  1:243785/1‑151 (MQ=255)                                                                                                                         ggAATGTTGAAGCAAGATCGGC‑GGGTGCCTTTTCTTCTTCAATGTTGC‑GT‑A‑CC‑GCGGGAAACTTCACGCCAAAGCGCTCGGCCACGGCCTGCTGTCTAAGTCCTTGTTTTTCGCGGGCATTACGCCGGCGAGCGCCGGGAGGAAGTGGtt                           >  1:258809/1‑150 (MQ=255)                                                                                                                         ggAATGTTGAAGCAAGATCGGC‑GGGGGCCTTCTCTTTTTCTATGTTGG‑GTAA‑CG‑GGGGAAACTTCCAGGCAAA‑TTACCCGGCAACAGGCCGCTGGCTAAATGCTAATTTTTCTCCCGCATTTCTCAGGGGAGCGCCgggcgggagtggt                            >  2:297665/1‑138 (MQ=255)                                                                                                                                 gAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGt                  <  2:791688/151‑1 (MQ=255)                                                                                                                                 gAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGt                  <  2:866520/151‑1 (MQ=255)                                                                                                                                                  gTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGt  <  2:909939/151‑1 (MQ=255)                                                                                                                                                   tGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGt  >  1:312882/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                          AAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGC‑GGGTGCCTTATCTTCTTCAATGTCGC‑GT‑A‑CC‑GTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGT  >  NC_000913/2642544‑2642838

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