Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	50.0%	78.8 / 21.1	64	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (6/26);  new base G (32/0);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 3.01e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAA  >  NC_000913/925751‑926043                                                                                                                                                        |                                                                                                                                                      gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGa                                                                                                                                                       <  2:109401/150‑1 (MQ=255)   gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTc                                                                                                                                                    <  1:400533/151‑1 (MQ=255)     tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGtt                                                                                                                                                     >  2:179223/1‑148 (MQ=255)     tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑g                                                                                                                                                  <  1:32106/150‑1 (MQ=255)            cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTC‑T‑GAccccc                                                                                                                                          >  2:176517/1‑150 (MQ=255)             gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGa                                                                                                                                        <  2:610092/151‑1 (MQ=255)                    cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTg                                                                                                                                 <  1:359674/151‑1 (MQ=255)                      gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTG‑G‑AAcccccccacgttgcc                                                                                                                                >  1:620586/1‑141 (MQ=37)                        ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑AGGGGGGGGTC‑G‑GACCCCCGACACGTTTcc                                                                                                                               >  2:210366/1‑149 (MQ=37)                         tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑AAccccccccacgtttcc                                                                                                                               >  1:116774/1‑137 (MQ=37)                         tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCg                                                                                                                              <  2:745738/149‑1 (MQ=255)                          aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑AAACCCCCCAACGTTGTCg                                                                                                                              >  1:463686/1‑148 (MQ=37)                          aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTC‑G‑AACCCCCCATACGtt                                                                                                                                  >  1:871431/1‑144 (MQ=255)                          aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGta                                                                                                                            <  1:220766/150‑1 (MQ=255)                           gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGtata                                                                                                                          <  1:996817/151‑1 (MQ=255)                               ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCCAGACGTTGcc                                                                                                                               >  2:642130/1‑142 (MQ=255)                                  acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACtt                                                                                                                  <  1:627533/151‑1 (MQ=255)                                             aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTg                                                                                                       <  1:172832/150‑1 (MQ=255)                                                 cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCt                                                                                                  <  2:570795/151‑1 (MQ=255)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTC‑G‑AAACCCCGATACGTTGGCGTATACAA‑ACTTTCCAGGC‑GTGCTCCtt                                                                                                 >  2:364103/1‑150 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑AAACCCCCCTACGTTTCCGTATACAC‑ACCTTTCCGGC‑GGGCTCCtt                                                                                                 >  2:199162/1‑149 (MQ=18)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAg                                                                                              <  2:490858/151‑1 (MQ=255)                                                       cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGATTC‑G‑AACCCCCCATACGTTGCCGTATACAC‑ACCTTTCCGGC‑GTGCTCCTTCAGcc                                                                                            >  2:238735/1‑151 (MQ=255)                                                        ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTC‑GAAACCCCCCAACCGTT‑CCGTATATAC‑ACTTTCCATGC‑GTGCGcccttcgcc                                                                                            >  1:997869/1‑143 (MQ=18)                                                          accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑Ct                                                                                        <  2:796188/151‑1 (MQ=255)                                                              tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTC‑C‑GACCCCCCCCACGTTGCCGGTTACAC‑ACCTTTCCGGC‑GTGCTCCTTctacccctc                                                                                        >  1:503509/1‑140 (MQ=18)                                                                ttCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGCGA‑GGGGGGGGTTCTT‑AACCCCCG‑TACGTCGCCGCATAAAC‑ACATTCCCGGC‑GGGCTCCTTCAGc                                                                                             >  2:233662/1‑141 (MQ=18)                                                                   aCTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑AACCCCCCATACGTTGCCGTTTACAC‑ACTTTTCAGGC‑GTGGTCCCTCCGcccccccggcaccc                                                                                >  1:810164/1‑139 (MQ=17)                                                                          tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGAGA‑GGGGGGGGTTC‑G‑GAACCCCGGTACGTTTCCGTATACAC‑ACTTTCCAAGC‑GTGCTCCTTCAGCCA‑CC‑CCGGCACCTCAcc                                                                          >  1:796699/1‑148 (MQ=17)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGtttt                                                                <  1:455094/151‑1 (MQ=255)                                                                                  gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTg                                                               <  2:483857/151‑1 (MQ=255)                                                                                     aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑AAACCCCCATACGTTTCCGTTTAAAA‑ACCTTCCCGGC‑GGGGTCCTTCCGCCA‑CC‑CCGGCACCCCCCCCAAAtttttt                                                                >  1:275432/1‑147 (MQ=9)                                                                                      aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTC‑G‑GAACCCCCATACGTTTGCGCATACAC‑ACTTTTCAGGC‑GGGCCCCCTTAGCCC‑CC‑CCGGCACCCCACCCAAATGtttt                                                                >  2:390130/1‑146 (MQ=9)                                                                                         tGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTC‑G‑AAACCCCCATACGTTGGCGTATACAC‑ACTTTTCAGGC‑GGGCTCCCTTAGCCC‑CC‑CGGACACCCCACCAAAATGTTTTTCTAcc                                                          >  2:627435/1‑149 (MQ=14)                                                                                         tGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTT‑G‑AACCCCCGCTACGTTTCCGTATATAC‑ACTTTTCAAGC‑GTGCTCCTTCAGCCC‑CC‑CCGGCACCTCACCAAAATGTTTTTCTTcc                                                          >  1:1174376/1‑149 (MQ=11)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATg                                                <  2:713891/151‑1 (MQ=255)                                                                                                  aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTT‑G‑GACCCCCCCTACGGTGCCGGATAAAA‑ACCTTCCCGGCGGGGCTCCTTC‑GCCACCC‑CGGCCCCCCCCCAAAATTTTTTTTTTACCAAcccccttg                                                >  1:942683/1‑146 (MQ=2)                                                                                                   aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTC‑G‑AACCCCCCCTACGTTGCCGTATACAC‑ACTTTCCCGGC‑GGGCTCCTTCCGCCC‑CC‑CGGACACCCCACCAAAATTTTTTTTTACCCAACCCCCTggg                                              >  2:1088390/1‑151 (MQ=11)                                                                                                    aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTC‑G‑AACCCCCGGAACGGTGCCGTATACAC‑ACTTTCCAGGC‑GGGGTCCCTCAGCCA‑AC‑CGGGCAACCCAACAAAATGTTTTTGTTCCAAAcaccctggg                                              >  1:362542/1‑142 (MQ=2)                                                                                                    aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGt                                             <  2:997128/151‑1 (MQ=255)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGG‑GGGGGGGATTG‑T‑AACCCCC‑ATACGTTTGCGTATTCAC‑ACCTTTCTGGG‑GTGGTCCTTCCGCCA‑CC‑CGGACACCCCCCCCAAATGTTTTTCTTCCCAAcaccacagggggga                                         >  1:570686/1‑137 (MQ=2)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTC‑G‑AACCCCCCATACGTTGCCCGATACAC‑ACTTTTCAGGC‑GTGCTCCTTCCGCCA‑CT‑CGGGCACCCCCCCCAAATGTTTTTTTACCAAACCTCATtggtgg                                           >  2:867953/1‑148 (MQ=14)                                                                                                               tCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑GAACCCCCCCACGTTTCCGTATACAA‑ACCTTTCCGGC‑GGGCTCCCTCAGCCA‑CT‑CCGGCACCCCACCCAAATTTTTTTCTACCAAAACTCATGGGGGGGAACgggg                                   >  1:111758/1‑150 (MQ=9)                                                                                                               tCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTT‑G‑AACCCCCCCTACGTTTCCGTATAAAC‑ACTTTTCAGGC‑GTGGTCCTTCCGCCA‑CC‑CGGACACCCCCCCCAAATTTTTTTTTACCCAACCCCCTTGGTGGGAACgggg                                   >  2:483735/1‑150 (MQ=2)                                                                                                               tCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgggg                                   <  1:1038725/150‑1 (MQ=255)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGC‑GGTTA‑GGGGGGGGTTC‑G‑GACCCCCCCAACGTTTGCGTATACAC‑ACCTTCCAGGG‑GTGGTCCTTCCCCCA‑CC‑CGGACACCCCACCCAAATTTTTTTCTACCAAACACCAAGGGGGGGAACGGGgcg                                 >  1:853010/1‑151 (MQ=9)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgc                                  <  1:1114189/150‑1 (MQ=255)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGGGTC‑G‑AAACCCCCCCTCGCTTGCATCTAAACAACTTTCTA‑GC‑GCGGTCCTTCTTCCA‑CC‑CCGGCAACTCCCCAAAATGTTTTGCTTTCCAAACTCATAGGGGGGAAAggggggg                                >  1:1060117/1‑148 (MQ=2)                                                                                                                    cGAGCTCTTCTTTAAATATGGC‑GGTTA‑GGGGGGGGGTC‑G‑AAACCCCCCCACCTTTCCGTATAAAC‑ACATTTCAGGC‑GGGCTCCTTCAGCCA‑CC‑CCGGCACCCCACCCAAATGTTTTTTTAACAAAACCCCTGGGTGGGAACgggg                                   >  1:389669/1‑145 (MQ=2)                                                                                                                    cGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctac                             <  1:936631/151‑1 (MQ=255)                                                                                                                     gAGCTCTTCTTTAAATATGGC‑GGTGT‑GGGGGGGGGTC‑G‑AAACCCCCACACCTTGCCGTATAAAC‑ACTTTCCAAGG‑GTGCTCCTTCAGCCC‑CC‑CGGGCACCCCACCAAAATTTTTTTTTACCAAAACTCATGGGGGGCAACGGGGGGccac                             >  1:803486/1‑150 (MQ=9)                                                                                                                     gAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTC‑G‑AACCCCCCCTACGGTGCCGTATACAC‑ACCTTCCAAGC‑GCGCCCCTTCAGCCA‑CC‑CGGGCACCCCACCAAAATGTTTTTCTTCCAAACCTCATTGGGGGCAACGGGgggc                                >  2:587330/1‑147 (MQ=14)                                                                                                                      aGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑CAACCCCCATACGTTTGCGTATAAAC‑ACCTTCCAGGC‑GTGGTCCTTTCTCCA‑CTCCGGCAACCCCCCCAAATT‑TTTTTTTACCCAAACTCATGGGGGGGGACgggggg                                 >  1:780851/1‑143 (MQ=2)                                                                                                                             cttTAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑AACCCCCCCTACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CC‑CCGACACCTCACCAAATTGTTTTTCTACCAAACCTCCTGGGTGGCAAAGGGGCGCTACTATAggg                      >  2:749816/1‑149 (MQ=18)                                                                                                                              tttAAATATGGC‑GGTGA‑GGGGGGGGGTC‑G‑AAACCCCCATACGTTTTCGTATAAAC‑ACTTTTCAGGC‑GGGCTCCTTCAGCCA‑CT‑CCGACACCCCCCCAAATTGTTTTTCTACCAAACCTCCTGGGGGGCGACGGGGCGCTACTATAggggg                    >  1:193778/1‑148 (MQ=17)                                                                                                                               ttAAATATGGC‑GGTGA‑GGGGGGGGTTT‑G‑AACCCCCCCTACGCTGCCGTCTACAC‑ACTTTTCAGGC‑GGGGTCCCTCCGCCC‑CA‑CGGCGACCCCACCCAAATGTTTTTCTACCCAACCCCCTGGGTGGCCACGGGGCGCCACTATAgggggg                   >  1:942595/1‑147 (MQ=12)                                                                                                                               ttAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAg                    <  2:589515/149‑1 (MQ=255)                                                                                                                               ttAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt                   <  1:303960/150‑1 (MQ=255)                                                                                                                                  aaTATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGa               <  1:38072/151‑1 (MQ=255)                                                                                                                                        gC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa         <  2:1019612/151‑1 (MQ=255)                                                                                                                                           ggTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg       >  2:861536/1‑151 (MQ=255)                                                                                                                                           ggTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg       >  2:41304/1‑151 (MQ=255)                                                                                                                                              gA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt     <  1:269696/150‑1 (MQ=255)                                                                                                                                                 gggggggATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa  <  2:707099/151‑1 (MQ=255)                                                                                                                                                 gggggggATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa   <  2:605069/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTC‑G‑AACCCCCGATACGTTGCCGTATACAC‑ACTTTCCAGGC‑GTGCTCCTTCAGCCA‑CT‑CGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAA  >  NC_000913/925751‑926043

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