Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,581	0	A	G	47.5%	93.3 / 18.6	61	noncoding (72/88 nt)	serX	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (3/29);  new base G (29/0);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 2.70e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAA  >  NC_000913/1097445‑1097729                                                                                                                                           |                                                                                                                                                         cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                            <  2:757891/151‑1 (MQ=38)         gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGcc                                                                                                                                    <  1:782819/151‑1 (MQ=255)           ccAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGt                                                                                                                                  <  1:268808/151‑1 (MQ=255)              aCTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGtata                                                                                                                               <  1:991673/151‑1 (MQ=255)                tGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑AGGGGGGGTTCGAACCCCCCCTACATTTCCCTATac                                                                                                                              >  2:38861/1‑150 (MQ=11)                     gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑Ctt                                                                                                                      <  1:363646/151‑1 (MQ=255)                           aGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAACCCCCCCTACCTTTCCGGATACAC‑A‑CTTTTCAgg                                                                                                                >  1:603991/1‑151 (MQ=14)                              tGTCTAACCACCGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAgg                                                                                                                <  2:95589/148‑1 (MQ=255)                               gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGt                                                                                                             <  1:574829/150‑1 (MQ=255)                                 cGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGATCGAACCCCCCATAAGTTGCCGGATACAC‑A‑CTTTCCAcgcgctc                                                                                                           >  1:102289/1‑147 (MQ=17)                                  gAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GT‑GGGGGGGGTGCGGGCGACCGGGACGGTGTCGGATACAC‑A‑ATTTTTCGGCGTg                                                                                                            >  2:483462/1‑148 (MQ=1)                                     ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCtt                                                                                                      <  2:137906/151‑1 (MQ=255)                                       ccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCa                                                                                                    <  2:1023497/151‑1 (MQ=255)                                          gTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑TA‑GGGGGGGGTTCGAAACCCCCATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGGTCCCTCCGcc                                                                                                 >  1:1111633/1‑151 (MQ=18)                                                gCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAACCCCCCAAACGTTGCCGCATACAC‑A‑CTTTTCACGCGTGCTCCCTCAGcccccc                                                                                             >  2:354259/1‑145 (MQ=17)                                                gCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑gg                                                                                          <  2:12966/151‑1 (MQ=255)                                                  ttCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGAACCCCCCATACGTTTCCGTATACAC‑A‑CTTTCCACGCGGGCTCCTTCAGCCCCTC‑gg                                                                                          >  1:975035/1‑149 (MQ=17)                                                   tctcATCCCCCCGGTGTGTGCAATATACGAAAAACAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGa                                                                                         <  2:125786/149‑1 (MQ=255)                                                   tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑C‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGa                                                                                         <  1:1046505/149‑1 (MQ=255)                                                   tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGac                                                                                        <  1:251923/150‑1 (MQ=255)                                                   tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGa                                                                                         <  1:1077568/149‑1 (MQ=255)                                                     tcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GC‑GGGGGGGGTTTGCAACCCCCCCACCTTGGCGGATACAC‑C‑CTTTCCAAGACTGCTCCTTCTACCACCC‑CGACAcc                                                                                     >  1:903047/1‑151 (MQ=14)                                                      cATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGacac                                                                                      <  1:604299/149‑1 (MQ=255)                                                        tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAACCCCCCGCACGTTTCCCGATAAAC‑A‑CTTTTCCGGCGTGCTCCTTTAGCCACTC‑GGACAcccc                                                                                   >  1:939681/1‑148 (MQ=17)                                                          cccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑AGGGGGGAGTCGGACCCCCGATACGTTTTCGTATATAC‑A‑CTTTCCCGGGGTGCTCCCTTAGCCACTC‑GGGCACCCCA‑cc                                                                               >  2:189256/1‑151 (MQ=17)                                                              cGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑AGGGGGGGATCGGACACCCCATACGTTGTCCTATACAC‑C‑CTTTTCCGGGGGGGTCCTTCCACCCCCC‑AGCCACCTCA‑CCcaa                                                                            >  2:884569/1‑150 (MQ=9)                                                                gtgtgtGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTg                                                                         <  1:805076/151‑1 (MQ=255)                                                                 tgtgtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAAACCCCCATACGTTGCCGTATATAC‑A‑CATTTCAGGCGTGCCCCCTCAGCCACTC‑GGGCACCCCACCCAAAttt                                                                         >  1:9309/1‑150 (MQ=14)                                                                    gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTCAACCCCCCAGACGTTGCCGTATAAAC‑A‑CTTTTCCGGGGTGCTCCCTCCGCCACCC‑GGACACCTCC‑CCAAAtttt                                                                        >  1:16910/1‑146 (MQ=14)                                                                    gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCg                                                                      <  1:969107/150‑1 (MQ=255)                                                                       cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAACCCCCCATACGTTTCCGTATACAC‑A‑CTTTTCCGGCGTGGCCCTTCAGCCACCC‑CGGCACCCCC‑CCCAAATTTCGTTcc                                                                  >  1:657527/1‑151 (MQ=14)                                                                        aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCt                                                                 <  1:1086024/151‑1 (MQ=255)                                                                        aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCt                                                                 <  1:150574/151‑1 (MQ=255)                                                                             aCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GT‑GGGGGGGGTTTGAACCCCCCACACGTTGCCGGTTAAACAA‑CTTTCCA‑GCGTGCTCCTTTTACCCCCC‑CGACACCCCA‑CCAAATTATCGTTCCTGTCtt                                                            >  1:560351/1‑151 (MQ=9)                                                                             aCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCtt                                                            <  1:1003976/151‑1 (MQ=255)                                                                              cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑GGGGGGGGTTTGTAACCCCGCTACGTTGCCGTATACAA‑ACCTTTCCA‑GCGGGGTCCCTTCGCCACCC‑CGACACCCCA‑CCAAAATttccttcc                                                                  >  2:940419/1‑144 (MQ=2)                                                                               gAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑TA‑GGGGGGGGTTCCAACCCCCCCTAAGTTGCCCGATTAAA‑A‑CCTTCCAAGGGGGGGCCATCAACCACCC‑GGAACCCCCC‑CCAAATTTTTGTTcccgacct                                                            >  2:238399/1‑143 (MQ=2)                                                                                  aaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAAACCCCGCAACGTTTTCGTATACAC‑A‑ACTTCCCGGCGTGCCCCTTTCCCCACCC‑CGACACCTCA‑CCAAATTTTCGTTCCTGTCTTTCttg                                                       >  2:1018128/1‑149 (MQ=14)                                                                                     aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTTGGAACCCCC‑CTACGTTTCCGTATACAC‑A‑CTTTTCACGGGGGGTCCTTCAGCCACCC‑GGACCCCCCC‑CCCAAATTTCGTTCCTTTCTTTCTggga                                                     >  2:1075708/1‑148 (MQ=9)                                                                                        gCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACggg                                                 <  2:751849/151‑1 (MQ=255)                                                                                           cGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTTCGACCCCCCCCATCTTTTCC‑TATAAAC‑A‑CATTTCAGGGGGGCTCCTTCCGCCACCCGGGCCACCCCC‑CCAAATT‑TTGTTTCTGTCTTGTTGGGAAggg                                                 >  2:468852/1‑148 (MQ=2)                                                                                           cGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAACCCCCGCTACGTTGCCGTATACAC‑A‑CTTTTCCGGCGTGCTCCTTCCGCCCCTC‑GGGCACCTCA‑CCCAAATTTCGTTTCTGTCTTTCTGGGACggg                                                 >  1:298779/1‑148 (MQ=14)                                                                                            gTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCt                                             <  2:1093470/151‑1 (MQ=255)                                                                                              aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGAACCCCCCCCACGTTTGCCTATAAAC‑A‑CTTTACTGGCGGGGTCCTTCTGCCACCC‑GGGAACCCCA‑CCAAATTGTTGGTTCTGGCTTTCTGGAAAGGgggc                                              >  2:220845/1‑148 (MQ=2)                                                                                               cTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGAACCCCCCATACGTTTCCGTATAAAC‑A‑CTTTTCCGGCGTGGTCCCTCAGCCACCC‑GGACAACCCA‑CCAAAATTTCGTTTCTTTTTTTTTGGAACGGGGGCCaa                                           >  2:702666/1‑150 (MQ=2)                                                                                                tttCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGGACCCCCGATACGTTGCCGTATACAA‑A‑CTTTTCCGGCGGGCTCCTTCAGCCACTC‑GGGCACCTCA‑CCAAATTTTTGTTCCTTTCTTGCTGGAACGGGCGCTAAtt                                         >  1:392791/1‑151 (MQ=17)                                                                                                tttCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAACCCCCCATACCTTTCCGTATAAAC‑A‑CCTTCCAGGGGGGCTCCTTCCCCCACCC‑CGACACCCCA‑CCAAAATTTTGTTCCTTCCTTTCTTGAAAGGGGGCTaaat                                         >  2:531106/1‑149 (MQ=2)                                                                                                 ttCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAAt                                          <  1:738951/149‑1 (MQ=255)                                                                                                      aCGAGCTCTTCTTTAAATATGGCGGTGGA‑GGGGGGGGTT‑GAAAACCCCATACGTTGGAGTATAAAC‑A‑CATTTCAGGCGGGGGCCCTCTGACACAC‑GGGGACCACA‑CCAAATTTTTTTTTCTTTCTTTTTGGAAGGGGGGGCAATTTTggg                                    >  2:644934/1‑150 (MQ=0)                                                                                                      aCGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAggg                                    <  2:127527/150‑1 (MQ=255)                                                                                                       cGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGCACCCCCCCTACGTTTCCGTATACAC‑A‑CTTTCCCGGCGTGCTCCTTCCGCCACCC‑CGGCACCTCC‑CCAAATTTTTGTTCCTGTCTTTCTTGGAAGGGCGCCAATTTAggg                                    >  1:951449/1‑149 (MQ=14)                                                                                                       cGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGaa                                  <  2:825075/151‑1 (MQ=255)                                                                                                           ctctTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAACCCCCCATACGTTTCCGTATTAAC‑A‑CTTTTCAGGCGGGCTCCTTCAGCCCCTC‑GGACACCCCC‑CCCAATTGTCGTTCCTGTTTTGTTTGTACGGGCGCTAATTTTGGGGAAt                                >  2:549064/1‑149 (MQ=14)                                                                                                           ctctTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAAACCCCCCCACGTTTCCCTATAAAC‑A‑CTTTTCAGGCGGGGTCCTTCAGCCACTC‑GGACACCTCC‑CCAAAATTTCGGTCCTTGCTTTCTTGAACGGGCGCTAATTTTGGGaaa                                 >  2:582054/1‑148 (MQ=9)                                                                                                                  ttAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGAc                          <  1:643942/148‑1 (MQ=255)                                                                                                                     aaTATGGCGGT‑GA‑GGGGGGGGTTAGGACCCCCGATACGTTGCCGTATTAAC‑A‑CTTTTCAGGGGGGGTCCTTCAGCCACTC‑GGACACCTCC‑CCAAAATTTCGTTCCTGTTCTTCTGGGAAGGGGGCTAATTTTGGGAAATCATGACCTTAgg                    >  2:1121422/1‑151 (MQ=12)                                                                                                                                  a‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTCTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttt        >  2:434315/1‑151 (MQ=255)                                                                                                                                     ggggggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttttt      <  1:1044949/151‑1 (MQ=255)                                                                                                                                        gggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaa   <  2:560336/151‑1 (MQ=255)                                                                                                                                         ggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaa  <  2:58565/151‑1 (MQ=255)                                                                                                                                          gATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaa  >  2:978460/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                                         CACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTC‑GGACACCTCA‑CCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAA  >  NC_000913/1097445‑1097729

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