BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,702,500	C→G	5.2%	H90D (CAC→GAC)	add →	adenosine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,702,500	0	C	G	5.2%	413.9 / 10.1	129	H90D (CAC→GAC)	add	adenosine deaminase
Reads supporting (aligned to +/- strand): ref base C (67/55); new base G (5/2); total (72/57)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 130 total aligned reads displayed.

GGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATT  >  NC_000913/1702352‑1702649

ggAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCa                                                                                                                                                      >  2:697411/1‑150 (MQ=255)

     cacTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGCCCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGcc                                                                                                                                                 >  1:876702/1‑150 (MQ=255)

     cacTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGACACGGcc                                                                                                                                                 >  1:1010044/1‑150 (MQ=255)

        tGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGc                                                                                                                                              >  2:1270727/1‑150 (MQ=255)

         gATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCa                                                                                                                                             <  1:1712716/150‑1 (MQ=255)

              cccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATg                                                                                                                                        >  1:395951/1‑150 (MQ=255)

               ccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTATCTGACCAAACTTGACTGGGGCGTTAAAGGTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGt                                                                                                                                       <  1:1250936/150‑1 (MQ=255)

                  cGTTCAGGTCTTTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCCCGGCCTGCACTATGTCGa                                                                                                                                    <  1:1848827/150‑1 (MQ=255)

                  cGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGa                                                                                                                                    <  2:737528/150‑1 (MQ=255)

                   gTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAg                                                                                                                                   <  1:1621282/150‑1 (MQ=255)

                    ttCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGc                                                                                                                                  <  1:1252632/150‑1 (MQ=255)

                        ggTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCg                                                                                                                              <  2:596343/150‑1 (MQ=255)

                        ggTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGAAGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCg                                                                                                                              >  1:727732/1‑150 (MQ=255)

                          tCAATGCCAACGAACCCGATCTGGTCAGCTTTCTGCCCAACCTTGCCGGCGGGGTTAAAGTTCTCGCCTCTCTTGCTGCCGGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCCCGGCCTGCACTATGTCGAGCTGCGtt                                                                                                                            <  2:1184684/150‑1 (MQ=255)

                            aTTGCCAACGAACCCTATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCGCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGtttt                                                                                                                          <  2:22815/150‑1 (MQ=255)

                               gCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTcc                                                                                                                       >  2:1023628/1‑149 (MQ=255)

                               gCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTATCa                                                                                                                       >  2:891177/1‑150 (MQ=255)

                                ccAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCAc                                                                                                                      >  2:914615/1‑150 (MQ=255)

                                ccAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCAc                                                                                                                      >  1:538508/1‑150 (MQ=255)

                                 cAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCAcc                                                                                                                     >  1:1563437/1‑150 (MQ=255)

                                  aaccaaccCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGGTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCGGCGGTTTTCACCa                                                                                                                    >  2:669838/1‑150 (MQ=255)

                                     gAAACCGATCTGGAGAGCTTTCTGCCCAAACTTGACTGGTGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCACGTCACTGCCTGGACTATGTCGCGCTGAGTTTTTCACCAGGc                                                                                                                 >  2:1835240/1‑150 (MQ=255)

                                      aaCCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCt                                                                                                                >  2:215479/1‑150 (MQ=255)

                                           gATCTGGTGAGCTTTCTAACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatg                                                                                                           >  1:412445/1‑150 (MQ=255)

                                              cTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCGCTTGATGCCTGGCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggca                                                                                                        <  2:185297/150‑1 (MQ=255)

                                                ggTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGACGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaat                                                                                                      >  2:769431/1‑150 (MQ=255)

                                                   gAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggc                                                                                                   <  1:1415550/150‑1 (MQ=255)

                                                   gAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCATGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggc                                                                                                   <  2:1303009/150‑1 (MQ=255)

                                                    aGCTTTCTGACCAAACGTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggca                                                                                                  <  1:738638/150‑1 (MQ=255)

                                                      cTTTCTGACCAAACTTGACTGGGGCGTAAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCGTGCAATATGTCGAGATTCGTTTTTCACCACGATACATGGCAATGTcaca                                                                                                >  1:106659/1‑150 (MQ=255)

                                                             aCCAAACTTGACTGGGGCGTTAAAGTGCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTg                                                                                         <  1:1134885/150‑1 (MQ=255)

                                                              ccAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGc                                                                                        >  2:1558567/1‑150 (MQ=255)

                                                               cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGcc                                                                                       <  1:951514/150‑1 (MQ=255)

                                                               cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGcc                                                                                       <  2:804057/150‑1 (MQ=255)

                                                                   cTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTa                                                                                   >  2:212882/1‑150 (MQ=255)

                                                                   cTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTa                                                                                   >  1:58121/1‑150 (MQ=255)

                                                                   cTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCGACATGGCAATGGCACATCAGCTGCCTGTa                                                                                   <  1:121805/150‑1 (MQ=255)

                                                                     tGACTGCGGCGTTACCGTTCTCGCCTCGCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGCTGCACCCCGTCACGGCCTGCACTATGTCGACCTGCGTTTTTCACCAGGCTACTTGGCAATGGCACATCAGCTGCCTGTAGc                                                                                 <  2:148141/150‑1 (MQ=255)

                                                                          ggggCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGgtg                                                                            <  1:805815/150‑1 (MQ=255)

                                                                           gggCGTTAATGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGTTGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGgtgt                                                                           <  2:991473/150‑1 (MQ=255)

                                                                           gggCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGgtgt                                                                           >  2:678818/1‑150 (MQ=255)

                                                                           gggCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGgtgt                                                                           >  2:1115699/1‑150 (MQ=255)

                                                                             gCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGtgttg                                                                         >  1:1630716/1‑150 (MQ=255)

                                                                              cGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGtgttgt                                                                        <  2:729567/150‑1 (MQ=255)

                                                                              cGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGtgttgt                                                                        >  2:445697/1‑150 (MQ=255)

                                                                              cGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGtgttgt                                                                        >  1:1585169/1‑150 (MQ=255)

                                                                               gTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGGGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTc                                                                       <  1:507286/150‑1 (MQ=255)

                                                                                tttaaGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGTCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCg                                                                      <  1:603812/147‑1 (MQ=255)

                                                                                ttAAAGTTCTCGCCTCGCTTGAGGCCTGTCGCCGCGTGGCATTTGAAAACATTGACGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCg                                                                      <  1:132130/150‑1 (MQ=255)

                                                                                  aaaGCTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAATACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGaa                                                                    <  1:1596994/150‑1 (MQ=255)

                                                                                      ttCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCgg                                                                <  1:777744/150‑1 (MQ=255)

                                                                                      ttCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCgg                                                                <  1:1705392/150‑1 (MQ=255)

                                                                                       tctcGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGt                                                               >  1:176442/1‑150 (MQ=255)

                                                                                        ctcGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTg                                                              >  1:539872/1‑150 (MQ=255)

                                                                                          cGCCTCTCTTGATGCTTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGAt                                                            >  2:1112399/1‑150 (MQ=255)

                                                                                          cGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGAt                                                            >  1:1353701/1‑150 (MQ=255)

                                                                                           gCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATc                                                           >  1:1068766/1‑150 (MQ=255)

                                                                                            ccTCTCTTGATGCCTGTCGCCGCGTGGCATGTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCg                                                          >  1:737375/1‑150 (MQ=255)

                                                                                             ctctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGa                                                         >  2:1027662/1‑150 (MQ=255)

                                                                                              tctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCCGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGAt                                                        <  1:1110551/150‑1 (MQ=255)

                                                                                                tctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATgg                                                      <  1:1557775/150‑1 (MQ=255)

                                                                                                 ctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGc                                                     >  2:777327/1‑150 (MQ=255)

                                                                                                 ctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGGCGAGCTGCGGTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGc                                                     <  2:901905/150‑1 (MQ=255)

                                                                                                  ttGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCg                                                    <  1:1270727/150‑1 (MQ=255)

                                                                                                     aTGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTAc                                                 <  2:1294363/150‑1 (MQ=255)

                                                                                                      tGCCTGTCGCCGCGTGTCATTTGAAAACCTTGCAGATGCATCCCGTCACTGCCTTCACTATGTCGAGCTGCTTTTTTCACAAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTTTCGAAGCGGTGATCGATGGCGTACg                                                <  1:959631/150‑1 (MQ=255)

                                                                                                          tGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGaa                                            <  1:1792681/150‑1 (MQ=255)

                                                                                                           gTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCTCGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACTTCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAg                                           >  2:71387/1‑150 (MQ=255)

                                                                                                             cgccgcGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGt                                         >  1:1703890/1‑150 (MQ=255)

                                                                                                              gccgcGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGtt                                        >  1:1136315/1‑150 (MQ=255)

                                                                                                                 gcgTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGcc                                     <  2:1709200/150‑1 (MQ=255)

                                                                                                                   gTGGCATTTGAAAACCATGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCGGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGc                                   <  2:1210777/150‑1 (MQ=255)

                                                                                                                   gTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGc                                   <  2:534532/150‑1 (MQ=255)

                                                                                                                     ggCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACAGGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGACGCAc                                 >  1:223767/1‑150 (MQ=255)

                                                                                                                      gCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCAcc                                <  2:1068766/150‑1 (MQ=255)

                                                                                                                       cATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCt                               >  1:1786879/1‑150 (MQ=255)

                                                                                                                        aTTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCtt                              >  1:1353450/1‑150 (MQ=255)

                                                                                                                        aTTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCtt                              >  2:694042/1‑150 (MQ=255)

                                                                                                                           tGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTgg                           >  1:581170/1‑150 (MQ=255)

                                                                                                                             aaaaCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGTCTACATGGCAATGTCACATCAGCTGCCTGTAGCGGTGGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGgtg                         <  1:116890/150‑1 (MQ=255)

                                                                                                                               aaCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGtgtg                       <  2:539872/150‑1 (MQ=255)

                                                                                                                               aaCATTGAAGATGCAGCCCGCCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGtgtg                       <  2:187696/150‑1 (MQ=255)

                                                                                                                                aCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGc                      <  1:975884/150‑1 (MQ=255)

                                                                                                                                    tGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGc                  >  1:331860/1‑150 (MQ=255)

                                                                                                                                    tGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGATACATGGCAATGGCACATAAGCTGCCAGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTGTGGTGTGCAGGc                  >  2:1004800/1‑150 (MQ=255)

                                                                                                                                     gAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCg                 <  2:1742456/150‑1 (MQ=255)

                                                                                                                                     gAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCg                 >  1:789281/1‑150 (MQ=255)

                                                                                                                                      atagtGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGACGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGa                >  1:926651/5‑150 (MQ=255)

                                                                                                                                      aaGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGa                >  2:1725133/1‑150 (MQ=255)

                                                                                                                                        agtGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCAGTGATCGATGGCGTACGTGAAGGTTGCCGCAACTTTGGTGTGCAGGCGAAg              >  2:398483/3‑150 (MQ=255)

                                                                                                                                         ttGCAGCCTGTGACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTGTGGTGTGCAGGCGAAGc             <  1:222729/149‑1 (MQ=255)

                                                                                                                                         aTGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGc             <  1:1626102/150‑1 (MQ=255)

                                                                                                                                          tGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt            >  1:195436/1‑150 (MQ=255)

                                                                                                                                          tGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGAGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt            >  1:1385778/1‑150 (MQ=255)

                                                                                                                                          tGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTCTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt            >  1:1378148/1‑150 (MQ=255)

                                                                                                                                            cAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCCGCTGCCTGTAGCGGGGGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTa          >  1:855716/1‑150 (MQ=255)

                                                                                                                                            cAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTa          >  2:146789/1‑150 (MQ=255)

                                                                                                                                                ccGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCgg      >  1:661054/1‑150 (MQ=255)

                                                                                                                                                   tCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCAt   >  2:976859/1‑150 (MQ=255)

                                                                                                                                                    cACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACATGAAGGTTGCCGCACCGTTGGTGTGCAGGCGAAGCTTATCGGCAtt  >  2:1661232/1‑150 (MQ=255)

GGAAACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATT  >  NC_000913/1702352‑1702649

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