BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,226,832	T→G	21.9%	F78C (TTT→TGT)	ycgJ →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,226,832	0	T	G	21.9%	213.1 / 67.8	115	F78C (TTT→TGT)	ycgJ	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (51/38); new base G (18/7); total (69/46)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 119 total aligned reads displayed.

TGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAAT  >  NC_000913/1226733‑1226929

tgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCTGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAAtt                                                                                                    >  1:207484/1‑100 (MQ=255)

tgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAAtt                                                                                                    >  2:2173616/1‑100 (MQ=255)

 gAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAAttt                                                                                                   <  2:2529462/100‑1 (MQ=255)

  aaCTCAGATTGCATTTCCCCAGAGTTAACAAGGAAATACCTCGGCGAAAAAGCCGCTGAAAACTTACAATCACTACAAGGCTACGATCCCAGCGAATTTa                                                                                                  <  1:1512191/100‑1 (MQ=255)

  aaCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTa                                                                                                  <  2:1664624/100‑1 (MQ=255)

  aaCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTa                                                                                                  >  1:1530981/1‑100 (MQ=255)

  aaCTAATATGGCATTGCCCCAGGGTTAACAAGGAAATATCTCGGCGAAAAAGACGCTGAAAACTTACAATCATTACAAGGCTACGATCGCAGCGAAtt                                                                                                    <  1:1988074/98‑1 (MQ=255)

   aCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTa                                                                                                  <  1:1754554/99‑1 (MQ=255)

   aCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTAc                                                                                                 >  2:1123876/1‑100 (MQ=255)

    cTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACa                                                                                                >  2:200170/1‑100 (MQ=255)

        gATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTc                                                                                             >  1:712745/1‑99 (MQ=255)

        gATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCg                                                                                            <  1:1297339/100‑1 (MQ=255)

        gATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCg                                                                                            >  1:1848376/1‑100 (MQ=255)

         aTGGCATTTCCGCAGAGTTAACAAGGAAATATCTAGGTGAAAAGCCCGCTGAAAACCTAGAATCATTACCAGGGTACGATCCCAGCGAATTTACATTCg                                                                                            >  1:2454863/1‑99 (MQ=255)

           ggCAGTTCCCCAGAGTGAACAAGGACATATCTCGGCGAAAAAGCGGCTGGAATATTATAGGCATGACAAGGCTACGATCCCAGCGCGTATACAt                                                                                               <  1:1199487/94‑1 (MQ=255)

            gCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtgg                                                                                        <  1:1948022/100‑3 (MQ=255)

            gCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtgg                                                                                        <  2:823668/100‑3 (MQ=255)

            gCATTACCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtgg                                                                                        <  1:916712/100‑3 (MQ=255)

             cATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTTCAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAAc                                                                                       >  1:1808108/1‑100 (MQ=255)

             cATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtggg                                                                                       <  2:1962069/100‑4 (MQ=255)

             cATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtggg                                                                                       <  2:2114324/100‑4 (MQ=255)

               tttCACCGGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTATAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtgg                                                                                        <  1:177443/97‑3 (MQ=255)

                   cccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTa                                                                                 >  2:1634179/1‑100 (MQ=255)

                   cccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTa                                                                                 >  1:511358/1‑100 (MQ=255)

                   cccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTa                                                                                 >  1:2080311/1‑100 (MQ=255)

                   cccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTa                                                                                 >  2:543610/1‑100 (MQ=255)

                      agagGTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCAGCACAAGGCTACGAACCCAGTGAATGTAAATTCGCtggg                                                                                       <  1:1807017/91‑4 (MQ=255)

                         gTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGc                                                                           <  2:1998220/100‑1 (MQ=255)

                         gTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGc                                                                           >  2:644693/1‑100 (MQ=255)

                           tAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGCCTACGATCCCAGCGAATTTACATTCGCTAACGGCGTATTCTGCGa                                                                         >  1:2006045/1‑100 (MQ=255)

                              cAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGt                                                                      >  1:1212571/1‑100 (MQ=255)

                              cAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGt                                                                      <  2:658113/100‑1 (MQ=255)

                               aaGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGtt                                                                     >  2:1525284/1‑100 (MQ=255)

                                        tctcGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑aaaaa                                                           >  1:360355/1‑100 (MQ=255)

                                        tctcGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑aaaaa                                                           >  1:1657347/1‑100 (MQ=255)

                                         ctcGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAAt                                                          >  2:144265/1‑100 (MQ=255)

                                           cGGCGAAAAAGCCGCTGAAAACATACAATCATTACAAGGCTACGATCCCAGCGAATTTAAATTCGCTAACGGAGTATTTTGCGATGTTAAAA‑ACACATTa                                                        >  1:2050781/1‑100 (MQ=255)

                                            ggCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTAt                                                       >  1:708766/1‑100 (MQ=255)

                                            ggCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTAt                                                       <  2:1477947/100‑1 (MQ=255)

                                             gCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATg                                                      >  1:2142466/1‑100 (MQ=255)

                                                aaaaaGCCGCTGAAAACTTACAATCATTATAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGTTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCg                                                   <  1:1561361/100‑1 (MQ=255)

                                                aaaaaGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCg                                                   <  1:337421/100‑1 (MQ=255)

                                                 aaaaGCCGCCGAACACTTACAATCATTACAAGGCTACGATCCCAGTGGATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGt                                                  <  1:1784830/100‑1 (MQ=255)

                                                    aGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGtgat                                               <  2:360357/100‑1 (MQ=255)

                                                    aGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGtgat                                               <  2:1657356/100‑1 (MQ=255)

                                                    aGCCGCAGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGtgat                                               <  2:1086634/100‑1 (MQ=255)

                                                     gCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGtgatg                                              >  1:834343/1‑100 (MQ=255)

                                                     gCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGtgatg                                              >  1:834260/1‑100 (MQ=255)

                                                     gCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACAGTGTATTTTGCGATGTTAAAGAAAAAATTATGTCGtgat                                               <  1:1252651/100‑1 (MQ=255)

                                                        gCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATc                                           <  1:1844644/100‑1 (MQ=255)

                                                           gAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGtt                                        >  1:2126170/1‑100 (MQ=255)

                                                            aaaaCTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTa                                       <  2:278764/100‑1 (MQ=255)

                                                             aaaCTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTAt                                      <  2:2126184/100‑1 (MQ=255)

                                                             aaaCTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTAt                                      <  2:2445107/100‑1 (MQ=255)

                                                             aaaCTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTAt                                      >  2:1290408/1‑100 (MQ=255)

                                                                cTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTAtttt                                   <  1:1669683/100‑1 (MQ=255)

                                                                cTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTActtt                                   >  1:484977/1‑100 (MQ=255)

                                                                 ttACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTAtttt                                   >  1:1557990/1‑99 (MQ=255)

                                                                  tACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGCTCGTTATTTTgg                                 >  2:33707/1‑100 (MQ=255)

                                                                  tACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTgg                                 >  2:1923302/1‑100 (MQ=255)

                                                                      aTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGtgtg                             >  1:2549309/1‑100 (MQ=255)

                                                                      aTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGtgtg                             >  2:233117/1‑100 (MQ=255)

                                                                      aTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTCGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGtgtg                             >  2:233104/1‑100 (MQ=255)

                                                                       tCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AACAACTATGTCGTGATGATccgatcggaa                                  >  1:2464333/1‑85 (MQ=255)

                                                                         aTTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGAt                          <  1:2579722/100‑1 (MQ=255)

                                                                         aTTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGAt                          >  2:51642/1‑100 (MQ=255)

                                                                         aTTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGAt                          >  1:1056863/1‑100 (MQ=255)

                                                                          ttACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATg                         >  2:1488466/1‑100 (MQ=255)

                                                                           tACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATgg                        <  2:1199629/100‑1 (MQ=255)

                                                                            aCAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGt                       <  2:854272/100‑1 (MQ=255)

                                                                             cAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTa                      <  2:1535734/100‑1 (MQ=255)

                                                                               aGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAg                    >  1:1963265/1‑100 (MQ=255)

                                                                                 gCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCg                  <  1:1072344/100‑1 (MQ=255)

                                                                                   tccGATCCCAGCGAATTTAACTTCGCCAACGGTGTATTTTGCGCTGTTAAAG‑AAAAATTATGTCGTGATGATCGTTGTGTTGGTGTGGATGGTAAGCg                  >  1:1772313/3‑98 (MQ=255)

                                                                                   tacaccgttaGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGa                <  2:7256/91‑1 (MQ=255)

                                                                                   tACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGa                >  1:1774979/1‑100 (MQ=255)

                                                                                     gaccggtaGCGGATGTAAATTGGCTAACGGTGTATTTTGCGGTGTTAAAG‑AAAAATTATGTCGGGATGATCGTGATTTTGGTGTGGATGGTAAGCg                  >  1:513255/8‑96 (MQ=255)

                                                                                     caccgttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGTTCGTTATTTTGGTGTGGATGGTAAgaggag               >  1:1618030/8‑99 (MQ=38)

                                                                                     caccgttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGt              >  2:2048736/8‑100 (MQ=255)

                                                                                      accgttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTg             >  1:1977388/7‑100 (MQ=255)

                                                                                       ccgttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTgg            >  2:2354678/6‑100 (MQ=255)

                                                                                       ccgttaGCGAATGGAAATTCGCTAACGGTGTATTTTGCAATGTTAAAG‑AAAAATTATGTCGTGATTATCGTTATTTTGGTGTGGATGGTAAGCGGACTgg            >  1:1596047/6‑100 (MQ=255)

                                                                                        ngttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGa           >  2:1411579/5‑100 (MQ=255)

                                                                                         nttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaa          >  1:1338898/4‑100 (MQ=255)

                                                                                         gttcgcaaaTGTAAATTCGATAACGGTGTATTTTGCGATGTTAAAG‑AATAATTATGTAGTGATGATCGCTATTTTGGTGAGGATGGTAAGCGGAGTGGaa          >  1:2474903/8‑100 (MQ=255)

                                                                                         gttaGCGAATGTAAATTCGCTCACGGTGTACTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaa          >  1:633459/4‑100 (MQ=255)

                                                                                         gttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGCTGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGGGGaa          >  2:807153/4‑100 (MQ=255)

                                                                                         gttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaa          >  2:2053405/4‑100 (MQ=255)

                                                                                          ttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaa         >  2:1735721/3‑100 (MQ=255)

                                                                                          ccAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaa         >  2:741473/1‑100 (MQ=255)

                                                                                           taGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaaa        >  2:2509385/2‑100 (MQ=255)

                                                                                           taGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaaa        >  1:885332/2‑100 (MQ=255)

                                                                                           taGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaaa        >  2:796042/2‑100 (MQ=255)

                                                                                           cAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaaa        <  1:449927/100‑1 (MQ=255)

                                                                                            aGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaaaa       >  1:1122921/1‑100 (MQ=255)

                                                                                             gCGAATTTACATTCGGTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaat      <  2:1801605/100‑1 (MQ=255)

                                                                                             gCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaat      <  1:1291316/100‑1 (MQ=255)

                                                                                              cGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatc     <  2:2289059/100‑1 (MQ=255)

                                                                                                 aTTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatcaat  >  1:1845670/1‑100 (MQ=255)

                                                                                                 aTTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatcaat  <  2:207486/100‑1 (MQ=255)

TGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAG‑AAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAAT  >  NC_000913/1226733‑1226929

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