Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	486,849	C→T	5.9%	A105V (GCG→GTG)	mscK →	potassium dependent, small conductance mechanosensitive channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	486,849	0	C	T	5.9%	213.7 / 2.3	67	A105V (GCG→GTG)	mscK	potassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (37/26);  new base T (2/2);  total (39/28)
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.83e-01

ACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATA  >  NC_000913/486703‑486998                                                                                                                                                   |                                                                                                                                                      aCAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGAc                                                                                                                                                    <  2:709945/150‑1 (MQ=255)   aaaaaGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACcg                                                                                                                                                  >  1:388807/1‑150 (MQ=255)           cTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTa                                                                                                                                          <  1:1211935/150‑1 (MQ=255)           cTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTa                                                                                                                                          <  1:553495/150‑1 (MQ=255)                 gCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCa                                                                                                                                    >  2:748152/1‑150 (MQ=255)                      ggACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAg                                                                                                                               >  2:644707/1‑150 (MQ=255)                       gACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGGCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGc                                                                                                                              <  1:1127719/150‑1 (MQ=255)                         cAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGa                                                                                                                            <  2:986820/150‑1 (MQ=255)                           aaCTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATg                                                                                                                          >  2:892966/1‑150 (MQ=255)                             cTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTc                                                                                                                        <  2:1404344/150‑1 (MQ=255)                               ggTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGa                                                                                                                      <  1:1250623/150‑1 (MQ=255)                                 tGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATa                                                                                                                    <  1:512209/150‑1 (MQ=255)                                     gcagGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAacga                                                                                                                >  2:156262/1‑150 (MQ=255)                                     gcagGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAacga                                                                                                                >  1:202355/1‑150 (MQ=255)                                        ggATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAAcgacga                                                                                                             <  2:43075/150‑1 (MQ=255)                                            cTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACgaagaa                                                                                                         >  1:1028352/1‑150 (MQ=255)                                             tGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGaaa                                                                                                        >  2:1396362/1‑150 (MQ=255)                                             tGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGaaa                                                                                                        >  1:162084/1‑150 (MQ=255)                                                cAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAAcgc                                                                                                     <  1:756019/150‑1 (MQ=255)                                                cAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAAcgc                                                                                                     >  2:1163354/1‑150 (MQ=255)                                                cAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAAcgc                                                                                                     >  1:15433/1‑150 (MQ=255)                                                cAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAAcgc                                                                                                     >  1:679243/1‑150 (MQ=255)                                                   aTACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCa                                                                                                  <  2:1270180/150‑1 (MQ=255)                                                    tACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCaa                                                                                                 >  1:275608/1‑150 (MQ=255)                                                      cATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCaaaa                                                                                               >  1:283007/1‑150 (MQ=255)                                                      cATTAGCCACCCTCGATAAAATAGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCaaaa                                                                                               >  2:472733/1‑150 (MQ=255)                                                          aGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGTGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCt                                                                                           >  1:673083/1‑150 (MQ=255)                                                          aGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGTGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCt                                                                                           >  1:678103/1‑150 (MQ=255)                                                          aGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCt                                                                                           >  1:167969/1‑150 (MQ=255)                                                           gCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTg                                                                                          >  2:291103/1‑150 (MQ=255)                                                                ccTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCAc                                                                                     <  2:162084/150‑1 (MQ=255)                                                                 cTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACg                                                                                    <  2:46623/150‑1 (MQ=255)                                                                  tCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGc                                                                                   >  1:499199/1‑150 (MQ=255)                                                                      tAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTc                                                                               <  1:743248/150‑1 (MQ=255)                                                                         aaTCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGtt                                                                            >  1:69687/1‑150 (MQ=255)                                                                         aaTCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGtt                                                                            >  1:1409753/1‑150 (MQ=255)                                                                          atcgatcgCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTg                                                                           >  1:545322/1‑150 (MQ=255)                                                                              atcgCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGcc                                                                       <  2:1409753/150‑1 (MQ=255)                                                                                   aTAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTg                                                                  <  1:1396362/150‑1 (MQ=255)                                                                                    tAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTgg                                                                 >  1:1377021/1‑150 (MQ=255)                                                                                         gaagaGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACAACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACt                                                            >  1:280081/1‑150 (MQ=255)                                                                                           agagACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTcg                                                          >  1:431128/1‑150 (MQ=255)                                                                                           agagACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTAAGCACGCTGTCGTTGCGCCAGCTGGAAACTcg                                                          >  2:991022/1‑150 (MQ=255)                                                                                               aCAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGtt                                                      <  2:283007/150‑1 (MQ=255)                                                                                                 aGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGc                                                    >  1:492469/1‑150 (MQ=255)                                                                                                  gTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGcc                                                   >  1:1379106/1‑150 (MQ=255)                                                                                                  gTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGcc                                                   >  2:240339/1‑150 (MQ=255)                                                                                                   ttCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGccc                                                  <  2:21142/150‑1 (MQ=255)                                                                                                    tCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCa                                                 <  1:748152/150‑1 (MQ=255)                                                                                                     cAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCTTTGCCCAg                                                <  2:679243/150‑1 (MQ=255)                                                                                                       gCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGgc                                              >  1:348938/1‑150 (MQ=255)                                                                                                         tACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGc                                                                                   <  1:970187/111‑1 (MQ=255)                                                                                                         tACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGc                                                                                   >  2:970187/1‑111 (MQ=255)                                                                                                              cAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGAc                                       <  1:391529/150‑1 (MQ=255)                                                                                                              cAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGAc                                       <  2:1148815/150‑1 (MQ=255)                                                                                                                   aGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGAttt                                  <  2:69687/150‑1 (MQ=255)                                                                                                                      cGCTGAAGCGCCGGAAAAAATGCGCCAGGTGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCa                               <  2:673083/150‑1 (MQ=255)                                                                                                                      cGCTGAAGCGCCGGAAAAAATGCGCCAGGTGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCa                               <  2:678103/150‑1 (MQ=255)                                                                                                                      cGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCa                               >  2:1269716/1‑150 (MQ=255)                                                                                                                         tGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCaaaa                            >  1:1036910/1‑150 (MQ=255)                                                                                                                               gcCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGcaca                      <  1:892966/150‑1 (MQ=255)                                                                                                                                ccGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCATAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACaa                     >  1:521937/1‑150 (MQ=255)                                                                                                                                 cGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACaaa                    >  2:1163924/1‑150 (MQ=255)                                                                                                                                        aaTGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCt             >  1:235895/1‑150 (MQ=255)                                                                                                                                            cgcCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCg         >  2:1314589/1‑150 (MQ=255)                                                                                                                                             gcCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCgg        >  2:454742/1‑149 (MQ=255)                                                                                                                                                 ggCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTta    <  1:934005/150‑1 (MQ=255)                                                                                                                                                   cGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTtata  <  2:431128/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATA  >  NC_000913/486703‑486998

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