Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,465,167	0	G	T	42.9%	69.0 / 28.1	56	C29F (TGT→TTT)	gspM	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (9/23);  new base T (24/0);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 9.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAA  >  NC_000913/3465018‑3465314                                                                                                                                                       |                                                                                                                                                    ttCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CCCGGCGGTTTttgt                                                                                                                                                    >  2:613465/1‑151 (MQ=255)  tCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGGGGtttttttt                                                                                                                                                   >  1:277260/1‑151 (MQ=255)       tGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCC‑CTCGGCGGTTTTTGTTgggt                                                                                                                                               >  1:630135/1‑148 (MQ=255)         gAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCC‑CCCGGCGGTTTTTTtttgtttg                                                                                                                                             >  1:494894/1‑150 (MQ=255)          aaaaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTTTTTGGTTggg                                                                                                                                           >  2:991256/1‑151 (MQ=255)          aaaaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGGTTTTTTTTTGCTTggg                                                                                                                                           >  2:349284/1‑151 (MQ=255)            aaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGGTTggggg                                                                                                                                         >  1:848633/1‑150 (MQ=255)            aaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGt                                                                                                                                         <  2:997128/151‑1 (MQ=255)             aaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTg                                                                                                                                        <  2:248410/151‑1 (MQ=255)             aaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGGTCTCTCG‑CCCGGCGGTTTTTTTTGGGTTgggggg                                                                                                                                        >  1:649051/1‑149 (MQ=255)              aaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCG‑CTCGGCGTTTTTTTTTGGCTTgggggg                                                                                                                                        >  1:656090/1‑148 (MQ=255)              aaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCC‑CTCCGCGGTTTTTTTTTGGTTgggg                                                                                                                                          >  2:1039465/1‑148 (MQ=255)               aCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTtttgggtttggg                                                                                                                                          >  2:553527/1‑147 (MQ=255)                   aaaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCCGCGGTTTTTTTTGGCTTGGGGGGGTaaa                                                                                                                                   >  1:393556/1‑150 (MQ=255)                      aaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCg                                                                                                                               <  1:402210/151‑1 (MQ=255)                       aGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCC‑CCCGGCGGTTTTTTTTTGGTTGGGGGGGTaaaa                                                                                                                                  >  2:773532/1‑147 (MQ=255)                           ccGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgataga                                                                                                                          <  2:1045414/151‑1 (MQ=255)                            cGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgatagat                                                                                                                         <  2:118148/151‑1 (MQ=255)                             ggTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCCCCGGGGTTTTTTTTTGGTTTGGGG‑GGTAAAAACcataga                                                                                                                          >  2:596875/1‑149 (MQ=255)                                cAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATAcc                                                                                                                      <  1:1028635/150‑1 (MQ=255)                                cAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCt                                                                                                                     <  2:634469/151‑1 (MQ=255)                                    ggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACa                                                                                                                  <  2:684217/150‑1 (MQ=255)                                         aTGGACAGTAAGGAACTCAGGGAAATGTTTAAATCATGGCGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGGTTGGGGGGGGAAAACCggaaaataa                                                                                                                       >  2:491566/1‑132 (MQ=255)                                         aTGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CCCGGCGTTTTTTTTTTGCTTGGGGGGATAAAACCGATAGATACCCACAtt                                                                                                                >  1:340954/1‑147 (MQ=255)                                          tGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCCGCGGTTTTTTTTGGGTTGGGGGGGTAAAAACCATAGATAACTACAtt                                                                                                                >  2:369276/1‑146 (MQ=255)                                           ggACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCaga                                                                                                           <  1:896262/150‑1 (MQ=255)                                            gACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGc                                                                                                         <  2:485732/151‑1 (MQ=255)                                               aGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATc                                                                                                      <  1:845049/151‑1 (MQ=255)                                                gTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCG‑CTCGGCGGTTTTTTTTTGCTTTGGGGGATAAAACCCATAGATAACTACATTTCAGAGCATc                                                                                                      >  1:713222/1‑150 (MQ=255)                                                gTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCa                                                                                                     <  2:994063/151‑1 (MQ=255)                                                  aaGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTTTCG‑CCCGGGGTTTTTTTTTGGGCTGGGGTGGTAAAACCGATAGATACCCACATTTCTGAGCa                                                                                                        >  2:895902/1‑146 (MQ=255)                                                   aGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CCCGGCGTTTTTTTGTGGGTTTGGGGGGTAAAACCCATAGAAACCCACATTGCagag                                                                                                          >  2:15884/1‑143 (MQ=255)                                                       aCTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCG‑CCCGGCGGTTTTTTTTGGGTTGGGGGGGTAAAACCGATAAAAACCCACATTGCAGAGCATCAGGccc                                                                                                >  1:818640/1‑149 (MQ=255)                                                       aCTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCAc                                                                                              <  2:649161/151‑1 (MQ=255)                                                        cTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCC‑CCCGGCGGTTTTTTTTGGGTTGGGGGGGTAAAACCGATAGATACCTACATTGCAGAGCATCAGTccccc                                                                                              >  1:993584/1‑148 (MQ=255)                                                          cAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCAcgcg                                                                                           <  1:995984/151‑1 (MQ=255)                                                           aGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCG‑CTCGGCGTTTTTTTTTGGGTTTGGGGGGTAAAAACGATAGATACCTACATTTCAGAGGATCAGTCCCAcgc                                                                                            >  2:413730/1‑149 (MQ=255)                                                             ggAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCC‑CCCGGCGGTTTTTTTTGGCTTGGGGGGGTAAAAACCATAGATACCTTCATTGCAGAGCATCAGTCCCCcgccc                                                                                          >  2:214332/1‑147 (MQ=255)                                                                 aTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGGGGTAAAAACGATAGATACCTACATTGCAGAGCATCAGTCCCCCGCGCaaaaa                                                                                     >  2:772808/1‑150 (MQ=255)                                                                 aTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCaaaaaa                                                                                    <  1:137485/151‑1 (MQ=255)                                                                   gATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAAt                                                                                   <  2:639080/150‑1 (MQ=255)                                                                   gATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CCCGGGGGTTTTTGTTGGGTTGGGGGGGTAAAAACGAAAAATAACTACCTTTCAGAGGATCAGTCCCACCCGCaaaaaa                                                                                    >  1:798140/1‑149 (MQ=255)                                                                     ttAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCG‑CTCGGCGTTTTTTTTTTGGTTGGGGGGGTAAAAACGATAGATACCTACCTTTCAGAGCATCAGTCCCCCGCGCcaaaaaat                                                                                  >  2:185428/1‑147 (MQ=255)                                                                                    cAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATg                                                                 <  2:266343/151‑1 (MQ=255)                                                                                     aGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCC‑CCCGGCGGTTTTTTTTTGGTTGGGGGGGTAAAAACCATAGATACCTACATTTCAGAGCATCAGTCCCCCGCGCAAAAAATTAAAAAGGGCATCAAAAgg                                                                >  1:79532/1‑151 (MQ=255)                                                                                          aaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCa                                                           <  1:134696/151‑1 (MQ=255)                                                                                              cacaTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCG‑CTCGGCGGTTTTTTTTTGCTTGGGGTGGTAAAAACGATAGAGACCCACATTTCCGAGCATCCGTTCCCCCCGCACAAAAATAAAAAAGGCATCAAAATGGTGCAgg                                                         >  1:523458/1‑149 (MQ=255)                                                                                               acaTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CCCGGCGTTTTTTTTTGGCTTGGGGGGATAAAACCGATAGATAACTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAAAAGGACATCAAATGGAt                                                              >  1:492870/1‑143 (MQ=255)                                                                                                caTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGGTTTTTTTTGGCTTGGGGGGATAAAAACGATAGATACCTACATTGCAGAGCATCAGTGCCACGCGCAAAAAAATAAGAAGGACATCAAATGGATGCACGATc                                                      >  1:491028/1‑150 (MQ=255)                                                                                                  tCGGAAAAACAGATCGTTGCGTCAGTGGCGGTTCTCTCG‑CGCGGCGTGTTTTGGTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGTACATCAAATGGATGCAGGATCaa                                                    <  1:616489/150‑1 (MQ=255)                                                                                                                    gCGGCACTGGCGGTTCTCTCG‑CTCGGGGTTTTTTTTTGGCTTGGGGGGGTAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCCAAAAATTAAAAAAGAAATCCAATGGATGCAGGATCAAACGAGCCCCCCTTGttt                                   >  1:641972/1‑149 (MQ=255)                                                                                                                            ggCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTc                           <  1:615963/149‑1 (MQ=255)                                                                                                                                  tctctcG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTgc                    <  2:277328/150‑1 (MQ=255)                                                                                                                                       cg‑ctcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCtt                <  1:204736/149‑1 (MQ=255)                                                                                                                                        g‑ctcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACg             >  2:720708/1‑151 (MQ=255)                                                                                                                                           tcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGc            >  1:365604/1‑150 (MQ=255)                                                                                                                                            cgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCa           <  2:630240/150‑1 (MQ=255)                                                                                                                                                gTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGccc       >  2:141412/1‑150 (MQ=255)                                                                                                                                                    tttGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCa   >  2:205310/1‑150 (MQ=255)                                                                                                                                                     ttgttgGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCaa  >  1:304537/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCG‑CTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAA  >  NC_000913/3465018‑3465314

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