Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,518,197	C→G	6.2%	noncoding (35/76 nt)	alaW ←	tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,518,197	0	C	G	6.2%	284.9 / 8.6	81	noncoding (35/76 nt)	alaW	tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base C (45/31);  new base G (2/3);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

AAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAC  >  NC_000913/2518051‑2518343                                                                                                                                                   |                                                                                                                                                   aaGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                                  <  4:119554/149‑1 (MQ=255)   gCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatg                                                                                                                                                <  1:515778/149‑1 (MQ=255)        aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                           <  4:199887/149‑2 (MQ=255)        aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                           <  1:521919/149‑2 (MQ=255)        aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc                                                                                                                                           <  3:124902/149‑1 (MQ=255)           gAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGct                                                                                                                                        >  3:115674/1‑149 (MQ=255)            aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctc                                                                                                                                       <  4:775027/149‑1 (MQ=255)                 cTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGCGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCa                                                                                                                                  >  4:193244/1‑149 (MQ=255)                      ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                             <  1:856764/149‑1 (MQ=255)                      ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                             <  1:700480/149‑1 (MQ=255)                       tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                            >  3:538613/1‑149 (MQ=255)                       gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                            >  1:197402/2‑149 (MQ=255)                       gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                            >  4:788816/2‑149 (MQ=255)                             atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAg                                                                                                                      >  2:596750/1‑149 (MQ=255)                              tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCTCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCACTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc                                                                                                                     <  3:603033/149‑1 (MQ=255)                              tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc                                                                                                                     >  3:727343/1‑149 (MQ=255)                              tgccacgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc                                                                                                                     <  1:630882/149‑1 (MQ=255)                                 catgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcccc                                                                                                                  >  3:820719/1‑149 (MQ=255)                                 catgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcccc                                                                                                                  <  3:595438/149‑1 (MQ=255)                                      aaGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCAc                                                                                                             >  3:379202/1‑149 (MQ=255)                                      aaGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGGTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGGTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCAc                                                                                                             >  3:161244/1‑149 (MQ=255)                                        gcgcTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGt                                                                                                           >  2:750086/1‑149 (MQ=255)                                           ctctCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTaaa                                                                                                        <  1:215871/149‑1 (MQ=255)                                            tctcCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAg                                                                                                       >  3:414146/1‑149 (MQ=255)                                            tctcCCAGCTGAGCTATAGCCCCACGATGAGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAg                                                                                                       >  3:525339/1‑149 (MQ=255)                                              tcCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCt                                                                                                     >  3:200576/1‑149 (MQ=255)                                                ccAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTg                                                                                                   >  1:772469/1‑149 (MQ=255)                                                ccAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTg                                                                                                   >  2:595319/1‑149 (MQ=255)                                                ccAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCACTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTg                                                                                                   >  4:705277/1‑149 (MQ=255)                                                  aGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTc                                                                                                 <  3:193244/149‑1 (MQ=255)                                                   gCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCCATAGCCCCGTCACGTAAAGCTTGTCg                                                                                                >  3:440958/1‑149 (MQ=255)                                                      gAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGt                                                                                             >  2:261842/1‑149 (MQ=255)                                                         cTATAGCCCCACGCTGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCGAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCTAGCTGAGCTATAGCCCCGGCGCGTAAAGCTTGTCGAGTTGa                                                                                          >  2:1152/1‑149 (MQ=25)                                                         cTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCACTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGa                                                                                          <  3:705277/149‑1 (MQ=255)                                                              gCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggc                                                                                     <  1:595319/149‑1 (MQ=255)                                                               ccccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAGTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcg                                                                                    >  2:780582/1‑149 (MQ=255)                                                                cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcgg                                                                                   >  1:53714/1‑149 (MQ=255)                                                                cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcgg                                                                                   <  4:200576/149‑1 (MQ=255)                                                                 ccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACgg                                                                                       <  2:617940/144‑1 (MQ=255)                                                                 ccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACgg                                                                                       >  1:617940/1‑144 (MQ=255)                                                                 ccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcggc                                                                                  >  1:836082/1‑149 (MQ=255)                                                                 ccACGATGCGTTTACGTACCAAGTATGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcggc                                                                                  >  3:702692/1‑149 (MQ=255)                                                                  cACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca                                                                                 >  4:521066/1‑149 (MQ=255)                                                                   aCGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGcat                                                                                <  3:506571/149‑1 (MQ=255)                                                                     gATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGcatca                                                                              <  1:225275/149‑1 (MQ=255)                                                                       tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCata                                                                            <  4:87130/149‑1 (MQ=255)                                                                        gCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGGGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCatat                                                                           <  4:115674/149‑1 (MQ=255)                                                                          gTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGa                                                                         >  4:586775/1‑149 (MQ=255)                                                                           tttACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGaa                                                                        <  4:820719/149‑1 (MQ=255)                                                                             tACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCATGTCGAGTTGACGGGCGGCATCATATGAAtt                                                                      <  2:332385/149‑1 (MQ=255)                                                                                  aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGcc                                                                 >  4:234169/1‑149 (MQ=255)                                                                                  aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGcc                                                                 <  3:469844/149‑1 (MQ=255)                                                                                   ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGccc                                                                >  2:585491/1‑149 (MQ=255)                                                                                   ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGccc                                                                >  4:135595/1‑149 (MQ=255)                                                                                   ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCCAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGccc                                                                >  1:538698/1‑149 (MQ=255)                                                                                       gTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGGAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAAt                                                            <  4:497846/149‑1 (MQ=255)                                                                                       gTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAAt                                                            >  3:355545/1‑149 (MQ=255)                                                                                       gTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCATGTCGAGTCGACGGGCGGCATCATATGAATTCCGCCCGAAt                                                            >  4:719640/1‑149 (MQ=255)                                                                                              gggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCaa                                                     >  4:51602/1‑149 (MQ=255)                                                                                               ggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAAc                                                    <  3:151172/149‑1 (MQ=255)                                                                                                 tGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACgg                                                  >  2:883245/1‑149 (MQ=255)                                                                                                  gCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGc                                                 <  1:418896/149‑1 (MQ=255)                                                                                                    aaaaTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCaa                                               >  1:274069/1‑149 (MQ=255)                                                                                                     aaaTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCaaa                                              <  4:414146/149‑1 (MQ=255)                                                                                                        tttGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTg                                           >  4:489627/1‑149 (MQ=255)                                                                                                        tttGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAACTg                                           >  1:259305/1‑149 (MQ=255)                                                                                                                 gCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTa                                  <  1:780582/149‑1 (MQ=255)                                                                                                                     aGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGttt                              <  3:26915/149‑1 (MQ=255)                                                                                                                     aGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGttt                              <  1:379173/149‑1 (MQ=255)                                                                                                                       cGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTCGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCa                            >  3:788790/1‑149 (MQ=255)                                                                                                                         ggATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAAt                          >  4:86212/1‑149 (MQ=255)                                                                                                                           aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCg                        >  3:573701/1‑149 (MQ=255)                                                                                                                            tCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGc                       >  1:362112/1‑149 (MQ=255)                                                                                                                             cGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCt                      <  3:51602/149‑1 (MQ=255)                                                                                                                              gAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTg                     >  1:418104/1‑149 (MQ=255)                                                                                                                               aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGa                    <  2:772469/149‑1 (MQ=255)                                                                                                                                     tGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATc              >  3:34436/1‑149 (MQ=255)                                                                                                                                     tGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTCGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATc              <  4:788790/149‑1 (MQ=255)                                                                                                                                       aCCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAg            <  3:6774/149‑1 (MQ=255)                                                                                                                                         ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGc          >  3:839402/1‑149 (MQ=255)                                                                                                                                          gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAg            <  2:114219/145‑1 (MQ=255)                                                                                                                                          gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAg            >  1:114219/2‑146 (MQ=255)                                                                                                                                          gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa         >  3:179581/2‑149 (MQ=255)                                                                                                                                               catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGa    <  2:382436/149‑1 (MQ=255)                                                                                                                                                atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGTTGAAAAATCAGGCAAATGaa   <  4:417010/149‑1 (MQ=255)                                                                                                                                                 tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc  >  1:681189/1‑149 (MQ=255)                                                                                                                                                 tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc  >  3:850020/1‑149 (MQ=255)                                                                                                                                                 tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc  >  1:194505/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAC  >  NC_000913/2518051‑2518343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑