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.5 / 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)

TTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920821‑3921082                                                                                                                                  |                                                                                                                                     ttCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACa                                                                                                                             >  2:363113/1‑139 (MQ=255)     tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  1:1148631/1‑139 (MQ=255)     tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  2:216785/1‑139 (MQ=255)     tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  2:1148631/139‑1 (MQ=255)     tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  2:723874/1‑139 (MQ=255)          aacagaacGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                    <  1:723874/139‑1 (MQ=255)          aacagaacGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                    >  2:849347/1‑139 (MQ=255)                    acAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                >  1:645727/1‑117 (MQ=255)                    acAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                <  2:645727/117‑1 (MQ=255)                     cAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  2:153551/139‑1 (MQ=255)                      aGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaacc                                                                                                        >  2:597828/1‑139 (MQ=255)                       ggACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATGGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccca                                                                                                       >  1:952736/1‑139 (MQ=255)                         aCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccaga                                                                                                     >  2:680605/1‑139 (MQ=255)                                  aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            <  1:849347/139‑1 (MQ=255)                                  aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            <  2:952736/139‑1 (MQ=255)                                       ctggccGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGatcat                                                                                       <  1:363113/139‑1 (MQ=255)                                            cGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  2:738748/139‑1 (MQ=255)                                                    gttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                >  2:714443/1‑85 (MQ=255)                                                    gttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                <  1:714443/85‑1 (MQ=255)                                                        cGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                <  2:185789/81‑1 (MQ=255)                                                        cGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                >  1:185789/1‑81 (MQ=255)                                                             aGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAca                                                                                                             <  2:1214940/95‑1 (MQ=255)                                                             aGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAca                                                                                                             >  1:1214940/1‑95 (MQ=255)                                                                  cacaGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccca                                                                                                       >  2:1362165/1‑96 (MQ=255)                                                                  cacaGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccca                                                                                                       <  1:1362165/96‑1 (MQ=255)                                                                       cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       <  2:139054/123‑1 (MQ=255)                                                                       cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       >  1:528266/1‑123 (MQ=255)                                                                       cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       <  2:528266/123‑1 (MQ=255)                                                                       cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       >  1:139054/1‑123 (MQ=255)                                                                       cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATaa                                                       >  2:1303257/1‑139 (MQ=255)                                                                         caATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAg                                                     >  2:1447613/1‑139 (MQ=255)                                                                         caATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAg                                                     >  2:385674/1‑139 (MQ=255)                                                                          aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    >  2:1218106/1‑139 (MQ=255)                                                                          aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    <  1:385674/139‑1 (MQ=255)                                                                          aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    <  1:216785/139‑1 (MQ=255)                                                                            tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                  <  1:680605/139‑1 (MQ=255)                                                                            tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                  >  1:465038/1‑139 (MQ=255)                                                                            tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                  >  2:261256/1‑139 (MQ=255)                                                                                 ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              <  2:1445329/74‑1 (MQ=255)                                                                                 ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              >  1:1445329/1‑74 (MQ=255)                                                                                 ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATAACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGt                                             >  1:941783/1‑139 (MQ=255)                                                                                    tttAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACg                                          <  2:465038/139‑1 (MQ=255)                                                                                       aGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGaa                                                   >  1:41021/1‑127 (MQ=255)                                                                                       aGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGaa                                                   <  2:41021/127‑1 (MQ=255)                                                                                        gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATa                                                        >  1:1405836/1‑121 (MQ=255)                                                                                        gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATa                                                        <  2:1405836/121‑1 (MQ=255)                                                                                         ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATGAAGAACTGAGTACgcagca                                     <  1:1303257/139‑1 (MQ=255)                                                                                         ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagca                                     >  2:1508883/1‑139 (MQ=255)                                                                                          tAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGcagcag                                    <  1:261256/139‑1 (MQ=255)                                                                                           aaCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                                   <  2:941783/139‑1 (MQ=255)                                                                                           aaCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                                   <  1:953181/139‑1 (MQ=255)                                                                                              gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                        >  2:1503094/1‑99 (MQ=255)                                                                                              gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                        <  1:1503094/99‑1 (MQ=255)                                                                                                  tGATATTGCTCTTTAAATAAAAGCAATGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               <  2:468378/56‑1 (MQ=255)                                                                                                  tGATATTGCTCTTTAAATAAAAGCAATGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               >  1:468378/1‑56 (MQ=255)                                                                                                  tGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatca                            >  2:1151471/1‑139 (MQ=255)                                                                                                        tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcagatcag                      >  1:396091/1‑139 (MQ=255)                                                                                                        tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcagatcag                      >  2:378968/1‑139 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               <  1:999515/49‑1 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  1:1152626/55‑1 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  2:1152626/1‑55 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              <  2:1245197/50‑1 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              >  1:1245197/1‑50 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                             <  2:1592984/83‑1 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                             >  1:1592984/1‑83 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATc                                                                           >  2:911308/1‑85 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATc                                                                           <  1:911308/85‑1 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAgg                     >  1:492243/1‑139 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTAAGCGACAGCGAACATCACGTa                                                                                                               >  2:999515/1‑49 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  >  2:1416880/1‑77 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  1:1416880/77‑1 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  >  2:608421/1‑77 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  1:608421/77‑1 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       <  2:489194/88‑1 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       >  1:489194/1‑88 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    >  1:745523/1‑107 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    <  2:745523/107‑1 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcag                           <  1:645231/132‑1 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcag                           >  2:645231/1‑132 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt                    >  1:675242/1‑139 (MQ=255)                                                                                                              ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                                   <  2:1315911/120‑1 (MQ=255)                                                                                                              ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                                   >  1:1315911/1‑120 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             <  1:1145653/93‑1 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             >  2:1145653/1‑93 (MQ=255)                                                                                                                      aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg        >  2:769723/1‑139 (MQ=255)                                                                                                                       aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  >  1:233903/1‑64 (MQ=255)                                                                                                                       aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  2:233903/64‑1 (MQ=255)                                                                                                                           aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTc   <  2:396091/139‑1 (MQ=255)                                                                                                                            cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            >  1:327556/1‑49 (MQ=255)                                                                                                                            cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            <  2:327556/49‑1 (MQ=255)                                                                                                                            cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTcc  >  1:1265713/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                     TTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920821‑3921082

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