Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,375,123	0	T	G	43.3%	100.6 / 13.1	60	G85G (GGT→GGG)	dcuD	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (7/27);  new base G (26/0);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 8.34e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCTC  >  NC_000913/3374978‑3375272                                                                                                                                                  |                                                                                                                                                        tgtgCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATg                                                                                                                                                     >  2:313237/1‑150 (MQ=255)         gtCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACa                                                                                                                                           <  1:124943/151‑1 (MQ=255)                      tCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAg                                                                                                                               <  2:164391/150‑1 (MQ=255)                       cGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTG‑CCCGCCACAATGAGCGCATAgg                                                                                                                              >  2:911054/1‑150 (MQ=255)                        ggtggtGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGCGGGGGGGTTTG‑CCCCCCACATTGAGCGCATAgggg                                                                                                                            >  2:349683/1‑151 (MQ=255)                           ggtgTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCa                                                                                                                         <  2:392249/151‑1 (MQ=255)                           ggtgTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCa                                                                                                                         <  2:567115/151‑1 (MQ=255)                           ggtgTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCa                                                                                                                         <  2:748515/151‑1 (MQ=255)                            gtgtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATT‑CCCGCCACATGGGGCGCAAAgggg                                                                                                                            >  1:2637/1‑147 (MQ=255)                            gtgtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAg                                                                                                                        <  1:129504/151‑1 (MQ=255)                                   cGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTG‑CCCCCCACATGGGGGGCCTAAGGGGCAggcccg                                                                                                                   >  2:543942/1‑144 (MQ=255)                                    gATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGa                                                                                                                <  1:91509/151‑1 (MQ=255)                                      tCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTG‑CCCCCTACATGGGGGGCCTAGGGGGCAGTCGCGCGATg                                                                                                              >  1:374643/1‑151 (MQ=255)                                       cccACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATgg                                                                                                             <  2:593931/151‑1 (MQ=255)                                       cccACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATg                                                                                                              <  1:530700/150‑1 (MQ=255)                                        ccaccaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGt                                                                                                            <  2:200293/151‑1 (MQ=255)                                         caccaAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTg                                                                                                           <  1:573718/151‑1 (MQ=255)                                            caAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTT‑CCCGCCACATGGGGGGCATAGGGGGCAGTcgcgcgcgggggagg                                                                                                        >  1:737440/1‑142 (MQ=255)                                              aaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTG‑GCCGCCACAAGGAGGGCATAGGGGGCAGTCGCGCGCTGGGGAGcgt                                                                                                      >  1:895390/1‑149 (MQ=255)                                                   agcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTaa                                                                                                 <  2:362683/151‑1 (MQ=255)                                                      ggCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTTGGCGGGTATG‑CCCGCTACATGGGGCGCATAGGGGGCAGGGGCGCGGTGGGGAGCTTTTTTAgcc                                                                                              >  1:704991/1‑151 (MQ=255)                                                       gCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAgccg                                                                                             <  1:537389/151‑1 (MQ=255)                                                              gATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTa                                                                                     <  1:641854/151‑1 (MQ=255)                                                              gATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTa                                                                                     <  2:951367/151‑1 (MQ=255)                                                                  gtcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGTTCTT‑CCCGGGCCACGGAGCGAAAAAGGGCGCGGCCCGCGCTGGGGAAGTGGTTATTCCG‑CCCGGTaaaa                                                                                  >  2:961778/1‑150 (MQ=255)                                                                   tcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑CCCGCTACATGGAGCGGATAGGGGGCCGTCGCGCGGTGGTGAGGTTGTTTAGCCG‑CCCCTTtaaaa                                                                                 >  2:434558/1‑149 (MQ=255)                                                                      atatTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCAt                                                                             <  2:489053/151‑1 (MQ=255)                                                                          ttATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTT‑CCCCCCCCATGGAGCGCACAGGGGGCAGGCGCGCGGTGGGGGGGTTGTTTAGCCG‑cccccttaaa                                                                                  >  2:162411/1‑135 (MQ=255)                                                                            ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATtcgc                                                                         <  2:684886/149‑1 (MQ=255)                                                                              aaTGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑GCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGGTGGGGAGCTTGTTAAGCCG‑CCCCTTAAAACTCATtcgctcgc                                                                     >  1:101788/1‑151 (MQ=255)                                                                                  aaaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATT‑CCCGCTACATGGGGCGCCTAGGGGGCAGGCGCGCGATGGTGAGCTTGTTTAACCG‑CCCCTTAAAACTCCTTCCCTCGcc                                                                    >  1:666398/1‑148 (MQ=255)                                                                                  aaaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGta                                                                 <  1:895512/151‑1 (MQ=255)                                                                                   aaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGAGGCGGGCATGACCCCCGCCGTGG‑GCGGACATGGGCAAGTCGCGCCAGGGGGTGCGTGTGTAGCCG‑CCCTTTAAAAAACAATCACTCCCCGtct                                                                >  2:622618/1‑149 (MQ=255)                                                                                    aTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATG‑TCCGCCACATTGGGGGCCTAGGGGGCAGGCGCGCGGTGGGGAGGTTGTTTAACCG‑CCCGGTAAAAATCATTCTCTCGcc                                                                    >  2:319554/1‑146 (MQ=255)                                                                                    aTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑CCCGCTACCTGGAGCGCATAGGGGGCCGTCGCGCGCTGGGGAGCTTGTTTAACCG‑CCCCGTAAAACACATTCGCTcccccgat                                                                >  1:731725/1‑146 (MQ=255)                                                                                      ccTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTAtatt                                                             <  1:979351/151‑1 (MQ=255)                                                                                         gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGCGGTTATG‑CCCGCTAAATGGGGCGCCTAGGGGGCAGTCGCGCGGTTGTGAGGGTGTTTAGCCG‑CCCGGTAAAACTCCTTCGCTCGCCGTATattttt                                                          >  1:738/1‑151 (MQ=255)                                                                                             aTGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTG‑CCCGCTACATGGAGCGCCTAGGGGGCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTTTTCTGt                                                      >  1:48314/1‑151 (MQ=255)                                                                                             aTGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑GCCGCTACATGGAGCGCATAGGGGCCCGTCGCGCGATGGGGGGCTTGTTAAACCG‑CCCGTTAAAAATCATTTGCTCCCCGTATattatt                                                          >  1:292288/1‑147 (MQ=255)                                                                                               gCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑TCCCCTCCATGGAGGGCATAAGGGGCCGTCGCGCGGTGGGGGGCTTTTTTAGCCG‑CCCGTTAAAAATCATTCGCTCGCCGGATattttttttt                                                      >  1:108186/1‑145 (MQ=255)                                                                                                cTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑TCCGCCACATTGGGCGGCTAGGGGGCCGGCGCGCGGTGGGGAGGTTGTTTAGGCG‑CCCCCTAAAAATCATTCCCCCGCCGGATattatttttttgg                                                   >  2:290415/1‑144 (MQ=255)                                                                                                 tGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑CCCCCTACCTGGAGCGCATAAGGGCCAGGCGCGCGGGGGGGGGCGTGTTAAGACC‑CCCGGTTAAAATCATTCTCTCGCCCGATATTTTTTTTTCggg                                                  >  2:834262/1‑150 (MQ=255)                                                                                                  gTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTT‑GCCGCTACATTGGGCGCCCAGGGGGCAGGCGCGCGGTGGGGGGGTTGTTTAGCCGCCCCGTTAAAAC‑CATTCGCTCGCCCGATattttttttttgg                                                   >  2:375026/1‑142 (MQ=255)                                                                                                     cAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATG‑CCCCCCACATGGGGCGGATAGGGGGCAGTCGCGCGGTGGGGAGCTTGTTTAGCCG‑CCCGTTAAAACTCATTTGCTCGCCGGATATTATTTTTTCGGCCACt                                              >  1:752038/1‑151 (MQ=255)                                                                                                           cATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTTTG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTc                                        <  1:985064/151‑1 (MQ=255)                                                                                                           cATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTATG‑GCCCCCGCATGGAGAGGAAAGGGGGCGGTCGCGGGGGGGTGGACTTGTTTAACCG‑CCCGTTAAAAATCATTCCCTCCCCGGATTTTTTTCTGTCTGCGACtt                                             >  1:580366/1‑146 (MQ=255)                                                                                                            aTTGCCGGATTGGGGCTGTCGATTATGGCGGTTGGGGGGTATG‑CCCGCTAAATGGAGGGGATATGGGGCAGTGGGGCGGTGGGGAGGTTGTTTAGGCG‑CCCGTTAAAACCCATTTTCTCTCCCGATATTATTCTGTGGGGAAcct                                             >  1:439698/1‑143 (MQ=255)                                                                                                            aTTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCTGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGt                                         <  2:954746/149‑1 (MQ=255)                                                                                                               gCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCCACATGGAGCGGATAGGGGGCCGTCGCGCGGTGGTGAGCTTGTTAAGCCG‑CCCCTTAAAAATCAATTGCCCGCCGGATATTTTTTTGTCGGGAAATTACGTCATCg                                    >  1:714415/1‑151 (MQ=255)                                                                                                                  ggATTGGGGCTGTCGATTATGGCGGGGGGGGGGTATG‑CCCCCCACATGGAGCGCATAGGGGGCCGGCGCGCGGTGGGGGGGTTGTTTAACCG‑CCCCCTAAAAATCATTTGCCCCCCGGATATTTTTTTTTTGGCCACTTTAGTCATCGGcc                                 >  2:796848/1‑151 (MQ=255)                                                                                                                    aTTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACTTCATCGGCCa                                <  2:978766/150‑1 (MQ=255)                                                                                                                        gggCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCa                           <  2:939483/151‑1 (MQ=255)                                                                                                                        gggCTGTCGATTATGGCGGTGGGCGGGTATG‑CCCGCCACATGGAGCGCATAGGGGGCAGGCGCGCGGTGGGGGGCTTGTTTAGCCG‑CCCCTTAAAAATCATTTCCTCGCCGCATATTATTTTTTCGGCAACCTACGTCATCGGCCaaaac                            >  1:295998/1‑148 (MQ=255)                                                                                                                          gCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATg                         <  2:121455/151‑1 (MQ=255)                                                                                                                             gTCGATTATGGCGGTGGGGGGGTATG‑TCCCCCACATGGAGCGCATAGGGGGGAGACGCGCGATGGTGGGGGTGTTAAGCCG‑CCCGCTTAAACTCATTTGCCCGCCGGATATTATTTTTTTGGCAACTTACGTCATCGGCCAAATCa                           >  2:767236/1‑146 (MQ=255)                                                                                                                             gTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGgcg                      <  2:645606/151‑1 (MQ=255)                                                                                                                                       gCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTttatt             <  1:620597/150‑1 (MQ=255)                                                                                                                                         ggTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTAcc          >  2:307686/1‑151 (MQ=255)                                                                                                                                             ggCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAgcg      <  1:13478/151‑1 (MQ=255)                                                                                                                                                 gTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCTc  >  2:900679/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                        TGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATG‑CCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCG‑CCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCTC  >  NC_000913/3374978‑3375272

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: ‑