Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,372,447	A→G	5.2%	L43P (CTC→CCC)	nanT ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,372,447	0	A	G	5.2%	242.7 / 2.8	76	L43P (CTC→CCC)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base G (2/2);  total (39/37)
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.64e-01

TGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGG  >  NC_000913/3372301‑3372593                                                                                                                                                   |                                                                                                                                                   tGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGt                                                                                                                                                  <  3:344261/149‑1 (MQ=255)    tGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGGGt                                                                                                                                                  <  6:167139/146‑1 (MQ=255)    tGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGGGt                                                                                                                                                  >  5:167139/1‑146 (MQ=255)    tGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGGGt                                                                                                                                                  <  4:216895/146‑1 (MQ=255)    tGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGGGt                                                                                                                                                  >  3:216895/1‑146 (MQ=255)         aCCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGGCAGCCCGAAATCACCTTGTACTTCGGTGAGTACCAggg                                                                                                                                           <  8:394178/148‑1 (MQ=255)          ccATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCCCCGTCGTCCGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGc                                                                                                                                          >  6:95710/1‑148 (MQ=255)           cATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGa                                                                                                                                        >  1:56312/1‑149 (MQ=255)           cATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCATCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGa                                                                                                                                        >  5:49261/1‑149 (MQ=255)                 cAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTaaaa                                                                                                                                  <  2:56312/149‑1 (MQ=255)                 cAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCAGCCGAACCAGAGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTaaaa                                                                                                                                  >  6:138003/1‑149 (MQ=255)                  aGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAAc                                                                                                                                 <  6:58628/149‑1 (MQ=255)                   gacgacGGCCCTAGCGGTCACCCATAGCGCCGAGCATCAGGCAGCCGAACCAGCTAGAGATAAAGGCTGCAGAGAGCAGACTTGCCGCCTGCACCGTCGGCAGCCCGAATTAACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACg                                                                                                                                <  1:356171/149‑1 (MQ=255)                   gacgacGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAAc                                                                                                                                 >  1:54184/1‑148 (MQ=255)                   gacgacGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGAGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGGCTGCACCGTCGTCAGCCCTATTTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACg                                                                                                                                <  1:350683/149‑1 (MQ=255)                      gacgCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGaaa                                                                                                                             >  8:284422/1‑149 (MQ=255)                      gacgCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGaaa                                                                                                                             <  7:284422/149‑1 (MQ=255)                      gacgCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGaaa                                                                                                                             <  3:45606/149‑1 (MQ=255)                            cGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCaaaa                                                                                                                       <  5:237888/149‑1 (MQ=255)                                 cGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCATCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTc                                                                                                                  <  6:49261/149‑1 (MQ=255)                                    tCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAg                                                                                                               >  7:230868/1‑149 (MQ=255)                                      aCCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCa                                                                                                             <  2:97588/149‑1 (MQ=255)                                        ccATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAg                                                                                                               >  8:1044/1‑145 (MQ=255)                                        ccATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCat                                                                                                            >  3:172371/1‑147 (MQ=255)                                         cATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGa                                                                                                           >  5:264724/1‑148 (MQ=255)                                             gcgcCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATc                                                                                                       <  8:339190/148‑1 (MQ=255)                                             gcgcCGAGCATCAGGACGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCCTCGTCAGCCCGAATTCACCTTGAACTTAGGTGAGTACCAGGGCGATTAAAACGAAATCAAACCCGTCAAGCAGATATcc                                                                                                      >  7:316604/1‑149 (MQ=255)                                               gcCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGCTCAGACTTGCCGCCTGCACCGTCGTCCGCCCGAATTCCCCTTGTACTTCGGTGAGTCCCAGGGCGATAAAAACGAAATCAACACCGTCAAGCAGATATacc                                                                                                     >  2:235802/1‑148 (MQ=255)                                                ccGAGCATAAGTCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGAGCAGACTTGCCGCCTCCACAGTCGTCAGCCCGAATTCACCTTGTACTTCGTTGAGTACAAGGGCGAGTAAAACGAAATCAAAACCGTCAAGCAGATTTCCCaa                                                                                                   >  2:253180/1‑149 (MQ=255)                                                 cGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAAc                                                                                                  <  4:337044/149‑1 (MQ=255)                                                        cAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAg                                                                                           >  7:243134/1‑149 (MQ=255)                                                         aGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGc                                                                                          <  6:264724/149‑1 (MQ=255)                                                                  aaCCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATgc                                                                                 <  5:126192/149‑1 (MQ=255)                                                                    ccAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATgcgc                                                                               <  5:138003/149‑1 (MQ=255)                                                                              aTAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATc                                                                                                       >  8:342824/1‑115 (MQ=255)                                                                              aTAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATc                                                                                                       <  7:342824/115‑1 (MQ=255)                                                                               tAAAGGCTGCAGAGATCAGACTTGCCGCCTGCCCCGTCGTCAGCCCGAATTCCCCTTGTACTTCGGTGAGTACCAGGGCGCTTAAAACGAACTCAAAACCGTCAAGCAGATATCCCACCCAGGCAGCGGAAAATCCGCGCCATTCTGCa                                                                    >  7:314069/1‑149 (MQ=255)                                                                                aaaGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCAc                                                                   >  5:203306/1‑149 (MQ=255)                                                                                     cTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGtt                                                               >  8:387323/1‑148 (MQ=255)                                                                                     cTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGAACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGtt                                                               >  7:355054/1‑148 (MQ=255)                                                                                      tGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                             <  6:203306/149‑1 (MQ=255)                                                                                          gagaGCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGGGAGTACCAGGGCGATAAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATg                                                         <  1:235802/149‑1 (MQ=255)                                                                                           agaTCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTgaga                                                           <  6:322754/146‑1 (MQ=255)                                                                                           agaTCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTgaga                                                           >  5:322754/1‑146 (MQ=255)                                                                                            gaTCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGc                                                       >  7:330312/1‑149 (MQ=255)                                                                                             aTCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGc                                                       >  1:287412/1‑148 (MQ=255)                                                                                             aTCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGc                                                       >  5:331192/1‑148 (MQ=255)                                                                                              tCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGa                                                     <  8:243134/149‑1 (MQ=255)                                                                                                aGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                             >  2:390200/1‑139 (MQ=255)                                                                                                aGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                             <  1:390200/139‑1 (MQ=255)                                                                                                  aCTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTg                                                              >  7:273663/1‑136 (MQ=255)                                                                                                  aCTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTg                                                              <  8:273663/136‑1 (MQ=255)                                                                                                       ccgccTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCTAAACCGTCAAGCATATATCCCAACCAGGCTGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGTTACCACgg                                             <  8:211931/148‑1 (MQ=255)                                                                                                               aCCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTg                                    <  2:54184/149‑1 (MQ=255)                                                                                                                ccGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATtgt                                                                       <  6:216833/113‑1 (MQ=255)                                                                                                                ccGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATtgt                                                                       >  5:216833/1‑113 (MQ=255)                                                                                                                 cgtcgtCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTggg                                  <  5:84764/149‑1 (MQ=255)                                                                                                                  gtcgtcAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGt                                 >  1:34371/1‑149 (MQ=255)                                                                                                                    cgtcAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTg                               <  5:95710/149‑1 (MQ=255)                                                                                                                      tcAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTgta                             <  8:230868/149‑1 (MQ=255)                                                                                                                      tcAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACAGTCAAGCAGATATACCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTg                               >  3:278139/1‑147 (MQ=255)                                                                                                                          cccGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGt                                                                                                                   >  8:113463/1‑59 (MQ=255)                                                                                                                          cccGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTAc                         >  7:187327/1‑149 (MQ=255)                                                                                                                            cGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                        >  6:361133/1‑148 (MQ=255)                                                                                                                                 tCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATtgt                  >  8:298721/1‑149 (MQ=255)                                                                                                                                 tCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATtgt                  >  8:201566/1‑149 (MQ=255)                                                                                                                                   cgcTTGTACTTCGGCGAGTACCAGGGCGATTAAAACGAATTCAATACCGTCAAGCAGATATCCCTACCAGGCAGCGTATAATGCGCGCCATTGTGCACGGTTGTGATGGGGATACCAAGGGATATTCTGTGTTGTAGTACTCATtgtgt                <  8:316604/147‑2 (MQ=255)                                                                                                                                   aCCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATtgt                  >  2:223967/1‑147 (MQ=255)                                                                                                                                   aCCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGa                <  7:387323/149‑1 (MQ=255)                                                                                                                                     cTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTCTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGt              >  3:165946/1‑149 (MQ=255)                                                                                                                                     cTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTCTGCACGCTTGAGATGCCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGAGAGt              >  3:171486/1‑149 (MQ=255)                                                                                                                                        gTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGtct            <  6:331192/148‑1 (MQ=255)                                                                                                                                            ttCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCcgc       >  7:305306/1‑149 (MQ=255)                                                                                                                                                gTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTg   <  7:201566/149‑1 (MQ=255)                                                                                                                                                 tGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTgg  <  8:330312/149‑1 (MQ=255)                                                                                                                                                 tGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTgg  <  5:361133/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCAGAGATCAGACTTGCCGCCTGCACCGTCGTCAGCCCGAATTCACCTTGTACTTCGGTGAGTACCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGG  >  NC_000913/3372301‑3372593

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