Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,055,120	G→C	5.2%	L20L (CTG→CTC)	amn →	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,055,120	0	G	C	5.2%	282.0 / 3.2	85	L20L (CTG→CTC)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base G (38/42);  new base C (2/2);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCC  >  NC_000913/2054973‑2055265                                                                                                                                                    |                                                                                                                                                  atcatcCGATCTTCTCCTTACTCTTTTGTTTTAGATTTCTGTTTCTTTCTATCGTTTGTGCGGTTATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTGAGGCACTGGTTATTCTCGACGGGCTGt                                                                                                                                                  <  8:62146/149‑1 (MQ=255)      ccGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAg                                                                                                                                             >  4:203020/1‑149 (MQ=255)         aTCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCaa                                                                                                                                          >  4:362470/1‑149 (MQ=255)              cGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTAt                                                                                                                                                <  6:435200/138‑1 (MQ=255)              cGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTg                                                                                                                                      <  2:237074/148‑1 (MQ=255)                  ttACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTc                                                                                                                                 >  6:163048/1‑149 (MQ=255)                      acTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCAt                                                                                                                             <  2:264912/149‑1 (MQ=255)                          ttGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTAc                                                                                                                          <  7:464933/148‑1 (MQ=255)                             tttCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAAGCGCTGTATGAGCAATCTGTAGTCGCATTACGCaa                                                                                                                      >  3:369445/1‑149 (MQ=255)                              ttCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAAc                                                                                                                     >  6:131055/1‑149 (MQ=255)                              ttCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAAc                                                                                                                     >  8:289717/1‑149 (MQ=255)                               taacaTTTCTGTGACATACTATCGGATATGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAACTCAGGCACAGGATAAACTCAACGCGCTATATGATCAAACTGTAATAGCATTACGCaaag                                                                                                                    >  4:48299/3‑147 (MQ=255)                               tCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACg                                                                                                                    >  7:274321/1‑149 (MQ=255)                                cacaTTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGATCTATGAGCAATCt                                                                                                                                       <  6:47112/129‑1 (MQ=255)                                cacaTTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAACCCAGCTCATGCACTGGATAAACTCGACGCGCTCTATGAGCAATca                                                                                                                                       >  5:53753/1‑128 (MQ=255)                                cacaTTTCTGTGACATACTATCGGATGTGCGGGAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTCTATGAGCAATCt                                                                                                                                       >  5:47112/1‑129 (MQ=255)                                caactttATGTGACATACTATCGTATGTGCGGAAATTGTATGTAACATGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTCTATGAGCAATCt                                                                                                                                       <  6:53753/125‑1 (MQ=255)                                    tttCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCAtt                                                                                                               <  3:203020/149‑1 (MQ=255)                                        tgtgACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGc                                                                                                            >  4:226855/1‑148 (MQ=255)                                        tgtgACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGc                                                                                                            >  4:335456/1‑148 (MQ=255)                                          tgACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCaa                                                                                                          <  5:131055/148‑1 (MQ=255)                                          tgACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAAc                                                                                                         <  7:289717/149‑1 (MQ=255)                                            aCATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACta                                                                                                       >  5:132189/1‑149 (MQ=255)                                            aCATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACta                                                                                                       <  5:163048/149‑1 (MQ=255)                                              aTACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACtata                                                                                                     >  3:202324/1‑149 (MQ=255)                                                 cTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATAtt                                                                                                   <  7:237416/148‑1 (MQ=255)                                                  tATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTa                                                                                                  >  3:322091/1‑148 (MQ=255)                                                   aTCGGATGTGAGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTAc                                                                                                 >  8:110153/1‑148 (MQ=255)                                                    tCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTAc                                                                                                 >  7:312502/1‑147 (MQ=255)                                                        aTGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTgg                                                                                           <  4:202324/149‑1 (MQ=255)                                                         tgtgCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGc                                                                                          <  3:34765/149‑1 (MQ=255)                                                         ggtgCCGTAATTGTATGCAACAGGAGACACACATGAATAATATGGGCTCCGGTCTGACCCCAGCTCATGCACTGGCTAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGc                                                                                          <  8:455956/148‑1 (MQ=255)                                                            gCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTa                                                                                                                                    >  6:450003/1‑104 (MQ=255)                                                            gCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTa                                                                                                                                    <  5:450003/104‑1 (MQ=255)                                                             cGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAAt                                                                                      <  3:48299/149‑1 (MQ=255)                                                                  aTTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTAcc                                                                                  <  3:362470/148‑1 (MQ=255)                                                                     gTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCtga                                                                               <  8:274321/148‑1 (MQ=255)                                                                        tGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCATCTCAGGCACTGGATAAACTCTACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACTCAATTGGCAACTATATTACAAGTGTCGAATTACAtgatga                                                                            >  5:435709/1‑148 (MQ=255)                                                                        tGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTg                                                                                                              <  5:341703/114‑1 (MQ=255)                                                                        tGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTg                                                                                                              >  6:341703/1‑114 (MQ=255)                                                                        tGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCtgatga                                                                            >  6:137667/1‑148 (MQ=255)                                                                         ggAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGaa                                                                           <  8:282635/148‑1 (MQ=255)                                                                               ggAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGccc                                                                    <  6:132189/149‑1 (MQ=255)                                                                                  gaCACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGc                                                                  >  8:271987/1‑148 (MQ=255)                                                                                    cacacaTGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCaaa                                                               <  4:322091/149‑1 (MQ=255)                                                                                     acacaTGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAAc                                                              <  2:215211/149‑1 (MQ=255)                                                                                     acacaTGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAAc                                                              <  5:347199/149‑1 (MQ=255)                                                                                       acaTGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCTCTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACaa                                                            <  2:7450/149‑1 (MQ=255)                                                                                       acaTGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAAc                                                              >  5:426082/1‑147 (MQ=255)                                                                                       acaTGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACaa                                                            >  1:7450/1‑149 (MQ=255)                                                                                              taataaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCt                                                                                 >  8:451177/1‑121 (MQ=255)                                                                                              taataaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCt                                                                                 <  7:451177/121‑1 (MQ=255)                                                                                              taataaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCttt                                                     <  6:426082/149‑1 (MQ=255)                                                                                                ataaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCTATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGTAAACGCCCGCAAACAAGGTCtttt                                                    <  7:110153/148‑1 (MQ=255)                                                                                                ataaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCttttt                                                   <  3:335456/149‑1 (MQ=255)                                                                                                 taaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATGTGTAGTCGCATTGCGCAACGCCATTGGCAACTTTATTGCAAGTGGCGAATTACCTGATGAAAACGCCCGCAAGCCAGGtggtgttg                                                  >  4:173705/1‑142 (MQ=255)                                                                                                 taaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCttt                                                     <  3:210388/146‑1 (MQ=255)                                                                                                 taaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCttt                                                     >  4:210388/1‑146 (MQ=255)                                                                                                 taaGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTg                                                  >  5:71854/1‑149 (MQ=255)                                                                                                   aGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGt                                                 <  6:71854/148‑1 (MQ=255)                                                                                                    gggCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCt                                               <  5:137667/149‑1 (MQ=255)                                                                                                        tCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATTAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATc                                            >  8:436216/1‑148 (MQ=255)                                                                                                        tCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATc                                            >  8:8849/1‑148 (MQ=255)                                                                                                            gTCTGACCCCAGCGCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCa                                       <  6:435709/149‑1 (MQ=255)                                                                                                                 aCCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGAc                                  <  3:226855/149‑1 (MQ=255)                                                                                                                   cccAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCg                                >  5:434555/1‑149 (MQ=255)                                                                                                                   cccAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCg                                >  7:459725/1‑149 (MQ=255)                                                                                                                   cccAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCg                                >  7:186168/1‑149 (MQ=255)                                                                                                                     cAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGTGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGt                               >  8:158295/1‑148 (MQ=255)                                                                                                                        cTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAAcc                           >  3:196113/1‑149 (MQ=255)                                                                                                                           aGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCt                          <  6:434555/147‑1 (MQ=255)                                                                                                                               aCTGGATAAACTCAACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACg                    >  4:457688/1‑149 (MQ=255)                                                                                                                                   gATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTa                 <  4:196113/148‑1 (MQ=255)                                                                                                                                   gATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTa                 <  7:271987/148‑1 (MQ=255)                                                                                                                                    aTAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAg                >  7:2049/1‑148 (MQ=255)                                                                                                                                        aCTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGAc                                  <  7:350728/126‑1 (MQ=255)                                                                                                                                        aCTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGAc                                  >  8:350728/1‑126 (MQ=255)                                                                                                                                        aCTCGACGCGCTGTATGAGCAATCTGCAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACg                    <  3:117436/140‑1 (MQ=255)                                                                                                                                        aCTCGACGCGCTGTATGAGCAATCTGAAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACg                    >  4:117436/1‑140 (MQ=255)                                                                                                                                          tCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAAc          >  3:384818/1‑148 (MQ=255)                                                                                                                                           cGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGt                  >  5:189743/1‑139 (MQ=255)                                                                                                                                           cGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGt                  <  6:189743/139‑1 (MQ=255)                                                                                                                                            gACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCaa       <  7:436216/149‑1 (MQ=255)                                                                                                                                               gcgcTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATcc    <  8:312502/149‑1 (MQ=255)                                                                                                                                                 gcTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACGTGGGACGGTAGCACAACCAATcccc  <  8:186168/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACCAATCCCC  >  NC_000913/2054973‑2055265

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