Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,395,094	A→C	5.5%	Y277D (TAT→GAT)	ynaI ←	mechanosensitive channel protein, very small conductance

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,395,094	0	A	C	5.5%	198.4 / 2.7	67	Y277D (TAT→GAT)	ynaI	mechanosensitive channel protein, very small conductance
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base C (2/2);  total (34/33)
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.90e-01

GCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTGGGT  >  NC_000913/1394949‑1395240                                                                                                                                                  |                                                                                                                                                   gCAAAGTCCGTGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAAc                                                                                                                                                 <  5:184754/149‑1 (MQ=255) gCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAAc                                                                                                                                                 >  3:300870/1‑149 (MQ=255)  cAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAAcc                                                                                                                                                <  6:42437/149‑1 (MQ=255)   aaaGTACGCGCCGTGTGAATTTTCAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCa                                                                                                                                               <  1:241366/149‑1 (MQ=255)      gTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGt                                                                                                                                             >  6:14845/1‑148 (MQ=255)         cgcgCCGTGTGACTGTACAATATCGATAATCTTCAAATAACCGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGGTTTGGTAAAGCAATAAACAATAAGATTCAATGAAGAGGACGCAAAATGGGTGAAATAAACCCGTACgg                                                                                                                                         >  5:40090/1‑149 (MQ=255)          gcgcCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGt                                                                                                                                        <  8:147930/149‑1 (MQ=255)          gcgcCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGt                                                                                                                                        >  7:343672/1‑149 (MQ=255)           cgcCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAACCTGGTTGAAATAAACCAGTAAGGtt                                                                                                                                       <  5:217078/149‑1 (MQ=255)              cGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGt                                                                                                                                    <  7:344673/149‑1 (MQ=255)                       gTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTc                                                                                                                            >  3:148070/1‑148 (MQ=255)                          cAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCt                                                                                                                                >  6:299667/1‑141 (MQ=255)                          cAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCt                                                                                                                                <  5:299667/141‑1 (MQ=255)                            atatCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGc                                                                                                                      >  7:153805/1‑149 (MQ=255)                             tatCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGcc                                                                                                                     <  5:327646/149‑1 (MQ=255)                                    aaTCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATaaa                                                                                                                                                  <  1:336840/113‑1 (MQ=255)                                    aaTCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATaaa                                                                                                                                                  >  2:336840/1‑113 (MQ=255)                                        ttCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTATAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGAttttt                                                                                                          >  6:297422/1‑149 (MQ=255)                                         tCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGAttttt                                                                                                          >  6:107794/1‑148 (MQ=255)                                          cAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTc                                                                                                         <  8:95840/148‑1 (MQ=255)                                           aaaTAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCa                                                                                                        <  8:343672/148‑1 (MQ=255)                                                      ttGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTa                                                                                            >  2:22073/1‑149 (MQ=255)                                                        gCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAAACATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTa                                                                                            <  6:50118/147‑1 (MQ=255)                                                        gCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTGTTGGTAAAGCAATACACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGGCGATGGCCGGGTGATTTTTCAGCATCTCACGTa                                                                                            <  4:269383/147‑1 (MQ=255)                                                         cTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACa                                                                                          <  5:14845/148‑1 (MQ=255)                                                          tgtgCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGc                                                                                        >  8:208137/1‑149 (MQ=255)                                                            tgCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCtt                                                                                      <  1:137226/149‑1 (MQ=255)                                                              cagcaAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTc                                                                                     <  7:15168/148‑1 (MQ=255)                                                                 caAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTAAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGAc                                                                                  >  8:130935/1‑148 (MQ=255)                                                                   aGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAAt                                                                               <  5:152004/149‑1 (MQ=255)                                                                        ttCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAAt                                                                               <  2:27838/144‑1 (MQ=255)                                                                        ttCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAAt                                                                               >  1:27838/1‑144 (MQ=255)                                                                               cATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACttt                                                                   <  7:208137/149‑1 (MQ=255)                                                                                 tACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTg                                                                 >  8:187643/1‑149 (MQ=255)                                                                                 tACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAAACAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTg                                                                 >  8:323549/1‑149 (MQ=255)                                                                                   ccGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGcc                                                               <  3:315973/149‑1 (MQ=255)                                                                                   ccGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGcc                                                               <  8:225253/149‑1 (MQ=255)                                                                                        gTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCAt                                                           <  1:141860/148‑1 (MQ=255)                                                                                         ttttGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAg                                                                                                       >  6:182250/1‑103 (MQ=255)                                                                                         ttttGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAg                                                                                                       <  5:182250/103‑1 (MQ=255)                                                                                         ttttGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGAGTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATcc                                                         <  4:311786/149‑1 (MQ=255)                                                                                         ttttGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGGGATTGTTCAGCATCTGAGGTGCAGCTTCGACAATAACGGCCACATTTGCCGCATcc                                                         >  6:278330/1‑149 (MQ=255)                                                                                                          aCCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACCGCTTCGACAATAACGCCCACTTTTGCCGCCTCCTCCTAACGTAAACCAAt                                        >  2:220589/1‑149 (MQ=255)                                                                                                           ccATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATg                                       >  5:347451/1‑149 (MQ=255)                                                                                                             ataataTTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGt                                     >  8:70309/1‑149 (MQ=255)                                                                                                              taataTTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTc                                    >  8:208080/1‑149 (MQ=255)                                                                                                                atatTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGt                                  <  5:297422/149‑1 (MQ=255)                                                                                                                   ttCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATATTTCAg                                                                                                       >  3:120171/1‑77 (MQ=255)                                                                                                                      aaTGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAAt                            >  6:205430/1‑149 (MQ=255)                                                                                                                      aaTGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAACATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTa                              >  2:339895/1‑147 (MQ=255)                                                                                                                       aTGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGACCACTTTTGCCACATACTCATAACGTAAACCAATGGTCGTCGTAATg                           >  1:60199/1‑149 (MQ=255)                                                                                                                          cagtggaaGCAAACTGGTTGAACTATAGGAGTTTGGTTTGTCGCTGGTCGATGGCGGGGTGATTTTTCAGCATCTCACGTACATCTTAGACAATTACGCCCACTTTTGCCGCATCATCATAACGTTAACCAATGGTCTTGGTAATgcgg                        <  6:40090/142‑1 (MQ=255)                                                                                                                          aaGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATgcgg                        >  8:111348/1‑149 (MQ=255)                                                                                                                                cAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTg                  <  1:339895/149‑1 (MQ=255)                                                                                                                                  ccttaCTGGTTTATTTCAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGt                >  8:340989/5‑149 (MQ=255)                                                                                                                                   cAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGt                >  8:93592/1‑148 (MQ=255)                                                                                                                                   cAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGt                >  8:98918/1‑148 (MQ=255)                                                                                                                                    ttaCTGGTTTATTTCAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGAGTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTg                  <  5:257130/143‑1 (MQ=255)                                                                                                                                    ttaCTGGTTTATATCAACCAGTAAGGTTTGTCGCTGCTCGATGGCCGGGTGATTTTTCAGCATCTCACGAACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATTGACGTGGTAATGCGGCGGtagca                >  6:257130/3‑143 (MQ=255)                                                                                                                                     aaCTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCAt             <  5:205430/149‑1 (MQ=255)                                                                                                                                      aCTGGTTTATTTCAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCa              <  7:340989/147‑1 (MQ=255)                                                                                                                                      aCTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCAtt            >  5:352228/1‑149 (MQ=255)                                                                                                                                             gAAATAAACCTGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGCCAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTg     >  2:103615/1‑149 (MQ=255)                                                                                                                                             gAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCt      <  4:148070/148‑1 (MQ=255)                                                                                                                                              aaaTAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTgg    <  6:352228/149‑1 (MQ=255)                                                                                                                                              aaaTAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGGGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTgg    <  5:107794/149‑1 (MQ=255)                                                                                                                                               aaTAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTggg   >  6:125228/1‑149 (MQ=255)                                                                                                                                               aaTAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTggg   >  1:158511/1‑149 (MQ=255)                                                                                                                                                aTAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTGGGt  <  7:130935/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTCGACAATAACGCCCACTTTTGCCGCATCCTCATAACGTAAACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTGGGT  >  NC_000913/1394949‑1395240

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: ‑