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%	282.3 / NA	85	*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 (48/37);  total (48/37)

TGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920825‑3921082                                                                                                                              |                                                                                                                                     tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  2:1039365/1‑139 (MQ=255) tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  1:392666/139‑1 (MQ=255) tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  2:392666/1‑139 (MQ=255) tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  1:1865976/1‑139 (MQ=255) tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  1:155289/139‑1 (MQ=255) tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATCGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  1:1039365/139‑1 (MQ=255)           aacgaacAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               <  2:1865976/139‑1 (MQ=255)            acgaacAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              <  1:825834/139‑1 (MQ=255)                acAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagac                                                                                                          >  2:1529722/1‑139 (MQ=255)                     aCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATTCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccaga                                                                                                     >  2:1246872/1‑139 (MQ=255)                     aCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccaga                                                                                                     >  2:1817913/1‑139 (MQ=255)                          cgcgATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcaga                                                                                                           <  2:455083/128‑1 (MQ=255)                          cgcgATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcaga                                                                                                           >  1:455083/1‑128 (MQ=255)                          cgcgATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTAc                                                                                                >  1:707625/1‑139 (MQ=255)                              aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            >  2:808213/1‑139 (MQ=255)                                        cGAGGATTGTTGCGTTAAGATTCACAGCACAATTCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  1:1246872/139‑1 (MQ=255)                                           ggATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGatcat                                                                                       <  2:447053/131‑1 (MQ=255)                                           ggATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGatcat                                                                                       >  1:447053/1‑131 (MQ=255)                                                 ttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGTGAACATCACGTAcagacc                                                                                                         <  2:1382206/107‑1 (MQ=255)                                                 ttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGTGAACATCACGTAcagacc                                                                                                         >  1:1382206/1‑107 (MQ=255)                                                  tgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                        <  1:808213/139‑1 (MQ=255)                                                      ttAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  2:213244/86‑1 (MQ=255)                                                      ttAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  1:213244/1‑86 (MQ=255)                                                      ttAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGa                                                                    >  1:1614102/1‑139 (MQ=255)                                                              cacaGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                    >  1:301231/1‑99 (MQ=255)                                                              cacaGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                    <  2:301231/99‑1 (MQ=255)                                                                        tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  1:403777/1‑84 (MQ=255)                                                                        tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  2:403777/84‑1 (MQ=255)                                                                        tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             <  2:139550/128‑1 (MQ=255)                                                                        tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             >  1:139550/1‑128 (MQ=255)                                                                        tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                  >  2:1765444/1‑139 (MQ=255)                                                                             ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                >  2:1559407/1‑56 (MQ=255)                                                                             ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                <  1:1559407/56‑1 (MQ=255)                                                                             ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCa                                                                                                                   <  1:647059/69‑1 (MQ=255)                                                                             ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCa                                                                                                                   >  2:647059/1‑69 (MQ=255)                                                                              tATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTa                                            >  1:414443/1‑139 (MQ=255)                                                                                  tAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgca                                        >  2:1818684/1‑139 (MQ=255)                                                                                     ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  1:1232992/55‑1 (MQ=255)                                                                                     ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  2:1232992/1‑55 (MQ=255)                                                                                     ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                    <  2:253277/60‑1 (MQ=255)                                                                                     ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                    >  1:253277/1‑60 (MQ=255)                                                                                      tAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGAcc                                                                  >  2:266354/1‑109 (MQ=255)                                                                                      tAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGAcc                                                                  <  1:266354/109‑1 (MQ=255)                                                                                        aCGTTCTGATATTGCTCTTTAAATAAATGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                              <  1:1202446/95‑1 (MQ=255)                                                                                        aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcaga                                                                                                           <  2:512719/66‑1 (MQ=255)                                                                                        aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcaga                                                                                                           >  1:512719/1‑66 (MQ=255)                                                                                        aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                              >  2:1202446/1‑95 (MQ=255)                                                                                          gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGa                                >  1:993666/1‑139 (MQ=255)                                                                                          gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCACCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGa                                >  2:1917029/1‑139 (MQ=255)                                                                                           ttCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGaa                               <  1:1817913/139‑1 (MQ=255)                                                                                             cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatc                             >  1:475249/1‑139 (MQ=255)                                                                                                 tatTGCTCTTTAAATAAAAGCAACGCTTATTACGCGGCAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcagat                         <  2:1614102/139‑1 (MQ=255)                                                                                                 tatTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcagat                         <  2:993666/139‑1 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               >  2:47114/1‑48 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               <  1:47114/48‑1 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              <  1:1057943/49‑1 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              >  2:1057943/1‑49 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  >  1:1892269/1‑77 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  2:1892269/77‑1 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt                    >  2:367931/1‑139 (MQ=255)                                                                                                      ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCCATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt                    >  1:372625/1‑139 (MQ=255)                                                                                                          ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                              <  2:1656222/77‑1 (MQ=255)                                                                                                          ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                              >  1:1656222/1‑77 (MQ=255)                                                                                                          ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGAc                >  2:1133115/1‑139 (MQ=255)                                                                                                           taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  2:1678896/49‑1 (MQ=255)                                                                                                           taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  1:1678896/1‑49 (MQ=255)                                                                                                           taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            <  1:1782231/62‑1 (MQ=255)                                                                                                           taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            >  2:1782231/1‑62 (MQ=255)                                                                                                           taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGAcg               >  1:922330/1‑139 (MQ=255)                                                                                                           taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGAcg               >  1:82285/1‑139 (MQ=255)                                                                                                              ataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACgcgc            >  2:1545746/1‑139 (MQ=255)                                                                                                                 aaaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACg                                                           >  2:1896145/1‑89 (MQ=255)                                                                                                                 aaaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACg                                                           <  1:1896145/89‑1 (MQ=255)                                                                                                                  aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg        <  1:1529722/139‑1 (MQ=255)                                                                                                                  aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg        <  2:82285/139‑1 (MQ=255)                                                                                                                  aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg        <  1:367931/139‑1 (MQ=255)                                                                                                                   aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcgc       >  1:858634/1‑139 (MQ=255)                                                                                                                   aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcgc       >  2:1683860/1‑139 (MQ=255)                                                                                                                       aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacca                                                                      >  1:538370/1‑72 (MQ=255)                                                                                                                       aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacca                                                                      <  2:538370/72‑1 (MQ=255)                                                                                                                       aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    >  2:1279561/1‑90 (MQ=255)                                                                                                                       aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    <  1:1279561/90‑1 (MQ=255)                                                                                                                        cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTcc  >  2:1739057/1‑139 (MQ=255)                                                                                                                          cTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATaaa                                                      >  1:652150/1‑85 (MQ=255)                                                                                                                          cTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATaaa                                                      <  2:652150/85‑1 (MQ=255)                                                                                                                              |                                                                                                                                     TGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920825‑3921082

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