Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	43.3%	103.5 / 15.0	61	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/30);  new base G (26/0);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 7.84e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925754‑926050                                                                                                                                                    |                                                                                                                                                        ttCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGtt                                                                                                                                                       >  1:809233/1‑149 (MQ=255)  tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCg                                                                                                                                                     <  1:1033653/150‑1 (MQ=255)        ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAccccc                                                                                                                                              >  1:334662/1‑151 (MQ=255)        ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAAccccc                                                                                                                                              <  1:598805/151‑1 (MQ=255)            aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                           <  1:752960/150‑1 (MQ=255)             cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGAccccc                                                                                                                                              >  1:885738/1‑146 (MQ=255)              aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGCAccccc                                                                                                                                              >  2:743569/1‑145 (MQ=255)                              cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCATACGTGTCCGTATacaaa                                                                                                                          >  2:464841/1‑147 (MQ=25)                                 aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGTATACACAAttt                                                                                                                      >  1:849736/1‑150 (MQ=255)                                 aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTc                                                                                                                     <  1:674559/151‑1 (MQ=255)                                  ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTcc                                                                                                                    <  2:230177/151‑1 (MQ=255)                                  ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTcc                                                                                                                    <  1:884503/151‑1 (MQ=255)                                  gatAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTcc                                                                                                                    <  2:319906/149‑1 (MQ=255)                                   gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCa                                                                                                                   <  1:357339/151‑1 (MQ=255)                                    tAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGCTGCAGTATACAAACctttc                                                                                                                    >  2:1083630/1‑147 (MQ=25)                                     aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTCGGGACCCCCATACGTTGCCGGATACACACttt                                                                                                                      >  1:118425/1‑146 (MQ=25)                                       cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCCTACGTTTCCGGATACAAACCTTTCAGGCg                                                                                                               >  1:955958/1‑151 (MQ=21)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCt                                                                                                           <  1:234578/149‑1 (MQ=255)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCt                                                                                                           <  1:892602/149‑1 (MQ=255)                                               aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTA‑GGGGGGGGATTGGACCCCCCCTACGTTGCCGGATACACACTTTCCAGGCggggccc                                                                                                         >  1:123481/1‑143 (MQ=21)                                                    cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGcc                                                                                                 <  2:408714/151‑1 (MQ=255)                                                           tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCCATACGTTTGCCTATAAACACTTTTCCGGCGGGCCCC‑TTCCTCCAccc                                                                                             >  2:194928/1‑146 (MQ=18)                                                           tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTCGAACCCCCCCTACGTTGCCCTATACCCACTTTTCCCGCGTGCTCC‑CTccacccccc                                                                                             >  2:1092636/1‑140 (MQ=18)                                                              tCACTGCCGTACAGACAGATAAAATGCGGAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGaca                                                                                        <  2:662214/150‑1 (MQ=255)                                                                   gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGACCCCCCCTACGCTTCCGGATACACACCTTTCAGGCGTGCTCC‑CTCCGCCCCCCGGACACCTCAc                                                                                  >  2:408951/1‑151 (MQ=17)                                                                   gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTc                                                                                    <  2:620893/149‑1 (MQ=255)                                                                          agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCAAACGTTGCCGCATAAACACTTTCCAGGCGGGCTCC‑TTCAGCCACCCCGACACCCCCCCaaaat                                                                            >  1:530538/1‑148 (MQ=14)                                                                            acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGtt                                                                         <  2:334663/151‑1 (MQ=255)                                                                                      tGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCGAAAAGTTGCCGTTTACAAAATTTTCCGGGGTGCTCC‑TTCAACCACTCCGGCACCTCACCAAAATTTTTTTTTAcc                                                                 >  1:166095/1‑149 (MQ=11)                                                                                         gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAAcc                                                            <  1:528125/151‑1 (MQ=255)                                                                                           aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGGACCCCCCCTACGTTTCCGTATAAACACTTTCCAGGGGGGGTCC‑TTCAGCCACCCGGACACCCCCCCCAAATTTTTTTTTACCAAAcccc                                                          >  2:355898/1‑149 (MQ=2)                                                                                            aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCa                                                         <  2:523238/151‑1 (MQ=255)                                                                                              aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCGAGAAGTTGTCGTTTACACACTTTTCAGGCGGGCTCC‑TTCAGCCACCCGGACACCCCACCCAAATGTTTTTTTACCAAACCTCATg                                                       >  1:99558/1‑151 (MQ=14)                                                                                                 aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGA‑GGGGGGGGTTTGTACCCCCCGCACGGTTCCGTATACACACTTTCCAGGGGGGGTCC‑TCCCTCCCCCCCGGGACCTCCCCCAATTTTTTTTCTTCCTAAcccccttgg                                                     >  2:456387/1‑143 (MQ=2)                                                                                                    tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGAAACCCCCATAACCTGTCGTATATACACCTTTCAGGGGGGGTCCTTTCCCCCACCCGG‑CACCCCCCCCAAATTTTTTTCTTCCAAAACTCATgggggg                                                  >  2:57292/1‑150 (MQ=2)                                                                                                     cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATACACACTTTTCAGGCGTGCTCC‑TTCCGCCACCCCGGCACCCCCCCAAAATGTTTTGCTACCAAACCTCATTGgggg                                                  >  1:953706/1‑149 (MQ=14)                                                                                                      gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTCGAA‑CCCCGCCACGTTGGCGTAAAAAAACTTTTCAGGGGGGCTCC‑TTCCGCCCCCCGGGGACCCCCCCCAAATGTTTTTTTACCCAACCCCCTgggggg                                                  >  2:592751/1‑148 (MQ=2)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGACCCCCCATACGTTGCCGTATAACCACCTTCCCGGCGGGCTCC‑CTCAGCCACCCCGACACCCCACCCAAATGTTTTTCTACCAAACCTCATGGGGGGCAACgggg                                          >  2:296249/1‑149 (MQ=17)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACGGTGCCGTATAAACACTTTCCCAGCGGGGTCC‑TTCAGCCACTCCGGCACCCCACCAAAATGTTTTGCTACCCCACCCCCTTGGGGGCAAAAGGgcg                                        >  1:239964/1‑151 (MQ=2)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg                                        <  2:809259/151‑1 (MQ=255)                                                                                                                aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGACCCCCCGCACGTTGTCGCTTATAAACTTTTCATGCGTGGTCC‑CTCAGCCACCCGGACACCCCACCCAAATCTTTTTTTTCCAAACCACATGGGCGGGAACgagg                                          >  1:122268/1‑146 (MQ=2)                                                                                                                 cGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctac                                    <  2:938738/151‑1 (MQ=255)                                                                                                                  gAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCATACGTTGCCGTATACACACTTTTCAGGCGGGCTCC‑TTCAGCCACCCGGACACCCCACCCAAATTTTTTTCTACCAAAACTCATTGGTGGGAACGGGGGGctac                                    >  1:240747/1‑150 (MQ=14)                                                                                                                    gCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACtat                                 <  1:695352/151‑1 (MQ=255)                                                                                                                         tcttTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCGCTACGTTGCCGTATATACACCTTTCCGGCGTGCCCC‑TTCAGCCACTCCGGCACCCCACCAAATTGTTTTGCTACCAAAACTCATTGGTGGCAACGGGGGGCCACTATAggg                             >  1:661754/1‑150 (MQ=17)                                                                                                                         tcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                               <  1:345649/148‑1 (MQ=255)                                                                                                                          cttTAAATATGGCGGTGA‑GGGGGGGGGTTGGACCCCCCATAAGTTTCCGTATAAACACTTTTCCGGCGGGGCCC‑TTCTGCCACCCCGGCCCCCCCCCAAAATTTTTTTCTTCCAAACCTCCTGGGTGGGCAAGGGGGGCGACTATAgg                              >  2:45260/1‑148 (MQ=9)                                                                                                                              aaaTATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTg                        <  1:446679/150‑1 (MQ=255)                                                                                                                                atatGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                     <  2:995291/151‑1 (MQ=255)                                                                                                                                atatGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                     <  1:283675/151‑1 (MQ=255)                                                                                                                                  atGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACCCCCTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                    <  1:743551/150‑1 (MQ=255)                                                                                                                                  atGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                   <  1:408956/151‑1 (MQ=255)                                                                                                                                  atGGCGGTGA‑GGGGGGGATTCGAACCCCCCATAAGTTGCCGTATACACACTTTTCAGGCGGGCTCC‑TTCAGCCCCTCGGACACCCCACCCAATTGTTTTGCTACCCAACCTCATGGGTGGCAACGGGGCGCGACTATAGGGGGGTgg                       >  2:1071722/1‑147 (MQ=18)                                                                                                                                     gCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa                 >  2:854236/1‑150 (MQ=255)                                                                                                                                       ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg              >  2:518655/1‑151 (MQ=255)                                                                                                                                        gTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg             <  1:863054/151‑1 (MQ=255)                                                                                                                                          gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           <  1:686328/151‑1 (MQ=255)                                                                                                                                          gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           >  2:938646/1‑151 (MQ=255)                                                                                                                                               gggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga       <  1:717819/151‑1 (MQ=255)                                                                                                                                                ggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  2:1013444/151‑1 (MQ=255)                                                                                                                                                 gggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag     <  2:279228/151‑1 (MQ=255)                                                                                                                                                    aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:118434/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCC‑TTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925754‑926050

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