Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	42.9%	122.0 / 25.2	77	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (5/39);  new base G (33/0);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 7.98e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050                                                                                                                                                      |                                                                                                                                                            gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGa                                                                                                                                                             <  1:567141/150‑1 (MQ=255)       aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGaa                                                                                                                                                       >  1:1082770/1‑150 (MQ=255)         tGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAacc                                                                                                                                                    >  2:925652/1‑151 (MQ=255)            cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTTGGAcccccc                                                                                                                                                 >  2:101744/1‑150 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGAGTTGAAACCCCGATa                                                                                                                                              >  2:1092636/1‑150 (MQ=255)                 aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACg                                                                                                                                            <  1:807674/151‑1 (MQ=255)                      gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGcc                                                                                                                                       <  2:749887/151‑1 (MQ=255)                      gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGcc                                                                                                                                       <  1:690033/151‑1 (MQ=255)                        ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATCGAACCCCCCATACGTTTcc                                                                                                                                       >  1:1062735/1‑149 (MQ=255)                        ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGGGTTGAATCCCCCGTACGTTGCCGt                                                                                                                                     >  2:410336/1‑151 (MQ=37)                        ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGGGTCGGACCCCCCCAACGTTGcc                                                                                                                                       >  1:928087/1‑149 (MQ=37)                         tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGt                                                                                                                                     <  1:227066/150‑1 (MQ=255)                          aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAGGGGGGGGTCGAAACCCCCATACGGTGTCg                                                                                                                                      >  2:1072680/1‑148 (MQ=37)                                tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCTACGTTGCCGtata                                                                                                                                  >  1:139774/1‑146 (MQ=255)                                   caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACttt                                                                                                                        <  2:887329/151‑1 (MQ=255)                                       tAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGTTTTGAcccccccaacgtt                                                                                                                                          >  1:1072643/1‑123 (MQ=25)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATACGGTGCCGGATAAA‑C‑CCttttccgg                                                                                                                   >  1:1038032/1‑143 (MQ=21)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGt                                                                                                                <  1:592045/151‑1 (MQ=255)                                             aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTGAACCCCCCGTACGGTTGCGTATACA‑C‑ACCTTTCAgg                                                                                                                   >  2:79639/1‑146 (MQ=21)                                                  aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAACCCCCCAAACGTTGCCGTATAAA‑C‑ACTTTCCAGGCGTGCTCCtt                                                                                                         >  1:384224/1‑151 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTc                                                                                                        <  1:969405/151‑1 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCATACGTTGCCGTATACA‑C‑ACTTTTCAGGCGGGCTCCtt                                                                                                         >  1:62468/1‑149 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCa                                                                                                       <  1:475479/151‑1 (MQ=255)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGGAACCCCCATACGTTTCCGGATACA‑C‑ACCTTTCAGGCGGGGTccct                                                                                                         >  1:25267/1‑146 (MQ=18)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTTGAAACCCCGCTACGATTCCGCATACA‑C‑ACCTTCCAcgcggggccctt                                                                                                         >  1:557053/1‑140 (MQ=18)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAg                                                                                                      <  1:392221/151‑1 (MQ=255)                                                           ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                                <  1:759886/151‑1 (MQ=255)                                                               gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGa                                                                                             <  1:135069/150‑1 (MQ=255)                                                                   aCTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCt                                                                                       <  1:1042495/151‑1 (MQ=255)                                                                      gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCa                                                                                     <  1:909910/150‑1 (MQ=255)                                                                       ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑cc                                                                                  <  1:1150326/151‑1 (MQ=255)                                                                        cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGTAACCCCCATACGTTGCCGTATAAA‑C‑ACTTTTCCGGCGTGGTCCTTTCGCCCCCCGGAA‑ACCTCA‑cc                                                                                  >  1:587605/1‑150 (MQ=17)                                                                        cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACCCCCGAAACGTTGCCGGATAAA‑C‑ACCTTCCCGGGGGGGTCCTTCCACCACCCGGGC‑ACCCCA‑cc                                                                                  >  2:814649/1‑150 (MQ=17)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAAt                                                                              <  2:990409/151‑1 (MQ=255)                                                                            cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGATTCGAAACCCCCATACGTTGCCGTATAAA‑C‑ACTTTCCAGGCGGGCTCCCTCCGCCCCCCCGAC‑Accccccccaa                                                                                >  2:785460/1‑141 (MQ=14)                                                                             agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTCCAAACCCCCATACGTTGGCGTATACACC‑ACTTT‑CAGGCGTGCTCCTTTAGCCCCCCCGGC‑ACCTCC‑CCcaa                                                                               >  2:715440/1‑148 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCATACGGTGGCGTATACA‑C‑ACCTTTCAGGCGGGCTCCTTCACCCACTCGGAC‑ACCCCA‑CCAAAATGtttt                                                                        >  1:1112162/1‑151 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAACCCCCATAAGATGTCGTATAAA‑C‑ACTTTCCAGGGGGGGTCCTTTCGCCCCTCGGAC‑ACCTCC‑CCCAAATGtttt                                                                        >  2:923269/1‑151 (MQ=14)                                                                                  gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTg                                                                       <  2:693387/151‑1 (MQ=255)                                                                                  gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTg                                                                       <  2:367669/151‑1 (MQ=255)                                                                                   aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTg                                                                       <  2:810721/150‑1 (MQ=255)                                                                                    tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTTAACCCCCCCAACGCTGACGTATAAA‑A‑ACTTTCCAGGGGGGCTCCCTCCGCCACCCGGAG‑ACCCCC‑CCAAATTGttttttt                                                                     >  2:704106/1‑148 (MQ=14)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCa                                                                 <  2:322216/151‑1 (MQ=255)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCTACGTTTCCGTATACA‑C‑ACTTTTCAGGCGTGCTCCTTCCGCCACACCGGC‑ACCTCC‑CCCAATTGTTTTTCTACCaa                                                                >  2:359689/1‑150 (MQ=14)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCaaa                                                               <  1:231265/151‑1 (MQ=255)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCaa                                                                <  2:867690/150‑1 (MQ=255)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCaa                                                                <  2:166615/150‑1 (MQ=255)                                                                                             aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCGATACGTTGCCGGATACA‑C‑ACTTTCCAGGCGGGGTCCCTCAGCCCCTCGGGC‑ACCCCA‑CCCAAATGTTTTTCTACCcaaac                                                             >  1:19290/1‑148 (MQ=14)                                                                                              aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGGAACCCCCATACGTTGCCGTATACA‑C‑ACTTTCCCGGCGTGCTCCTTCAGCCCCTCGGGC‑ACCTCA‑CCAAAATTTTTTGTTACCaca                                                               >  2:813778/1‑145 (MQ=11)                                                                                               aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCTACGTTTCCGTATAAA‑A‑ACTTTCCTGGCGGGGTCCCTCCGCCACCCGGAC‑ACCCCA‑CCCAAATTTTTTTTTACCAAAcc                                                             >  2:837034/1‑148 (MQ=9)                                                                                               aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCa                                                          <  2:244668/151‑1 (MQ=255)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑g                                                       <  1:967163/151‑1 (MQ=255)                                                                                                    aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGt                                                    <  2:887832/151‑1 (MQ=255)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGAACCCCCCCTACGTTTCCCGATACA‑C‑ACTTTCCCGGCGGGCTCCTTCCGCCACCCGGAC‑CCCTCC‑CCCAAATGTTTTTCTACCAAACCCCCT‑gggggg                                                  >  1:756455/1‑148 (MQ=9)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCaa                                               <  2:19295/151‑1 (MQ=255)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCCCCTTCAGCCACCCCGAC‑ACCCCC‑CCAAAATGTTTTTCTACCAAACCTCAT‑GGGGGGCaa                                               >  1:683919/1‑149 (MQ=14)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCaa                                               <  2:592072/149‑1 (MQ=255)                                                                                                             tttCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTTGGAACCCCCCTACGTTTTCGTTTAAA‑CAACTTTCC‑GGCGGGCTCCCTTAACCACCCGGGC‑ACCCCACCCAAATCTTTTTGTTA‑CAAAACTCATAGGGGGGCAAAgg                                            >  2:828370/1‑151 (MQ=2)                                                                                                                          cttcttTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATAAGTTGTCGTATACA‑A‑ACTTTTCAGGGGGGCTCCTTCCGCCACCCCGAC‑ACCCCA‑CCAAAATGTTTTTTTACCAAAACCCAT‑GGGGGGCAACGGGGCGCTACcat                                 >  2:425981/1‑148 (MQ=12)                                                                                                                          cttcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACtata                                <  1:1015518/149‑1 (MQ=255)                                                                                                                            tcttTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCTACGTTGCCGTATACA‑C‑ACTTTCCAGGGGGGGTCCCTCAGCCACTCCGCC‑ACCACC‑CCAAATTGTTTTTTTACCCAACCTCCT‑CGGGGGGAACGGGGCGCTCCTATAggg                             >  2:139488/1‑150 (MQ=9)                                                                                                                            tcttTAAATATGGCGGTGAGGGGGGGGCTGGAACCCCCCGTCCATTGCCATCTACA‑C‑ACTTTTCATGCGGGCTCCTTCTGCCACCCGGACAACCCCA‑CCAACAT‑TTTTGGTACCAAAACCCTT‑GGGGGGCGACGGGGGGCTCCTCTAggg                             >  2:575607/1‑150 (MQ=2)                                                                                                                            tcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAggg                             <  2:258698/150‑1 (MQ=255)                                                                                                                              tttAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTTGCGTATACA‑C‑ACCTTTCCGGCGGGGTCCTTCAACCACCCGGAC‑ACCCCC‑CCCAAATGTTTTTCTACCAAACCTCAT‑GGGGGGGAACGGGGGGCTACTATAggg                             >  1:541600/1‑148 (MQ=12)                                                                                                                              tttAAATATGGCGGTGAGGGGGGGGGTTGTACCCCCGATAAGTTGCTGTATAAA‑C‑ACTTTCCAGGCGGGCTCCTTCAGCCACTCGGCC‑ACCTCC‑CCCAATTGTTTTTCTCCCAAACCTCAT‑GGGGGGCAACGGGGCGCTACTATAggg                             >  2:337786/1‑148 (MQ=18)                                                                                                                               ttAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGt                          <  2:1035598/150‑1 (MQ=255)                                                                                                                                tAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  1:72179/150‑1 (MQ=255)                                                                                                                                 aaaTATGGCGGTGAGGGGGGGGTTTTGACCCCCCATACGGTGTCCTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCCCC‑CCCAATTTTTTTTCTACCAAACCTCCT‑TGGGGGCAACGGGGCGCCACCATAGGGAGtt                         >  2:461336/1‑149 (MQ=17)                                                                                                                                  aaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGa                      <  2:937305/151‑1 (MQ=255)                                                                                                                                     atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                    <  2:111448/150‑1 (MQ=255)                                                                                                                                     atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                   <  2:533622/151‑1 (MQ=255)                                                                                                                                       ggCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCATCGGGGCGCTACTATAGGGAGTTGGAGTaaa                 <  2:1119621/151‑1 (MQ=255)                                                                                                                                       ggCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                  <  1:169581/150‑1 (MQ=255)                                                                                                                                         cGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  2:556976/151‑1 (MQ=255)                                                                                                                                         cGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               >  1:814898/1‑151 (MQ=255)                                                                                                                                          ggTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg              >  1:893072/1‑151 (MQ=255)                                                                                                                                             gAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt            <  2:567531/150‑1 (MQ=255)                                                                                                                                                   gggATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag     <  2:609560/151‑1 (MQ=255)                                                                                                                                                      aTTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  >  1:628910/1‑151 (MQ=255)                                                                                                                                                      aTTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:453608/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                            GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACA‑C‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGAC‑ACCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050

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