Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,882	0	T	G	42.6%	104.3 / 11.9	61	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (4/31);  new base G (26/0);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCA  >  NC_000913/2434736‑2435029                                                                                                                                                    |                                                                                                                                                     gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCGGGCAGGGGaa                                                                                                                                                    >  1:172340/1‑149 (MQ=255)      tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATTAAACCCGGCAGGGGAATACAcc                                                                                                                                              >  1:264687/1‑150 (MQ=255)       atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAGATCAATCACGTCACGGTAATTCAcccc                                                                                                                                            >  1:350333/1‑149 (MQ=255)       atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGGAGATCATACCCGGCAGGGGAATCCAcccc                                                                                                                                            >  2:919487/1‑149 (MQ=255)          tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGTAGATTATACCGGGCAGGGGAATTCACcgc                                                                                                                                            >  1:304690/1‑148 (MQ=255)          tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGAc                                                                                                                                         >  2:307649/1‑151 (MQ=255)            ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCa                                                                                                                                       <  2:362645/151‑1 (MQ=255)                tGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGACCCGa                                                                                                                                     >  2:253782/1‑149 (MQ=255)                 gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATAc                                                                                                                                  <  1:501608/151‑1 (MQ=255)                      gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGccc                                                                                                                             <  2:908840/151‑1 (MQ=255)                       tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGGAGATCATAACCGGCAGGGGAATCCACCCCGACCCGAAACAAccc                                                                                                                             >  2:461404/1‑150 (MQ=255)                        aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCt                                                                                                                            <  1:109752/150‑1 (MQ=255)                        aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGGTTTTGGAGATCATAACGGTCAGGGGAATCCACCCCGACCAGATACAGCCCtt                                                                                                                           >  1:377033/1‑151 (MQ=255)                           cGTCAACTTTTAGGCTGTAACTACGTTGGCTACAGCCTGAAAGGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑TGTTTTGGAATATCATACCGGGCAGGGTAATCCACCGCGAGCAGATAGAGCCCTTc                                                                                                                          <  1:602922/149‑1 (MQ=255)                                tCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATAACGGTCAGGGGAATCCACCCCCCCCAGATAAAGCCCCTCCGC‑tttt                                                                                                                  >  1:527938/1‑151 (MQ=255)                                 cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCGCGACCCGATACAACCCCTCCGC‑tttt                                                                                                                  >  2:247735/1‑150 (MQ=255)                                   ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGc                                                                                                                <  1:652422/150‑1 (MQ=255)                                    tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGcc                                                                                                               <  1:176065/150‑1 (MQ=255)                                      aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTTcc                                                                                                               >  2:443367/1‑148 (MQ=255)                                          tGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAA‑AACATACCCGTGCGGGGAATCCACCCCGACCAGAAACAGGCCCTCCGC‑TTTTGCCGt                                                                                                             >  1:218306/1‑146 (MQ=255)                                           gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTtgct                                                                                                         <  2:167661/149‑1 (MQ=255)                                           gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTtgct                                                                                                         <  1:713351/149‑1 (MQ=255)                                             aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAGATTATAACCGGCAGGGGAATCCAACCCGACCAGAAAAAGCCCCTCCCC‑CTTTTCCCTTGCTgcc                                                                                                      >  1:908437/1‑150 (MQ=255)                                              aCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCGC‑TTTTGCCGTTGCTgccgc                                                                                                    >  2:309000/1‑151 (MQ=255)                                                   gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCg                                                                                               <  1:777748/151‑1 (MQ=255)                                                   gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCg                                                                                               <  2:390617/151‑1 (MQ=255)                                                    ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCCCGACCAGATACAACCCTTCCCC‑TTTTTCCGTTGCTgccccc                                                                                                   >  2:491967/1‑146 (MQ=255)                                                      ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGtt                                                                                             <  1:559761/150‑1 (MQ=255)                                                      ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTc                                                                                            <  1:468306/151‑1 (MQ=255)                                                       gCTACAGCCTGAAATGTTCCGAGCGTTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCt                                                                                           <  2:375183/151‑1 (MQ=255)                                                          aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTg                                                                                          <  1:505153/149‑1 (MQ=255)                                                               cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAAAGCAAACCGGGCAGGGGAATTCACCGCGGCCAGAAACAACCCCTCCGC‑CTTTTCCGCTGCCGCCCCCCACGCTCTTTcc                                                                                       >  2:149802/1‑147 (MQ=255)                                                                     tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCGCGACGAGATACAGCCCTTCCGC‑CTTTTCCGTTGCTGCCGCCCGCGCTCTGTCCTTTtccgccc                                                                             >  1:114941/1‑150 (MQ=255)                                                                      gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAg                                                                            <  1:607558/151‑1 (MQ=255)                                                                      gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAg                                                                            <  1:43917/151‑1 (MQ=255)                                                                      gTTCCGACCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                             <  1:295951/150‑1 (MQ=255)                                                                              cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAGGTCATAACGGCCAGGGGAATTTACCGCGCCCAGCCACAACCCTTCCGCGTTTTTCCCTT‑CTGCCCCCCGCGGTCTGGCCttttc                                                                                  >  1:16805/1‑135 (MQ=255)                                                                               attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAGAACATTACCGACAGGGGAAATCACCGCGACCCGATAAAGGCCTTCCCC‑CTTTGCCGTTTTTGGCGCCAGCCGTCCTTCCCTTTCCCCCCGCAACCCTg                                                                    >  1:585406/1‑150 (MQ=255)                                                                               attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGC‑CTTTGCCGTTGCTGCCCCCCGCGTTCTGTCCTTTGCCGCCAGCaa                                                                         >  2:1024595/1‑145 (MQ=255)                                                                                  attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTAt                                                                <  1:655414/151‑1 (MQ=255)                                                                                    tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATc                                                               <  2:864829/150‑1 (MQ=255)                                                                                          tAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATAAAGCCCTTCCGC‑TTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCCGCAACTCTGCTATcc                                                              >  2:641112/1‑145 (MQ=255)                                                                                            gTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCGGTCAGGGGAAACCACCGCGGCCAGATACAGCCCCTCCCC‑TTTTTCCGTTGCTGCCCCCCGCGTTCTGTCCTTTTCCCCCAGCCACTCTGCCATCCCGCTCTcc                                                      >  1:96096/1‑151 (MQ=255)                                                                                            gTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                                      <  2:252054/151‑1 (MQ=255)                                                                                             tCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGGAGATGCTACCGGTCAGGGGAATCTACCCCGCCCCAAGACAACCCCTCCGC‑CTTTGCCGTTGCTGCCCCCAGCGGTCTGTTCTTTGTCCCCAAGAACTCTGCTATCCatcac                                                         >  2:169439/1‑143 (MQ=255)                                                                                             tCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGAAGATCATACCCGGCAAGGGAATCCACCCCGAACCGATACAGCCCTTCCGC‑TTTTTCCGTTGCTGCCGCCCCCGGTTTGTTCTTTTCCGCCCGCAAATCTGCTAt                                                                >  1:355401/1‑140 (MQ=255)                                                                                              ccGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTTGTAGTTTCTACCGGGTCGGGGAATGCACACCGACCAGATAAAGCCCCTCCCC‑TTTTTCCGCGGCTGCCCCCAACGTTTTGTTCTTTCCCGCCCGCCACTTTGCCCTCCAGCTCTcc                                                      >  1:162542/1‑149 (MQ=255)                                                                                                 ccAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                                                 <  1:194126/151‑1 (MQ=255)                                                                                                        tAGCGGGCCCATTGGC‑GGTTTTTGAAGATCAAACCGGTCAGGGGAATCCACCGCGCCCCGATACAGCCCTTCCCC‑TTTTTCCGTTTCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTTTCCCGCTCTCCCGCCGGTtgtg                                           >  1:911888/1‑150 (MQ=255)                                                                                                           cGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                                       <  1:43930/151‑1 (MQ=255)                                                                                                             ggCCCATTGGC‑GGTTTTTGGAGATTATAACGGTTAGGGTAATACACCGCCACCAGATACAGCCCTTCCGC‑TTTTTCCGTTGCTGGCCCCCGCGTTCTGTCCTTTGCCGCCCGCAAACCTGCTCTCCAGCTCTCCGGCCGGGTGTGGGCAcc                                     >  1:156454/1‑151 (MQ=255)                                                                                                             ggCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCAcc                                     <  2:668630/151‑1 (MQ=255)                                                                                                                  aTTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACtt                                <  1:575865/151‑1 (MQ=255)                                                                                                                    tGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCCCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                              >  1:810294/1‑151 (MQ=255)                                                                                                                     ggC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                              <  2:886164/150‑1 (MQ=255)                                                                                                                                  aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                 <  1:43398/151‑1 (MQ=255)                                                                                                                                      cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat             <  2:784123/151‑1 (MQ=255)                                                                                                                                      cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat             <  1:589828/151‑1 (MQ=255)                                                                                                                                          ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc         >  1:674177/1‑151 (MQ=255)                                                                                                                                          ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc         <  2:943957/151‑1 (MQ=255)                                                                                                                                                 gggTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCa  <  1:48052/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGC‑GGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGC‑TTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCA  >  NC_000913/2434736‑2435029

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