Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,333,757	0	T	G	42.6%	114.6 / 13.1	61	G16G (GGT→GGG)	ispB	octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/34);  new base G (26/0);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 2.22e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCA  >  NC_000913/3333616‑3333906                                                                                                                                               |                                                                                                                                                           tttttCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑Gc                                                                                                                                                  <  1:383157/150‑1 (MQ=255)      cAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCaa                                                                                                                                             <  2:608905/150‑1 (MQ=255)        gTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGCG‑GTTTAT‑ACGGCAAtac                                                                                                                                          >  1:949351/1‑149 (MQ=255)          aCAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCtt                                                                                                                                        <  2:97771/151‑1 (MQ=255)          aCAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTTAT‑GCggga                                                                                                                                              >  1:549335/1‑143 (MQ=255)           cAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTg                                                                                                                                       <  2:887319/151‑1 (MQ=255)            aaTCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTAAT‑gggg                                                                                                                                                >  2:130806/1‑141 (MQ=255)               cACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGagc                                                                                                                                    <  1:826414/150‑1 (MQ=255)                 cccTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTAAT‑TAGGCAATCCTTGagc                                                                                                                                    >  2:354128/1‑148 (MQ=255)                     atatTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑tt                                                                                                                             <  1:26262/149‑1 (MQ=255)                     atatTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑tt                                                                                                                             <  1:391749/149‑1 (MQ=255)                      tatTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTGAT‑TCGGCCAACCTTTagg                                                                                                                                    >  2:610459/1‑142 (MQ=255)                        ttCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTTAT‑ACGGGAAACCTTTAGCAG‑C‑TTAAtt                                                                                                                         >  2:37150/1‑150 (MQ=255)                          ccTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTcc                                                                                                                       <  1:809914/150‑1 (MQ=255)                          ccTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTTAT‑TAGGGAATCCTTTTGCAG‑C‑CTAAtttc                                                                                                                       >  2:1011616/1‑148 (MQ=255)                           cTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTTAT‑ACGGGAAACCTTTAGGAG‑CGTTAATTcgc                                                                                                                      >  2:903939/1‑149 (MQ=255)                             aaCCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGa                                                                                                                     <  1:479773/149‑1 (MQ=255)                             aaCCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACg                                                                                                                   <  2:751292/151‑1 (MQ=255)                                    aaCCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACt                                                                                                           <  2:376939/151‑1 (MQ=255)                                       ccTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGat                                                                                                        <  1:271893/151‑1 (MQ=255)                                       ccTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGat                                                                                                        <  1:272468/151‑1 (MQ=255)                                          aaGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGatcaa                                                                                                     <  2:468452/151‑1 (MQ=255)                                          aaGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTAAT‑GAGGCAAACCTTTAGCAG‑C‑ATTATTTCGACG‑Tccca                                                                                                             >  2:399382/1‑141 (MQ=255)                                              tGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTTAT‑GAGGGAAACCTTTAGCAG‑C‑CTTATTTCGACG‑GCCAACTGatcaaaaa                                                                                                  >  2:543779/1‑147 (MQ=255)                                               gCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGt                                                                                                <  1:346426/151‑1 (MQ=255)                                               gCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTAAT‑GCGGCGATCCTTTTTCAG‑C‑CTAATTCCGACG‑GCCAACTGat                                                                                                        >  1:379900/1‑143 (MQ=255)                                                   ttGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAg                                                                                             <  1:845283/150‑1 (MQ=255)                                                       tCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTAAT‑TCGGGAAACCCTTTGCAA‑C‑TTAATTTCCGCG‑GCCAACTGATTAAAtaggtagggt                                                                                          >  2:429481/1‑147 (MQ=255)                                                          caGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTAc                                                                                     <  1:305668/151‑1 (MQ=255)                                                          caGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTTAT‑GGGGCAATCCTTGAGCAG‑C‑CTTATTTCGACG‑TCCAACTGGTAAATAAATTAGGCTAtt                                                                                       >  2:464634/1‑149 (MQ=255)                                                           aGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACa                                                                                    <  1:712671/151‑1 (MQ=255)                                                             tAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTGAT‑TCGGCAATCCTTTAGCAG‑A‑TTAATATTCACG‑TCCAAATGCTCAAACAGTGATGCTAt                                                                                        >  1:788508/1‑145 (MQ=255)                                                                 gTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCa                                                                              <  1:250404/151‑1 (MQ=255)                                                                   aaTCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgc                                                                            <  1:101141/151‑1 (MQ=255)                                                                      cGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgcggc                                                                         <  2:860334/151‑1 (MQ=255)                                                                       ggggCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGG‑GTTATA‑GGGGAAAACCTT‑TGCAGCC‑TTAATTCCAACG‑CCAACCCGATCAAACAG‑TAGGGTAtt                                                                                       >  1:64158/1‑136 (MQ=255)                                                                       ggggCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgcggcg                                                                        <  1:524045/151‑1 (MQ=255)                                                                             aaaaGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTa                                                                  <  1:300257/151‑1 (MQ=255)                                                                               aaGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTaaa                                                                <  1:610363/151‑1 (MQ=255)                                                                                 gCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGG‑GTAAAG‑GGGGAAACCCTTGGGC‑G‑C‑CTAATTTCCGCG‑TCCCACTGATCAATTAGTTAGGGTATTACATCGTCAGCGGGGGCGGGaaa                                                                >  2:1008651/1‑149 (MQ=255)                                                                                 gCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACg                                                              <  1:9113/151‑1 (MQ=255)                                                                                  cccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGG‑GTTAAT‑‑CGGGGATCCCTTAGCAG‑C‑CTAATTTCCACG‑GCCCACTTATCAAAAAGTTAAGGTCTTTCATCCTCCGcggcgggggggaa                                                                >  2:554764/1‑144 (MQ=255)                                                                                  cccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTTAT‑GAGGGAAACCTTTTGCAG‑C‑CTAATTTCCAAG‑TCCCACCGATCAAACAAATAGGATTTTTCATCCTCAGCGGCGGGGGGGAACgg                                                             >  1:928824/1‑150 (MQ=255)                                                                                   ccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTTAT‑TTGGCAATCCTTTAGCAG‑C‑CTAATTTTGACG‑CCCCAATGATAAAACAATTAGGGTATTTAAATGTCAGCCGCGGGGGGaaaa                                                               >  2:273897/1‑147 (MQ=255)                                                                                       cTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTAAT‑GCGGCAAACCTTTAGCAG‑C‑CTAATTTCCACC‑TCCAACTGATCAATAAATTTGGCTTTTTAATTCTCAGcggggggggggaaagaatt                                                          >  2:165128/1‑139 (MQ=255)                                                                                        ttttGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTTAT‑GCGGGAAACCTTTTGCAG‑C‑CTAATTCCCACGCTCCAA‑TGATCAAAAAATTTTGCTATTTCATCTTCAGCCGCGGGGGGAAAAGGAtt                                                          >  2:226831/1‑148 (MQ=255)                                                                                         tttGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGt                                                       <  1:261845/150‑1 (MQ=255)                                                                                           tGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTcc                                                     <  2:493248/150‑1 (MQ=255)                                                                                            gCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGTGGTTAAT‑‑CGGCAAATCTTTAGCAC‑C‑CTAATTTCCAAC‑TCCCACCGATCAATCACATAGGGTATTTAATCGTTCGCGGCGGCGGGGAACGTATTTGt                                                       >  1:662777/1‑147 (MQ=255)                                                                                            gCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCga                                                   <  1:369921/151‑1 (MQ=255)                                                                                             cGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCgat                                                  <  2:299643/151‑1 (MQ=255)                                                                                              gATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCgatg                                                 <  1:628497/151‑1 (MQ=255)                                                                                                tGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTTAT‑GCAGGAAAACTTTAGGAG‑C‑TTAATTTCGACG‑TCCCACTGATAAAACAGTTATGGTATTATATCGTCAAAGGCGGCGGGGAAAGTATTCGttccctg                                                 >  1:878551/1‑143 (MQ=255)                                                                                                 gAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTAAA‑GAGGGAAACCTTTAGCAG‑C‑CTTATTTCCACG‑TCCAAATTATCAATAAAATATGCTATTATATTTTCAGCGGCGGCGGGAAAAGAATTTTTCCGATGAtt                                              >  2:353963/1‑151 (MQ=255)                                                                                                    tttAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GTTACTATGGGGAAACCTT‑TGCAG‑C‑CTAATTCCGACG‑GCCAACTGATCAATCAGTTAGGGTCTTAACTCCTCAGCGGCGGCGGGAAAcacatttttcc                                                     >  2:966115/1‑131 (MQ=255)                                                                                                     ttAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCt                                           <  2:679863/150‑1 (MQ=255)                                                                                                       aGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTAAA‑TAGGCAAACACTTAGAAA‑C‑CCAATTTCGTCC‑GCCCACCGATCAATCAGTTAGGCTATTACATCCTCCGCGGCCGCGGGGAACGAATTTGTCCCTTGATTGCTGt                                         >  1:79466/1‑150 (MQ=255)                                                                                                              aTCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTg                                 <  1:197543/151‑1 (MQ=255)                                                                                                              aTCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTg                                 <  2:50501/151‑1 (MQ=255)                                                                                                              aTCAATGAGTTAACCGCGCAAGATATGGC‑GGGG‑GGTAAT‑GCGGGAAAACTTTAGGAG‑C‑TTTATTTCCCCG‑TCCAACTCATCAACCAGATAGGGTATTTTATCGTCAGCGGGGGCGGGAAACGTATTTGTCCCATGATTGTTGTTCTggg                                   >  2:792438/1‑148 (MQ=255)                                                                                                                    gAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGa                            <  2:579721/150‑1 (MQ=255)                                                                                                                                              gt‑gtTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCa  <  1:609981/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                           TTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGC‑GGGT‑GTTAAT‑GCGGCAATCCTTGAGCAG‑C‑TTAATTCCGACG‑TCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCA  >  NC_000913/3333616‑3333906

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