BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,649,933	C→A	P47Q (CCG→CAG)	ydfD →	Qin prophage; lysis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,649,933	0	C	A	100.0%	560.7 / NA	155	P47Q (CCG→CAG)	ydfD	Qin prophage; lysis protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (62/93); total (62/93)

Only 100 of 159 total aligned reads displayed.

AATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGT  >  NC_000913/1649797‑1650069

aaTTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGt                                                                                                                                       <  1:908193/140‑1 (MQ=255)

 aTTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTc                                                                                                                                      <  1:127256/140‑1 (MQ=255)

   tCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGa                                                                                                                                    <  2:567610/140‑1 (MQ=255)

   tCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGa                                                                                                                                    <  1:2051630/140‑1 (MQ=255)

   tCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGa                                                                                                                                    <  2:1163371/140‑1 (MQ=255)

    cAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGCCGAt                                                                                                                                   >  2:798254/1‑140 (MQ=255)

      gCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATaa                                                                                                                                 <  1:278996/140‑1 (MQ=255)

      gCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATaa                                                                                                                                 >  2:1506603/1‑140 (MQ=255)

      gCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATaa                                                                                                                                 <  2:1916820/140‑1 (MQ=255)

        aTTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAg                                                                                                                               >  2:1787242/1‑140 (MQ=255)

          ttGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCGACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGtt                                                                                                                             >  1:2229673/1‑140 (MQ=255)

          ttGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGtt                                                                                                                             >  1:1427767/1‑140 (MQ=255)

           tgtgCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTa                                                                                                                            <  1:435431/140‑1 (MQ=255)

           tgtgCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTa                                                                                                                            <  2:156897/140‑1 (MQ=255)

           tgtgCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTa                                                                                                                            <  1:435354/140‑1 (MQ=255)

           tgtgCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTa                                                                                                                            <  1:979514/140‑1 (MQ=255)

             tgCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGTAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTAtt                                                                                                                          <  2:1101432/140‑1 (MQ=255)

             tgCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTAtt                                                                                                                          <  1:1904364/140‑1 (MQ=255)

               cTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCa                                                                                                                        >  2:2020611/1‑140 (MQ=255)

                ttGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCAc                                                                                                                       <  1:1506603/140‑1 (MQ=255)

                ttGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCAc                                                                                                                       <  2:538102/140‑1 (MQ=255)

                ttGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCAc                                                                                                                       <  1:1911275/140‑1 (MQ=255)

                 tgctCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTAtt                                                                                                                          >  2:874877/4‑136 (MQ=255)

                 tGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTAtt                                                                                                                          <  1:874877/136‑1 (MQ=255)

                 tGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCAcc                                                                                                                      >  2:245850/1‑140 (MQ=255)

                  gTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGa                                                                                                                                    <  1:2098398/125‑1 (MQ=255)

                  gTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGa                                                                                                                                    >  2:2098398/1‑125 (MQ=255)

                   ttCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCAc                                                                                                                       >  1:599009/1‑137 (MQ=255)

                    tCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAgg                                                                                                                   >  1:1564791/1‑140 (MQ=255)

                         cAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGataat                                                                                                              >  1:1274629/1‑140 (MQ=255)

                            tttttCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATc                                                                                                           <  1:1998936/140‑1 (MQ=255)

                              tttCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGa                                                                                                         <  2:1713204/140‑1 (MQ=255)

                              tttCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGa                                                                                                         <  1:1187983/140‑1 (MQ=255)

                                tCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGaaa                                                                                                       <  1:879430/140‑1 (MQ=255)

                                tCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGaaa                                                                                                       >  1:883781/1‑140 (MQ=255)

                                      ttCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCgg                                                                                                 <  2:165343/140‑1 (MQ=255)

                                         cGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAg                                                                                              <  1:712398/140‑1 (MQ=255)

                                         cGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAg                                                                                              <  2:1206359/140‑1 (MQ=255)

                                         cGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAg                                                                                              <  2:1207956/140‑1 (MQ=255)

                                         cGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAg                                                                                              <  2:1225129/140‑1 (MQ=255)

                                         cGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAg                                                                                              <  1:1002531/140‑1 (MQ=255)

                                           gagaGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACTAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGt                                                                                            <  2:1412163/140‑1 (MQ=255)

                                           gagaGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGt                                                                                            <  1:883412/140‑1 (MQ=255)

                                                   gTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTaaaa                                                                                    <  2:903828/140‑1 (MQ=255)

                                                   gTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTaaaa                                                                                    <  1:1864675/140‑1 (MQ=255)

                                                       atatTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGt                                                                                <  2:24114/140‑1 (MQ=255)

                                                        tatTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGtt                                                                               <  1:130095/140‑1 (MQ=255)

                                                        tatTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGtt                                                                               <  2:1564791/140‑1 (MQ=255)

                                                         atTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAg                                                                                              <  2:1507540/124‑1 (MQ=255)

                                                         atTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTc                                                                              >  1:809409/1‑140 (MQ=255)

                                                          ttGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTcc                                                                             >  1:1599075/1‑140 (MQ=255)

                                                             cagtcagtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGta                                                                          >  1:1895707/1‑140 (MQ=255)

                                                             cagtcagtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGta                                                                          >  1:1992882/1‑140 (MQ=255)

                                                               gtcagtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGtaat                                                                        <  2:2016846/140‑1 (MQ=255)

                                                                tcagtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGtaata                                                                       >  2:772999/1‑140 (MQ=255)

                                                                  agtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTaataaa                                                                     >  2:164909/1‑140 (MQ=255)

                                                                      ttCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCCCCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATc                                                                 <  2:2052229/140‑1 (MQ=255)

                                                                         acaGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCgg                                                              >  1:1764911/1‑140 (MQ=255)

                                                                          caGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGt                                                             <  1:1601451/140‑1 (MQ=255)

                                                                              aCAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCa                                                         >  2:1714578/1‑140 (MQ=255)

                                                                                aaTGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCAtt                                                       >  1:2043590/1‑140 (MQ=255)

                                                                                        aGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATg                                               <  1:1494231/140‑1 (MQ=255)

                                                                                          tgtATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCAt                                                        <  1:1921962/129‑1 (MQ=255)

                                                                                           gtATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCa                                            >  2:1274624/1‑140 (MQ=255)

                                                                                            tATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGttt                                                           >  2:1966818/1‑124 (MQ=255)

                                                                                            tATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAg                                           >  1:1045738/1‑140 (MQ=255)

                                                                                            tATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAg                                           >  1:1466323/1‑140 (MQ=255)

                                                                                            tATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAg                                           >  2:2126097/1‑140 (MQ=255)

                                                                                            tATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAg                                           >  1:2020894/1‑140 (MQ=255)

                                                                                            tATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAg                                           >  1:2084931/1‑140 (MQ=255)

                                                                                              tGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCa                                         >  1:2225273/1‑140 (MQ=255)

                                                                                                 cTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATg                                      >  1:216506/1‑140 (MQ=255)

                                                                                                 cTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATg                                      >  2:261262/1‑140 (MQ=255)

                                                                                                   gcagcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGc                                    >  2:1989765/1‑140 (MQ=255)

                                                                                                   gcagcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGc                                    >  2:1591116/1‑140 (MQ=255)

                                                                                                   gcagcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGc                                    <  2:433899/140‑1 (MQ=255)

                                                                                                    cagcgACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCa                                   <  2:2229673/140‑1 (MQ=255)

                                                                                                    cagcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCa                                   >  2:186300/1‑140 (MQ=255)

                                                                                                    cagcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCa                                   <  2:442821/140‑1 (MQ=255)

                                                                                                    cagcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCa                                   >  1:259287/1‑140 (MQ=255)

                                                                                                    cagcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGCTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCa                                   <  2:1923646/140‑1 (MQ=255)

                                                                                                     agcaACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAg                                  <  1:1591116/140‑1 (MQ=255)

                                                                                                       caACCGAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAggg                                >  1:334939/1‑140 (MQ=255)

                                                                                                            gAACAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCa                                   <  1:636352/132‑1 (MQ=255)

                                                                                                             aaCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGt                          <  1:2132151/140‑1 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         <  1:2126097/140‑1 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         >  2:611574/1‑140 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         <  2:1274629/140‑1 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         <  2:1222715/140‑1 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         <  1:1519195/140‑1 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         <  1:1535479/140‑1 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         >  1:724749/1‑140 (MQ=255)

                                                                                                              aCAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGtt                         <  1:1714578/140‑1 (MQ=255)

                                                                                                               cAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTc                        <  2:2050844/140‑1 (MQ=255)

                                                                                                               cAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTc                        <  2:64770/140‑1 (MQ=255)

                                                                                                               cAGAAAATTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTc                        <  2:1143222/140‑1 (MQ=255)

                                                                                                                    aaTTCCCGGTAACTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCt                   >  1:1211814/1‑140 (MQ=255)

                                                                                                                                cTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCAt       <  2:1203393/140‑1 (MQ=255)

                                                                                                                                cTGTTACCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCAt       <  2:2084931/140‑1 (MQ=255)

                                                                                                                                     aCCAGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGt  <  2:334939/140‑1 (MQ=255)

AATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGT  >  NC_000913/1649797‑1650069

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