Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,920,950	C→T	*80* (TAG→TAA)	atpE ←	F0 sector of membrane‑bound ATP synthase, subunit c

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,920,950	0	C	T	100.0%	296.3 / NA	91	*80* (TAG→TAA)	atpE	F0 sector of membrane‑bound ATP synthase, subunit c
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (49/42);  total (49/42)

CATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTC  >  NC_000913/3920823‑3921081                                                                                                                                |                                                                                                                                    cATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGc                                                                                                                          >  2:1251599/1‑139 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                        >  1:349865/1‑139 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                        >  2:374996/1‑139 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                        >  2:566411/1‑139 (MQ=255)     cagaacagaaCGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                      <  1:566411/139‑1 (MQ=255)        aacagaacGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                   >  2:550272/1‑139 (MQ=255)                 aacAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcaga                                                                                                          <  2:349865/139‑1 (MQ=255)                                aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                           >  2:889783/1‑139 (MQ=255)                                aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                           >  2:1366752/1‑139 (MQ=255)                                             ggATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAg                                                                              >  1:1034159/1‑139 (MQ=255)                                                 tgttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATc                                                                          <  2:1034159/139‑1 (MQ=255)                                                  gttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAg                                                                                             >  2:599965/1‑119 (MQ=255)                                                  gttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAg                                                                                             <  1:599965/119‑1 (MQ=255)                                                    tgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                       <  1:374996/139‑1 (MQ=255)                                                      cGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 >  1:840964/1‑127 (MQ=255)                                                      cGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 <  2:840964/127‑1 (MQ=255)                                                          aaGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGAcc                                                                 <  2:20165/139‑1 (MQ=255)                                                            gATTCTCAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCa                                                               <  1:889783/139‑1 (MQ=255)                                                               tCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                            <  1:1251599/139‑1 (MQ=255)                                                                    gCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATa                                                       >  1:1097920/1‑139 (MQ=255)                                                                     cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATaa                                                      >  1:193753/1‑139 (MQ=255)                                                                     cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATaa                                                      >  1:1177835/1‑139 (MQ=255)                                                                      acaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATaaa                                                     <  2:193753/139‑1 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACTCTTATTACGCGACAGCGa                                                                                                                        >  1:1212695/1‑68 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                        <  2:1372498/68‑1 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                        >  1:1372498/1‑68 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                        <  2:1212695/68‑1 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                 >  2:1030058/1‑139 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                 >  2:1383342/1‑139 (MQ=255)                                                                               ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                   <  1:266654/68‑1 (MQ=255)                                                                               ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                   >  2:266654/1‑68 (MQ=255)                                                                               ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGt                                            >  2:650182/1‑139 (MQ=255)                                                                                tATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTa                                           <  1:650182/139‑1 (MQ=255)                                                                                  tttAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACg                                         <  1:1383342/139‑1 (MQ=255)                                                                                  tttAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACg                                         <  2:1177835/139‑1 (MQ=255)                                                                                    tAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 <  1:881997/97‑1 (MQ=255)                                                                                    tAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 >  2:881997/1‑97 (MQ=255)                                                                                     aGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccag                                                                                                     >  2:22026/1‑76 (MQ=255)                                                                                     aGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccag                                                                                                     <  1:22026/76‑1 (MQ=255)                                                                                      gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACa                                                                                                                            <  1:654735/52‑1 (MQ=255)                                                                                      gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACa                                                                                                                            >  2:654735/1‑52 (MQ=255)                                                                                      gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      <  2:1200419/106‑1 (MQ=255)                                                                                      gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      >  1:1200419/1‑106 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGc                                                                                                                          <  2:990384/53‑1 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGc                                                                                                                          >  1:990384/1‑53 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      <  1:1063784/105‑1 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      >  2:1063784/1‑105 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagca                                    >  2:1009993/1‑139 (MQ=255)                                                                                         aaCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                                  <  2:294154/139‑1 (MQ=255)                                                                                         aaCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                                  <  1:1366752/139‑1 (MQ=255)                                                                                          aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      <  1:664183/102‑1 (MQ=255)                                                                                          aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTa                                           >  2:89110/1‑129 (MQ=255)                                                                                          aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTa                                           <  1:89110/129‑1 (MQ=255)                                                                                          aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCagag                                 >  1:749171/1‑139 (MQ=255)                                                                                          aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAAAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      >  2:664183/1‑102 (MQ=255)                                                                                            gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGa                               >  2:455165/1‑139 (MQ=255)                                                                                              tCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCa                                                                                                                  >  1:106894/1‑54 (MQ=255)                                                                                              tCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCa                                                                                                                  <  2:106894/54‑1 (MQ=255)                                                                                                tGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 >  2:1341740/1‑85 (MQ=255)                                                                                                tGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 <  1:1341740/85‑1 (MQ=255)                                                                                                tGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatca                           >  1:218739/1‑139 (MQ=255)                                                                                                 gATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcag                          >  2:738313/1‑139 (MQ=255)                                                                                                      tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGcagcag                                   <  1:187151/125‑1 (MQ=255)                                                                                                      tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGcagcag                                   >  2:187151/1‑125 (MQ=255)                                                                                                       gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGAACAACGGGATAGAATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGAACAgg                    >  2:1383756/1‑139 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              <  2:816642/48‑1 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              >  1:816642/1‑48 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                        <  2:799381/54‑1 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                        >  1:799381/1‑54 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt                   >  1:128842/1‑139 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt                   >  2:46793/1‑139 (MQ=255)                                                                                                          ctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTg                 <  1:1009993/139‑1 (MQ=255)                                                                                                            ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACt                                                <  1:256840/106‑1 (MQ=255)                                                                                                            ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACt                                                >  2:256840/1‑106 (MQ=255)                                                                                                             taaataaaAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                      >  2:466692/1‑35 (MQ=255)                                                                                                             taaataaaAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                      <  1:466692/35‑1 (MQ=255)                                                                                                             taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacca                                                                     >  2:1153384/1‑84 (MQ=255)                                                                                                             taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacca                                                                     <  1:1153384/84‑1 (MQ=255)                                                                                                                  aaaaGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                        <  1:217446/44‑1 (MQ=255)                                                                                                                  aaaaGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                        >  1:1350793/1‑44 (MQ=255)                                                                                                                  aaaaGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                        >  2:217446/1‑44 (MQ=255)                                                                                                                  aaaaGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                        <  2:1350793/44‑1 (MQ=255)                                                                                                                    aaGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              <  2:688783/36‑1 (MQ=255)                                                                                                                    aaGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              >  1:688783/1‑36 (MQ=255)                                                                                                                    aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg       <  1:707409/139‑1 (MQ=255)                                                                                                                     aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                   >  2:520918/1‑46 (MQ=255)                                                                                                                     aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                   <  1:520918/46‑1 (MQ=255)                                                                                                                     aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcgc      >  2:454886/1‑139 (MQ=255)                                                                                                                         aCGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                   <  2:423194/42‑1 (MQ=255)                                                                                                                         aCGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                   >  1:423194/1‑42 (MQ=255)                                                                                                                         aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTc  <  1:1030058/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                    CATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTC  >  NC_000913/3920823‑3921081

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