Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,481,585	T→A	6.9%	M226K (ATG→AAG)	ynbD →	putative phosphatase inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,481,585	0	T	A	6.9%	211.0 / 4.1	69	M226K (ATG→AAG)	ynbD	putative phosphatase inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/36);  new base A (3/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCATCACAGGGCTGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTACCGCCCGCCG  >  NC_000913/1481442‑1481728                                                                                                                                                |                                                                                                                                                tGTCATCACAGGGCTGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGc                                                                                                                                            <  4:236789/149‑1 (MQ=255)    catcaCAGGGCTGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTaa                                                                                                                                         >  5:6835/1‑149 (MQ=255)     atcaCAGGGCTGGCGGTAGGTATGTTGATTGACTGTATGGTGAACGTCTACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGAACTGCCCTATGTCGTAGGAGCGGGCTCGTGAATTGTATTGATGGAGCTAat                                                                                                                                        >  1:201213/1‑149 (MQ=255)      tcaCAGGGCTGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAat                                                                                                                                        >  1:17/1‑148 (MQ=255)        acaGGGCTGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAatgat                                                                                                                                     >  3:354818/1‑149 (MQ=255)          aGGGCTGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAatgatga                                                                                                                                   <  8:147449/149‑1 (MQ=255)               tGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATGAGTAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGAtt                                                                                                                              <  7:270255/149‑1 (MQ=255)                 gCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGAtt                                                                                                                              <  2:251570/147‑1 (MQ=255)                    gtaggtaTGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGtt                                                                                                                         <  8:371455/149‑1 (MQ=255)                     taggtaTGTTGATTGACTGGATGGTGCCGGTCGTCCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATGGATGGAGCTAATGATGATGATTCAGTTa                                                                                                                        <  6:347324/149‑1 (MQ=255)                        gtaTGTTGATTGACTGGATGGTGCGCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTTCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTAtg                                                                                                                      <  1:278660/148‑1 (MQ=255)                         taTGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTAtggt                                                                                                                    <  6:6835/149‑1 (MQ=255)                           tGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTAtggtgg                                                                                                                  >  4:137964/1‑149 (MQ=255)                            gTTGATTGACTGTATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTTATGGAGCTAATGATGATGATTCAGTTATggtggt                                                                                                                 >  3:313153/1‑149 (MQ=255)                             ttgattgaCTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATggtggt                                                                                                                 <  4:348514/148‑1 (MQ=255)                             ttgattgaCTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTACATTGTATTGATGGAGCTAATGATGATGATTCAGTTATggtggt                                                                                                                 >  7:231058/1‑148 (MQ=255)                               gattgtCTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAACCCTGATCAACGTCGATTCTAAATAGCACTGCCCTTTGTCGTAGGCGCGGGCTCGTGTATTGTATTTATGGAGCTAATGATGATGATTGAGTTATggtggtga                                                                                                               <  6:280690/148‑3 (MQ=255)                                 ttgaCTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTc                                                                                                            <  4:354818/149‑1 (MQ=255)                                        gATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTaa                                                                                                                                         <  6:355416/113‑1 (MQ=255)                                        gATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTaa                                                                                                                                         >  5:355416/1‑113 (MQ=255)                                        gATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGtt                                                                                                       <  8:348830/147‑1 (MQ=255)                                        gATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTAt                                                                                                     >  1:339994/1‑149 (MQ=255)                                             tGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTgg                                                                                                >  5:28401/1‑149 (MQ=255)                                              gCCCGTCGACCGTCGTTGTAATTATCAGAAACATGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTAAGTTATGGTGGTCAGTCTGGTTATGTTTGc                                                                                               >  3:300231/1‑149 (MQ=255)                                                ccGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCa                                                                                             <  4:245427/149‑1 (MQ=255)                                                  gTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCa                                                                                             <  6:229776/147‑1 (MQ=255)                                                  gTCGACCGTAGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCa                                                                                             >  5:229776/1‑147 (MQ=255)                                                       ccGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGGGGGCAGTCTGGTTATGTTGGCCAGtattat                                                                                      <  4:273003/149‑1 (MQ=255)                                                          tcgtTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCTCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATggtggt                                                                                                                 >  5:342599/1‑119 (MQ=255)                                                             ttGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTAc                                                                                <  4:154940/149‑1 (MQ=255)                                                                 aaTTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGTAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCa                                                                             >  8:305532/1‑148 (MQ=255)                                                                  aTTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCAt                                                                            >  7:346872/1‑148 (MQ=255)                                                                      tCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGAAGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTTTTATCGCTACTCATAATTg                                                                       >  3:230690/1‑149 (MQ=255)                                                                        aGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGc                                                                     >  6:164332/1‑149 (MQ=255)                                                                        aGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGc                                                                     <  3:137964/149‑1 (MQ=255)                                                                           aaCCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTACATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGt                                                                  <  8:231058/149‑1 (MQ=255)                                                                             ccTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTgg                                                                <  2:17/149‑1 (MQ=255)                                                                             ccTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCACGGGCTCGTGCATTGTATTGATGGTGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTgg                                                                <  5:371676/149‑1 (MQ=255)                                                                             ccTGATCAACCTCGAATCATAATAGCACTGCCCTTTGTCGTAGGCGGTGGCTAGTGCTTTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTgg                                                                <  4:230690/149‑1 (MQ=255)                                                                              cTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTggg                                                               >  5:129995/1‑149 (MQ=255)                                                                                gATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTa                                                             >  7:70655/1‑149 (MQ=255)                                                                                    aaCGTCGAATCAAAATAGCACTGACCTATGTCGTAGGCGCGGTCTCGTGCATTGTATTGATTGAGCTAATGATGATTATTCAGTTATGGTGGTCAGTATGGTTATGTTGGCCAGTA‑TATCGCTACTCATAATTGGCCGTGTGTACGGt                                                         <  2:201213/148‑1 (MQ=255)                                                                                     aCGTCGAATCAAAATAGCTCTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTg                                                        <  8:369126/149‑1 (MQ=255)                                                                                             tCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGgc                                                >  8:302547/1‑149 (MQ=255)                                                                                               aaaaTAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGgcg                                               <  6:129995/148‑1 (MQ=255)                                                                                                     gCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaa                                         <  5:164332/148‑1 (MQ=255)                                                                                                     gCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGGCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaa                                         <  3:80821/148‑1 (MQ=255)                                                                                                       aCTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGAAGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaacaa                                      >  4:284173/1‑149 (MQ=255)                                                                                                       aCTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGAAGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaaca                                       >  4:289082/1‑148 (MQ=255)                                                                                                           cccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGtattat                                                                                      >  4:188999/1‑97 (MQ=255)                                                                                                           cccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGtattat                                                                                      <  3:188999/97‑1 (MQ=255)                                                                                                           cccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGTTTGGCGCGATAacaaca                                    <  5:369533/147‑1 (MQ=255)                                                                                                           cccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAacaaca                                    <  1:70121/147‑1 (MQ=255)                                                                                                           cccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAgg                                  >  1:135658/1‑149 (MQ=255)                                                                                                                tGTCGTAGGCGCGGGCTCGTGCATTGTATTGAAGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGaa                               <  3:284173/147‑1 (MQ=255)                                                                                                                tGTCGTAGGCGCGGGCTCGTGCATTGTATTGAAGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGaa                               <  3:289082/147‑1 (MQ=255)                                                                                                                  tCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGa                            >  3:285940/1‑148 (MQ=255)                                                                                                                  tCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGa                            <  2:339994/148‑1 (MQ=255)                                                                                                                     tAGTCGCGGGCTCGTGCAGTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGGTATGTTGGCCAGTATTATCGCTAATCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGt                        >  4:294140/1‑149 (MQ=255)                                                                                                                     tAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGt                        <  8:346872/149‑1 (MQ=255)                                                                                                                          gcgGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAggg                   >  3:319263/1‑149 (MQ=255)                                                                                                                            gggCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTGGGCGCGATAACAACAGGGAAAGAAAGTCAGGGGa                 <  8:70655/149‑1 (MQ=255)                                                                                                                              gCTCGTGCATTGTATTTATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGaaa               <  4:313153/149‑1 (MQ=255)                                                                                                                                  gTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTa            <  6:28401/148‑1 (MQ=255)                                                                                                                                    gCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGa                            <  4:19221/130‑1 (MQ=255)                                                                                                                                    gCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGa                            >  3:19221/1‑130 (MQ=255)                                                                                                                                        tGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGAACACTATTATCCCTACTCATAATTGGCCGTCGGTACGGTTCTCTAGGCGGGATAACAACAGCGAAAGAAAGTCACAGGAAACTAccgcccg     >  8:380888/1‑149 (MQ=255)                                                                                                                                         gTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTAccgcccgc    >  2:18484/1‑149 (MQ=255)                                                                                                                                           aTTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTACCGCccgccg  >  2:363109/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                TGTCATCACAGGGCTGGCGGTAGGTATGTTGATTGACTGGATGGTGCCCGTCGACCGTCGTTGGAATTATCAGAAACCTGATCAACGTCGAATCAAAATAGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTACCGCCCGCCG  >  NC_000913/1481442‑1481728

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