Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,080,837	0	G	T	40.4%	53.7 / 29.2	57	intergenic (+24/‑519)	putP/efeU	proline:sodium symporter/ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base G (11/23);  new base T (23/0);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 5.59e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCA  >  NC_000913/1080690‑1080983                                                                                                                                                    |                                                                                                                                                      cAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGtttttttt                                                                                                                                                    >  1:853469/1‑151 (MQ=255)       tGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCtt                                                                                                                                               >  2:1103246/1‑150 (MQ=255)        gggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTc                                                                                                                                              <  2:191128/150‑1 (MQ=255)          gATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTcagc                                                                                                                                           <  2:1196411/151‑1 (MQ=255)          gATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTcagc                                                                                                                                           <  2:998557/151‑1 (MQ=255)          gATTGTAGTGTTTAGTTTGCTGGGTAAAGCACCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTcagc                                                                                                                                           <  2:793394/151‑1 (MQ=255)            ttGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGtttttttttggtt                                                                                                                                               >  1:777033/1‑140 (MQ=255)               tAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGGTTcagccgggg                                                                                                                                      >  1:209305/1‑146 (MQ=255)                aGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑ggg                                                                                                                                   <  2:1385083/151‑1 (MQ=255)                    tttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCCCTATCATTCGGGTCCGCCGTCACGGTTGCAGGAAAGCTAAGTGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGtt                                                                                                                                 <  1:982381/149‑1 (MQ=255)                    tttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTTAACA‑AC‑ggg                                                                                                                                   >  1:623379/1‑147 (MQ=255)                     ttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCA‑GC‑ggg                                                                                                                                   >  2:1235907/1‑146 (MQ=255)                        gTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGATTTAGCA‑GCAGGGTTggccc                                                                                                                            >  2:367959/1‑147 (MQ=255)                         tttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCA‑GC‑GGGGTGCG‑CTcc                                                                                                                         >  2:902783/1‑151 (MQ=255)                           tGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCCGCA‑GC‑GGGTTGGG‑CTcccc                                                                                                                       >  2:1104824/1‑150 (MQ=255)                            gCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTAGCA‑GC‑GGGGTGCG‑Cccccct                                                                                                                      >  2:353858/1‑149 (MQ=255)                                  tAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTGTTTGTCCTCTCCA‑GC‑GGGGTGCG‑CCCCCtt                                                                                                                      >  1:43988/1‑145 (MQ=255)                                    aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCA‑GC‑GGGGTGCG‑CTCCCCTAAtgtg                                                                                                                >  2:74720/1‑149 (MQ=255)                                    aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CCCCCTTAATGTGcc                                                                                                              >  2:695147/1‑151 (MQ=255)                                    aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGcc                                                                                                              <  2:499386/151‑1 (MQ=255)                                       cgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGGA‑GC‑ggggttcggccccctta                                                                                                                     >  1:644892/1‑128 (MQ=255)                                        gcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTc                                                                                                            <  1:1327321/149‑1 (MQ=255)                                         ccGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTTTTTTGCTTCTGCC‑GC‑GGGGTGGGCCCCCCTTAAT‑TGGCTCGcc                                                                                                         >  1:989537/1‑151 (MQ=255)                                          cGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTTCTTCCGCA‑GC‑GGGGTGGG‑CTCCCTTTAAATGCCCCGccc                                                                                                        >  1:404643/1‑150 (MQ=255)                                           gTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGGTGCG‑CCCCCCTTATTTGCCTcccc                                                                                                         >  2:978747/1‑149 (MQ=255)                                            tCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTTTTTTGCTTTCGCA‑GC‑GGGGTTCG‑CCCCCCTTATGTGcc                                                                                                              >  1:682476/1‑143 (MQ=255)                                             cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCata                                                                                                      <  1:33134/150‑1 (MQ=255)                                               gcggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATaa                                                                                                   <  2:779258/151‑1 (MQ=255)                                                  gcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAAtt                                                                                                <  2:475270/151‑1 (MQ=255)                                                    gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGa                                                                                              <  1:1402504/151‑1 (MQ=255)                                                     aTGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCCGCA‑AC‑GGGGTTCG‑CCCCCTTAATGGGCCTCGCCCAATaaaat                                                                                                >  2:517254/1‑146 (MQ=255)                                                     aTGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGGTCTACA‑GA‑GGGGTGGG‑CGCCCCTAAAGTGGCCTGCCCTATaaa                                                                                                  >  1:1366367/1‑146 (MQ=255)                                                      tGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCA‑GC‑GGGGTGCG‑CCCCCTTAATGTGCCTCGCCCTATAAATTTaa                                                                                             >  2:12559/1‑150 (MQ=255)                                                        cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATgg                                                                                          <  2:671670/151‑1 (MQ=255)                                                        cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATgg                                                                                          <  1:1383551/151‑1 (MQ=255)                                                        cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTGGTGTCAGCA‑GG‑GGGGAGCG‑CTCCCCTTAAATGTCCCTCCACATAAATTTaa                                                                                             >  1:756837/1‑148 (MQ=255)                                                                          gCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGCTCAGCA‑GC‑GGGGTTCG‑CCCCCTTTATTTTGCCCGCCCTATTAAATGGATGGGGCAGGGGGCGCGCCggg                                                                        >  1:850239/1‑151 (MQ=255)                                                                             gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCA‑GC‑GGGGTGCG‑CTCCCTTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAggggg                                                                      >  2:147712/1‑150 (MQ=255)                                                                             gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCA‑GC‑GGGGTGCG‑CCCCCTTTATGTGGCTCGCCATATAAATTGAAAGGTGCAGGGGGGGCGCAggggg                                                                      >  1:762179/1‑150 (MQ=255)                                                                             gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                     <  2:1302452/151‑1 (MQ=255)                                                                              aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCA‑GC‑GGGGTGCG‑CTCCCTTTATGTGCCCCCCCCTATAAATTTAATGGGGCAGGGGGCGCGCAGGGggcg                                                                    >  2:932369/1‑151 (MQ=255)                                                                              aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                     <  1:451461/150‑1 (MQ=255)                                                                                  gcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCATCA‑GC‑GGGGTGCG‑GCCCCTTAATGTGTCCCGCCATATTAAATGAAATGTGCAGGGAGCGCGCAGCgggggg                                                                   >  2:502132/1‑148 (MQ=255)                                                                                    aCTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGGTTTACCC‑GC‑GGGGTGCG‑CCCCCTTAATTAGGCTCCCCCCATTAATTTAATGGTGGAggg                                                                                   >  1:1325471/1‑130 (MQ=255)                                                                                         cATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATcgcc                                                          <  2:287713/150‑1 (MQ=255)                                                                                             cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAACA‑GC‑GGGGTGCG‑CTCCCTTAATGTGCCTCGCCCTATAAAATGGATGGTGCAGGGAGCGCGCAGGGGGCGGCCAAAcgccccc                                                       >  2:533477/1‑149 (MQ=255)                                                                                             cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCA‑CG‑GGGGTGGG‑CGCCCCTGATGGGGCTCCCCCTATTAATTGAATGATGCCGGGAGCGCGGAGGGGGCGGCCAACcgcggccg                                                      >  1:278491/1‑150 (MQ=255)                                                                                             cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGAGAGGGCGCAGGGGGCGGCCAATcgccgccg                                                      <  2:180379/150‑1 (MQ=255)                                                                                                   cgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CCCCCTTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGGGGCCCATCGCCGCCGccccc                                                 >  2:579799/1‑149 (MQ=255)                                                                                                     ccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCctgct                                             <  1:991564/151‑1 (MQ=255)                                                                                                      cgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTTAGCA‑GC‑GGGGTGGG‑CTCCCTTTAACTGCCCCGCCCCATAAAATGAAAGGTGGCGGGGGCGGGCAGGGGGCGGCCCAACGCCGCCGccccc                                                 >  2:604061/1‑146 (MQ=255)                                                                                                          aCGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTccc                                        <  1:365369/151‑1 (MQ=255)                                                                                                                 cAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGGTGCG‑CTCCCCTTATGTGCCTCGCCCTATAAATTGAATGGTGCAAGGGGGGCGCCGGGGGGGGGCAATCGCCGCCGCCCCCTGCCGGCCCCGcc                                    >  2:1076577/1‑148 (MQ=255)                                                                                                                    gAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCggg                              <  2:717629/151‑1 (MQ=255)                                                                                                                        gCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCACCC‑GGGTTCCG‑‑TCCCTTAATGTGCCTCGCCCTATAAAATGAATTGTGCAGGGGGAGCGCAGGGGGGGGCCAATCGCCGCCGCCCCCCGCTGTCCCCGCCTTCgggg                             >  2:618952/1‑148 (MQ=255)                                                                                                                            aGGGACTTAGCCTGCGGCGGTTTTGTTTTGGTTTAGCA‑GC‑GGGGTGCG‑CTCCCCTAATTTGGCTCGCCATATAAATTGAATGGGGCCGGGGGGGCGCCAGGGGGGGCCAAACGCCCCCGCCCCCCGCCGTTCCGGCCCTCgggg                             >  1:1221169/1‑144 (MQ=255)                                                                                                                               gACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAg                   <  1:596537/151‑1 (MQ=255)                                                                                                                                       cTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGAtttt            <  2:228118/150‑1 (MQ=255)                                                                                                                                              ggTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccac    >  1:24758/1‑151 (MQ=255)                                                                                                                                                 tttgtttgGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacca  >  1:797258/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCA  >  NC_000913/1080690‑1080983

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