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

CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925756‑926050                                                                                                                                                 |                                                                                                                                                          cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGa                                                                                                                                                      <  1:348165/150‑1 (MQ=255)   cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAc                                                                                                                                                    <  1:995670/150‑1 (MQ=255)   cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAc                                                                                                                                                    <  1:316236/150‑1 (MQ=255)     gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAccc                                                                                                                                                  <  1:39457/150‑1 (MQ=255)       cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAacccc                                                                                                                                                >  1:962728/1‑150 (MQ=255)       cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCg                                                                                                                                               <  2:227223/151‑1 (MQ=255)        gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAccccc                                                                                                                                                >  2:94966/1‑149 (MQ=255)                aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATACGtt                                                                                                                                        >  2:574514/1‑149 (MQ=255)                     aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGta                                                                                                                                  <  2:1012998/150‑1 (MQ=255)                             acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACtt                                                                                                                         <  1:303183/151‑1 (MQ=255)                              caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAGTCGAACCCCCCCTACGTTGTCGTATACACACttt                                                                                                                        >  1:786189/1‑151 (MQ=255)                                ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCCACGGTGCCGta                                                                                                                                  >  1:349756/1‑139 (MQ=25)                                ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtat                                                                                                                                 <  2:349784/140‑1 (MQ=255)                                 gTAACATACTCCCACCGCCCGCCAATTTCACTGCCGTACAGACAGATAAAATCCGAAAAAAAAGCTCGCCCTTTCGTACGAGCTCTTCTTTAAATATCGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACTCTTTcc                                                                                                                      <  1:822978/150‑1 (MQ=18)                                   aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAAACCCCCCCACCTTGCCCGATACACACttt                                                                                                                        >  1:491989/1‑146 (MQ=25)                                   aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAg                                                                                                                    <  2:492026/150‑1 (MQ=255)                                      aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTGGAAACCCCACTAAGTTGCCGTTTTTACAATTTCCAGACg                                                                                                                 >  2:117966/1‑150 (MQ=21)                                         cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCCAACGTTGCCGTCTTCACAACTTCCAcgcg                                                                                                                 >  1:666555/1‑147 (MQ=21)                                          tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAAACCCCCGTACGTTTCCGGATACACACTTTCCAGGc                                                                                                                  >  2:580442/1‑145 (MQ=21)                                            cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCACACGTTTCCGTATAAACACTTTTCCGGCGTGCTcc                                                                                                           >  2:621260/1‑150 (MQ=21)                                            cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGAGGGGGGTTTGAAACCCCCCCAAGTTGCCGTATAAACACCTTCCCAGCGTggccc                                                                                                           >  1:72015/1‑146 (MQ=21)                                             aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGAACCCCCCAACAGGTGCCGGATACAAACTTTCCAGGCggg                                                                                                               >  1:874766/1‑143 (MQ=21)                                               ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCa                                                                                                       <  1:538329/151‑1 (MQ=255)                                               ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATACGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCa                                                                                                       <  2:954922/151‑1 (MQ=255)                                                    caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCTACGTTTCCGCATACACACCTTTCAGGCGGGCTCCTTCAGcccc                                                                                                  >  1:998971/1‑149 (MQ=18)                                                    caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                                  <  2:893729/151‑1 (MQ=255)                                                        aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑gg                                                                                             <  1:945987/151‑1 (MQ=255)                                                          gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGGGGGGGGGGATCGTACCCCCCGAACGGTGTCCTGTACACACCTTTCCCGCGTGCTTCTTTTGCCAcaccg                                                                                              >  2:746494/1‑145 (MQ=12)                                                           ttCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACACCCCATAAGTTGCCGTATACACACTTTCTCGGCGTGCTCCTTCATCCACTc                                                                                                >  1:627188/1‑146 (MQ=17)                                                                 gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTc                                                                                     <  2:702480/149‑1 (MQ=255)                                                                   cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAAACCCCCACACGTTGCCGTATAAACACTTTTCCGGGGGGCCCCCTCACCCCCCC‑GGGCA‑ccccccc                                                                                  >  1:475732/1‑145 (MQ=17)                                                                   cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAGGGGGGGGTCGAAACCCCGCGACCTTGCCGTTTACACACCTTTCCCGCGTGCTCCTTCACACACAC‑CGGCA‑CCTCAcc                                                                                  >  2:307044/1‑150 (MQ=17)                                                                    gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATCGAAACCCCGCTACGTTGACGTATACACACTTTTCAGGCGGGGTCCTTCCGCCACTC‑GGGCA‑CCTCAcc                                                                                  >  2:174285/1‑149 (MQ=17)                                                                      acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGGCA‑CCTCACCcaa                                                                               >  2:824929/1‑150 (MQ=255)                                                                       cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTCAAACCCCCCTACGTTGGCCGATAAACACTTTTCCGGCGGGCCCCCTCAGCCCCTC‑CGGCA‑ccccccccaa                                                                               >  1:175772/1‑141 (MQ=9)                                                                           cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGttt                                                                         <  1:206293/151‑1 (MQ=255)                                                                              aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATAAGTTGCCGGATACACACTTTTCAGGCGTGCTCCTTCAGCCACTC‑CGACA‑ACTCACCCAAATGtttt                                                                        >  2:806050/1‑149 (MQ=14)                                                                                aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTa                                                                    <  2:1019004/151‑1 (MQ=255)                                                                                 aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCGCTACGGTGCCGTTTACACACTTTTCAGGCGGGCTCCCTTAGCCACCC‑CGGCA‑CCCCCCCAAATTTTTTTTCt                                                                     >  2:198915/1‑149 (MQ=14)                                                                                 aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGAGTGGTACCCCCAATACACTGGCGAATCAACACTTTTCAATGGGGGTCCTTCCACCACTC‑GGGAA‑CCCCACCAAAATAttttt                                                                       >  2:437662/1‑146 (MQ=9)                                                                                     gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCGTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACTCTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCaaa                                                               <  1:265854/151‑1 (MQ=255)                                                                                        aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCATACGTTGCCGTATACACACCTTCCAGGCGGGCTCCTTCAGCCACTC‑CGACA‑CCTCCCCAAATTGTTTTTCTACCAAAcc                                                             >  2:578465/1‑150 (MQ=14)                                                                                        aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCATAAGTTGCCGTATATACACTTTTCAGGCGTGGTCCTTTAGCCACCC‑GGGCA‑CCCCACCAAAATGTTTTTTTACCAAAcc                                                             >  2:321870/1‑150 (MQ=14)                                                                                          aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCa                                                          <  2:666580/151‑1 (MQ=255)                                                                                               aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGt                                                    <  2:365280/151‑1 (MQ=255)                                                                                                 cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATAAGTTTCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCCCCC‑CGACA‑CCTCACCAAAATGTTTTTCTACCAAACCTCA‑TGGGTg                                                   >  2:370207/1‑150 (MQ=17)                                                                                                  tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTTTCACCCCCCACACGTTTTCGTATACAAACTTTTCAGGGGTGGTCCCTCAGCCCCCC‑CGGCA‑CCCCACCAAAATTTTTTTTTACCAAAACCCC‑Tgggggg                                                  >  1:334814/1‑150 (MQ=2)                                                                                                    gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCTACGTTGCCGTATAAACACTTTCCAGGCGGGGTCCCTCAGCCCCTC‑GGGCA‑CCCCCCCCAATTTTTTTTCTTCCCAAACCCA‑TGGGTgg                                                  >  2:426081/1‑148 (MQ=2)                                                                                                    gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCaa                                               <  2:311867/151‑1 (MQ=255)                                                                                                     cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAAc                                              <  2:17905/151‑1 (MQ=255)                                                                                                         ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACggg                                           <  2:168354/150‑1 (MQ=255)                                                                                                            gTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTTGAAACCCCCCTAATTTTAGTTCTACACACATTCTCGGCGGGCTCCTTCTGCCACCCGGGACA‑CC‑CACCCAAATGTTTTTATAACTAACCTCACAGGGTGG‑GACGGGgcg                                        >  1:797750/1‑150 (MQ=2)                                                                                                              aCGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGAAACCCCGATACGGTGCCGGATAAACACCTTTCAGGCGTGTTCCTTCAGCCCCCC‑CGGCACCCCCCCCAAATGTTTTTTTTACCAAACCCCA‑TGGGGGGAAAgggg                                           >  2:411206/1‑146 (MQ=2)                                                                                                                gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGctact                                   <  1:425341/151‑1 (MQ=255)                                                                                                                  gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGctacta                                  <  1:715104/150‑1 (MQ=255)                                                                                                                  gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACtat                                 <  2:572573/151‑1 (MQ=255)                                                                                                                  gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACtat                                 <  2:275183/151‑1 (MQ=255)                                                                                                                            aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTg                        <  1:621208/150‑1 (MQ=255)                                                                                                                            aaaTATGGCGGTGAGGGGGGGAGTCGAACCCCCCCAACCGTTCCGTATACACACTTTCCCGGGGGGCTTCTTCCGCCCCCC‑CGACA‑CCCCCCCCAAATTTTTTTCTTCCCAACCCCA‑TTGGGGGGGAAGGGGCGCCCCCATAAGGAGGTgg                       >  1:344366/1‑151 (MQ=2)                                                                                                                              atatGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                     <  1:805990/151‑1 (MQ=255)                                                                                                                                    cGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  1:450649/151‑1 (MQ=255)                                                                                                                                            gggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          >  2:148001/1‑148 (MQ=255)                                                                                                                                                 aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:344393/151‑1 (MQ=255)                                                                                                                                                 aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:558658/151‑1 (MQ=255)                                                                                                                                                 aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAAGTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:882857/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                          CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCA‑TGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925756‑926050

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